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Charles W. Walker et al. 


Abstract 

The human P53 tumour suppressor protein is inactivated in many cancers and 
is also a major player in apoptotic responses to cellular stress. The P53 protein 
and the two other members of this protein family (p63, P73) are encoded by 
distinct genes and their functions have been extensively documented for 
humans and some other vertebrates. The structure and relative expression 
levels for members of the P53 superfamily have also been reported for most 
major invertebrate taxa. The functions of homologous proteins have been 
investigated for only a few invertebrates (specifically, P53 in flies, nematodes 
and recently a sea anemone). These studies of classical model organisms all 
suggest that the gene family originally evolved to mediate apoptosis of dam¬ 
aged germ cells or to protect germ cells from genotoxic stress. Here, we have 
correlated data from a number of molluscan and other invertebrate sequencing 
projects to provide a framework for understanding P53 signalling pathways in 
marine bivalve cancer and stress biology. These data suggest that (a) the two 
identified P53 and P63/73-like proteins in soft shell clam ( Mya arenaria ), blue 
mussel ( Mytilus edulis) and Northern European squid ( Loligo forbesi) have 
identical core sequences and may be splice variants of a single gene, while 
some molluscs and most other invertebrates have two or more distinct genes 
expressing different P53 family members; (b) transcriptional activation 
domains (TADs) in bivalve P53 and P63/73-like protein sequences are 
67-69% conserved with human P53, while those in ecdysozoan, cnidarian, 
placozoan and choanozoan eukaryotes are <33% conserved; (c) the Mdm2 
binding site in the transcriptional activation domain is 100% conserved in all 
sequenced bivalve P53 proteins (e.g. Mya, Mytilus, Crassostrea and Spisula) 
but is not present in other non-deuterostome invertebrates; (d) an Mdm2 
homologue has been cloned for Mytilus trossulus-, (e) homologues for both 
human P53 upstream regulatory and transcriptional target genes exist in 
molluscan genomes (missing are ARF, CIPi and BH3 only proteins) and (f) 
P53 is demonstrably involved in bivalve haemocyte and germinoma cancers. 
We usually do not know enough about the molecular biology of marine 
invertebrates to address molecular mechanisms that characterize particular 
diseases. Understanding the molecular basis of naturally occurring dis¬ 
eases in marine bivalves is a virtually unexplored aspect of toxicoproteo- 
mics and genomics and related drug discovery. Additionally, increases in 
coastal development and concomitant increases in aquatic pollutants have 
driven interest in developing models appropriate for evaluating potential 
hazardous compounds or conditions found in the aquatic environment. Data 
reviewed in this study are coupled with recent developments in our under¬ 
standing the molecular biology of the marine bivalve P53 superfamily. 
Taken together, they suggest that both structurally and functionally, bivalve 
P53 family proteins are the most highly conserved members of this gene 
superfamily so far identified outside of higher vertebrates and invertebrate 
chordates. Marine bivalves provide some of the most relevant and best 
understood models currently available for experimental studies by biomed¬ 
ical and marine environmental researchers. 
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Introduction 


In humans, the p53 gene family includes three isoforms, p53, p63 and 
p73 that are located on different chromosomes (p53, chromosome 17pl2.1; 
p63, chromosome 3q27-29; p73, chromosome lp36.3). Distinct p53-, 
63- and 73-related genes are also found in other mammalian and non- 
mammalian vertebrates (Le Bras etal, 2003). All identified family members 
produce splice variants that are uniquely expressed in different cell types, both 
spatially and temporally (Murray-Zmijewski et al., 2006; Khoury and 
Bourdon, 2010). While p53-related genes are absent from plant and fungal 
genomes, both p53- and p63/73-like genes have been detected throughout 
the non-vertebrate animal world (Wahl and Carr, 2001; Nedelcu and Tan, 
2007; Pankow and Bamberger, 2007; Rutkowski etal., 2010). Distinct p53 
and p63/73-like family members can be identified in completed genome 
projects for invertebrate deuterostomes, lophotrochozoan, ecdysozoan, cni- 
darian, placozoan and choanozoan eukaryotes, some of which include splice 
variants (Fig. 1.1). Comparative evidence from available eukaryotic genome 
projects suggests that vertebrate p53, p63 and p73 genes form a monophy- 
logenetic group distinct from invertebrate p53 superfamily members 
(Nedelcu and Tan, 2007; Fernandes and Atchley, 2008; Rutkowski etal., 
2010). These studies also indicate that an ancient p63-like gene was ancestral 
to all p53- and p73-related genes in the major lines of animal evolution and 
existed from 500 million to 1 billion years ago (Belyi etal., 2010). This 
ancestral p63 gene developed in the germ line and subsequently gave rise to 
the alternate forms of genes in the p53 superfamily now found in both germ 
and somatic cells of widely separated phylogenetic taxa (Lu and Abrams, 
2006; Nedelcu and Tan, 2007; Belyi et al., 2010; Dotsch et al., 2010; 
Rutkowski etal., 2010). Functional data for several nematodes and flies 
indicate that a p53-related protein maintains genomic stability after geno- 
toxic challenge and prevents germ cell proliferation upon induced expres¬ 
sion (Akdemir et al., 2007; Belyi etal., 2010). Additionally, p63 protects the 
mammalian female germ line during meiotic arrest (Suh et al., 2006). Similar 
results have been obtained for a p63/73-like protein in the sea anemone 
Nematostella vectensis that is involved in selective germ cell death (Akdemir 
etal., 2007; Pankow and Bamberger, 2007). With the exception of the 
sea anemone, genes in the cellular senescence, DNA editing and repair 
pathways governed by the transcriptional activation domain (TAD) are 
not present in these organisms, perhaps because their adult somatic cells 
do not divide by mitosis and their TAD is only 2% ( Drosophila ) and 33% 
(Caenorhabditis ) conserved respectively relative to human p53. 

In the vertebrates, transcriptional target genes for the p53 protein are 
focused in four areas: down-regulation ofp53, cellular senescence, DNA 
editing/repair and apoptosis. These functions are important in both 
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Molluscs 

Non-vertebrate 
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Figure 1.1 Phylogenetic tree of the p53 superfamily. Sequences corresponding to 
the conserved DNA binding domain were used to generate an Unrooted 
Phylogenetic Tree using the Splits Tree 4 program (Huson and Bryant, 2006). 
Abbreviations used are Ag, Anopheles gambiae ; Ap, Acyrthosiphon pisum ; Bf, 
Branchiostoma Jloridae (Amphioxus); Cbr, Caenorhabditis brenneri ; Cbri, Caenorhabditis 
briggsae ; Ce, Caenorhabditis elegans ; Ci, Ciona intestinalis; Co, Capsaspora owczarzaki ; 
Cp, Capitella sp.; Cq, Culex quinquefasciatus ; Cre, Caenorahbditis remanei ; Cs, Ciona 
savignyi ; Dp, Daphniapulex; Dm, Drosophila melanogaster ; Dr, Danio rerio (Zebra fish); 
Hr, Helobdella robusta; Hs, Homo sapiens; Lf, Loligo forbesi; Lg, Lottiagigantea; Ma, My# 
arenaria; Mb, Monosiga brevicollis; Me, Mytilus edulis; Nav, Nasonia vitripenes; Nv, 
Nematostella vectensis; Ph, Pediculus humanus; Pp, Pristionchus pacijicus (Elephant 
shark); Le, Leuconraja erinacea (Little skate); Sp, Strongylocentrotus purpuratus; Ta, 
Trichoplax sp.; Tc, Tribolium castaneum ; Tp, Trichonella sp.; TET, tetramerization 
domain; SAM, sterile alpha motif; TS, tetramerization domain and SAM domains. 
Reprinted from Rutkowski et al. (2010) with permission from CSH Press. (See Plate 
no.l in the Color Plate Section). 

vertebrate cancer and stress biology, but the roles of homologous p53 
pathways in invertebrates have scarcely been investigated. Data mining 
molluscan genomes indicates that (1) conserved genes corresponding to 
human p53 upstream regulatory and downstream transcriptionally regu¬ 
lated cancer and stress pathways can be identified and (2) these genes are 
homologues for functional components in all four of the p53-related 
vertebrate transcriptional outcomes. 
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Despite the wide taxonomic distribution of p53-like family member 
genes in both vertebrate and invertebrate animals, virtually nothing is 
known about the activity of proteins resulting from the transcriptional 
activity of these genes in most invertebrates, especially those marine 
species for which sequences have recently become available. Expression 
data are available for several marine bivalves demonstrating correlations 
between mRNA levels and fluctuations in environmental parameters or 
experimental treatments (Jessen-Eller et al., 2002; Butler et al., 2004; 
Muttray etal., 2008; Holbrook etal., 2009). These studies do not provide 
an understanding of what p53-related proteins are actually doing as a 
result of these changes. Functions of p53 proteins have been more thor¬ 
oughly investigated in the soft shell clam Mya arenaria in cancerous clam 
haemocytes (CCH) (Kelley etal., 2001; Walker etal., 2006; Bottger etal., 
2008; Walker and Bottger, 2008) and in germinomas (Van Beneden etal., 
1997; Olberding etal., 2004). In earlier literature, CCH have also been 
described as clam leukaemia or systemic disseminated neoplasia (Barber, 
2004); germinomas have been called gonadal neoplasias or gonadal can¬ 
cers. The former disease is characterized by cytoplasmic sequestration and 
inactivation of clam p53 protein mediated by over-expression of the clam 
mitochondrial hsp70 protein, mortalin (Kelley etal., 2001; Walker etal., 
2006; Bottger etal., 2008; Walker and Bottger, 2008). The germinoma is 
characterized by lower levels of p53 protein in gonadal tissues relative to 
non-tumour bearing gonads (Olberding etal., 2004). With the exception 
of CCH and clam germinoma, no other invertebrate models have been 
shown to involve p53 family members in normal or abnormal adult 
somatic cell functions related to cancer. 

Little molecular evidence regarding stress biology and relationships to 
human biology is available in relation to the function of p53 superfamily 
members in bivalves. Experimental results for M. arenaria reflect those in 
the human cancers undifferentiated neuroblastoma and colorectal adeno¬ 
carcinoma and suggest that apoptosis can be induced both transcriptionally 
(by nuclear-directed p53) and non-transcriptionally (by mitochondria- 
directed p53) following treatment with agents that generate genotoxic 
stress (DNA-damage inducing) and non-genotoxically by competing with 
p53 for mortalin binding (Kelley etal., 2001; Walker etal., 2006, 2009; 
Bottger etal., 2008; Walker and Bottger, 2008). 

While we are at an early stage in our understanding of the function of 
the p53 gene superfamily in bivalve cancer and stress biology in particular, 
even less information exists for invertebrates in general. In this review, we 
combine the results of gene-specific cDNA and genomic cloning ofp53 
family genes from the bivalves, M. arenaria (the soft shell or longneck 
clam) and Mytilusedulis (the blue mussel) with information obtained from 
data mining to provide a framework for understanding p53 pathways in 
marine bivalves. The major focus of this review is the p53 gene family and 
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its structure, function and involvement in the cancer and stress biology of 
marine bivalves and the importance of these organisms as models for 
human cancer and environmental stress biology. In particular, we empha¬ 
size human health and disease-related issues and the highly neglected 
disciplines of marine bivalve cancer and environmental stressors. 



p53 in Vertebrate Cancer and Stress Biology 


Function of the p53 protein was originally thought to rely exclu¬ 
sively on its action as a transcription factor. Many transcriptional targets 
for p53 have been identified that yield proteins involved in down- 
regulation ofp53, cellular senescence, DNA editing/repair and apoptosis 
(Horvath et al., 2007; Jegga et al., 2008; Riley et al., 2008). In approxi¬ 
mately 60% of human cancers, transcription of p53 target sequences is 
perturbed by mutations in the DNA binding domain and can result in 
disregulation of all four of these transcriptional outcomes (Yee and 
Vousden, 2005). Importantly, the apoptotic function is retained when 
mutant p53 proteins lacking the transcriptional activation domain are 
employed (Vaseva and Moll, 2008). Efforts to account for this result have 
revealed an additional non-transcriptional function involving a cyto¬ 
plasmic pool of wild-type p53 that can also participate in a mitochondrial 
pathway driving apoptosis (Talos etal., 2005). This pathway is normally 
induced by severe stress (Moll and Zaika, 2001; Vaseva and Moll, 2008) 
and initially results in translocation of wild-type p53 to the mitochondria 
where it inactivates multi-domain antiapoptotic Bcl-2 proteins resulting 
in apoptosis (Tomita etal., 2006). 

Variations in mutational levels and differential expression, phosphor¬ 
ylation and acetylation of p53 are among the factors responsible for p53 
malfunction (Greenblatt etal., 1994; Shieh etal., 1997). Although cyto¬ 
plasmic sequestration has not yet been recorded for human leukaemias or 
lymphomas, a subset of naturally occurring human cancers do retain p53 
in the cytoplasm (e.g. undifferentiated neuroblastoma, breast, retinoblas¬ 
toma (Rb), colorectal adenocarcinoma and glioblastoma cancers) and this 
phenotype can be induced in mouse NIH 3T3 and human HeLa cells 
(Wadhwa etal., 2002). In all of these cases, wild-type p53 is inactivated 
because it is retained in the cytoplasm. A variety of molecular mechanisms 
have been linked to cytoplasmic sequestration in these human cancers, 
including truncation of the nuclear localization motif receptor protein 
importin Ot (Kim etal., 2000), overactive nuclear export by an Mdm2- 
dependent pathway (Rodriguez-Lopez et al., 2001) and cytoplasmic teth¬ 
ering by foreign viral (Zhao and Liao, 2003) or local cytoplasmic proteins 
(Wadhwa etal., 1998; Nikolaev et al., 2003; Kaul etal., 2007). In human 
undifferentiated neuroblastoma and colorectal adenocarcinoma cells, 
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mortalin protein is responsible for cytoplasmic tethering when the latter 
protein is over-expressed (Dundas et ah, 2005; Walker et ah, 2006). 
Dundas et ah (2005) also indicate that high levels of mortalin expression 
are correlated with poor clinical outcome for colorectal cancer. A recent 
study of primary and secondary glioblastomas suggests that mortalin and 
possibly other tethering molecules (e.g. cullin 7 or PARC) may also be 
responsible for cytoplasmic sequestration in these cancers (Nikolaev etai, 
2003; Nagpal et ah, 2006). 

The effects of environmental stressors on the biology of human cancer 
are addressed in a number of studies (Meisner etah, 1993; Nawrot et ah, 
2006; Soliman etah, 2007; Kampa and Castas, 2008). We will emphasize 
those that have parallels in marine bivalve haemocyte and germinoma 
cancers. 



Structure of p53 Family Member Proteins in Bivalves 


Approximately 17 kb of M. arenaria genomic DNA indicate that 
both soft shell clam p53 and p63/73 proteins are encoded by a single gene 
(GenBank Accession number FJ041332). Clam p53 (earlier called Map53 
by Kelley et ah, 2001) consists of 10 exons and clam p63/73 (earlier called 
Map73 by Kelley etah, 2001) is comprised of 12 exons. Of these, the first 
nine exons encode the core sequence for both proteins (Fig. 1.2a and b). 
The molecular organization of the clam gene is similar to its vertebrate 
counterparts: (a) exon 1 is non-coding with the translation initiation site 
contained within exon 2; (b) the presence of a relatively large intron at the 
5 ' end of the gene, located between exons 3 and 4 and (c) intron/exon 
boundaries are similar between clam p63/73 and other p53 gene family 
members. 

The distinction between p53- and p63/73-like sequences is based on 
the presence of a SAM (sterile alpha motif) domain in the carboxy 
terminus of the latter (Thanos and Bowie, 1999; Nedelcu and Tan, 
2007). In at least three molluscs, cDNAs for both proteins have been 
cloned and core sequence for the predicted p53- and p63/73-like proteins 
is virtually identical (Fig. 1.2b). In all three of these cases, the core 
sequence includes the following highly conserved regions: the transcrip¬ 
tional activation domain (the clam TAD is 67% conserved with human 
p53; it is only 42% conserved with human p73 and only 33% conserved 
with human p63), Mdm2 binding site (including all essential Mdm2 
binding residues—F, L, W, L), ATM phosphorylation site (S15), pro¬ 
line-rich domain, DNA binding domains (DBDs II—V), nuclear import 
(3 NLSs) and export (1 NES) signals and the tetramerization domain. In 
the soft shell clam, the core sequence is a structural mosaic of the corre¬ 
sponding human proteins, with the TAD closely resembling human p53 
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LEIAALLPQNQLQSLKQKQVEVQRQ 
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SSSIQAATSSARIPAIAATPATTYVQQGVTTTSSDGKQLTMPFNTQELVQVTSSDVSHDGAVPQPIKEETIQNDMQDNSVSTWLNALGLGAYIDGFH 

EQNLYSLLQLDDFSLDDLAKMKIGNSHRNKIWKSLLELRNQGFTTAESQDSLAKQASSTSTISLASQGSISQNSTYNPGFYEVTRHTFKHTISLTKE 
ERHVQTRS SASKAD 

Map53 Extension: 

WLTNVLAKEGKSRLIKKKHRPGKIIRHPLK 

c 

Map53/63/73 intron 3 ..ataaataaaaatgattaagtttgagagaactggcttcacaacttatagggtgagtacgc... 

, Mill II llllllllllllll II Mill 111111 
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MIKFERTGFTTYR 

Figure 1.2 (a) Genomic structure of p53/63/73 in M. arenaria. Boxes represent 

exons and V indicates intron regions (not drawn to scale). The S^UTR. is grey, 
the core sequence for both clam p53 and p63/73 proteins is green, the p53 extension 
is blue and the p63/73 extension is red. The yellow arrow denotes the location of a 
predicted DN isoform; (b) amino acid sequences for clam p53 and p63/73 with 
intron/exon boundaries marked (▼, clam p53; V, human p53). Functional 
domains are indicated by horizontal lines above the sequence. Conserved amino 
acids implicated in Mdm2 binding are underlined. Sequence highlighted in yellow 
denotes the predicted DN isoform. Genomic DNA was purified from flash-frozen 
gill tissue of M. arenaria as previously described (Barker et al., 1997). DNA template 
was prepared using the Universal Genome Walker Kit according to the 
manufacturer’s instructions (Clontech). Clam sequence-specific PCR. primers 
were designed based on the clam p53 and p63/73 coding sequences (GenBank 
Accession numbers AF253323 and AF253324, respectively) and used in 
conjunction with adapter primers (Clontech). Amplified products were cloned 
into the pCRII vector (Invitrogen) and sequenced in both directions at the 
University of Maine DNA Sequencing Facility. Contiguous sequences were 
assembled using Sequencher version 4.1.2 software (Gene Codes Corp., Ann 
Arbor, MI) and the open reading frames identified; sequence for genomic 
DNA was deposited under GenBank Accession number FJ041332 and 
(c) sequence for M. edulis DNp63/73 (DQ865151.1) was aligned to FJ041332 
using Spidey (http://www.ncbi.nlm.nih.gov/spidey/) to elucidate a sequence and 
intron/exon boundaries for a predicted DN isoform in M. arenaria , and was found to 
be 87% identical to M. edulis at the nucleic acid level and 100% at the protein level and 
located in intron 3. (See Plate no.2 in the Color Plate Section). 
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(>67%) and the DBDs resembling human p73 (75—88%). Clam proteins 
have either short (p53) or long (p63/73) C-terminal extensions. Clam 
genomic DNA sequence demonstrates that unlike humans, where p53, 
p63 and p73 proteins are generated from distinct genes on different 
chromosomes, clam p53 and p63/73 mRNAs result from splice variants 
of a single gene. M. edulis and Loligoforbesi p53- and p63/73-like proteins 
also have identical core sequences and may be splice variants of one gene 
as in M. arenaria (Fig. 1.2). However, Mytilus trossulus has distinct genes 
that code for unique p53 and p63/73-like sequences (Muttray etal., 2008; 
Vassilenko etal., 2010). 

Invertebrates the DeltaNotchp63/73-like (DN) splice variant ofp73 
is commonly expressed. The distinction between p73 and DeltaNotch 
p73-like (DN) sequences is based on the replacement of the long TA 
domain with a truncated ‘DN’ domain (Ishimoto etal., 2002). Among 
invertebrates, only M. trossulus and M. edulis have been found to express 
DNp63/73-like mRNAs (Muttray etal., 2007). Unfortunately, the lack 
of genomic p53 sequence data from Mytilus sp. limits our understanding 
of how this DNp63/73 isoform may be spliced from the gene. However, 
given the overall sequence homology between M. arenaria and Mytilus 
p63/73 isoforms, some inferences can be made by aligning the M. are¬ 
naria genomic p53 sequence with the Mytilus DNp63/73 mRNA 
sequence. It appears that the nucleotide sequence encoding the truncated 
DN domain from Mytilus is present in the long intron located between 
exons 3 and 4 of the M. arenaria genomic sequence and is 87% identical 
between Mya and Mytilus (Fig. 1.2c). It is therefore likely that an addi¬ 
tional exon encoding a DN isoform may also exist in Mya and that the 
DN isoform could be expressed under conditions that have not been 
explored by isoform cloning. Expression levels for the DN isoform are 
generally low in normal mussel haemocytes and variable in cancerous 
mussel haemocytes (Muttray etal., 2010). The function oftheDN isoform 
in mussels is unknown. 

Four residues (F19, L22, W23 and L26) that are critical for hydro- 
phobic binding interactions with Mdm2 are located in the TAD of human 
p53 (Kussie etal., 1996; Dobbelstein and Roth, 1998; Holbrook etal., 
2009). Protein modeling predicted homologous Mdm2 binding sites in 
Mya/Mytilus p53 (Kelley etal., 2001; Muttray etal., 2005; Holbrook etal., 
2009) (Figs. 1.3 and 1.4). These residues are also present in p53 sequences 
of the Pacific oyster Crassostreagigas (Farcy etal., 2008) and the surf clam 
Spisida solidissima (Cox et al., 2003). In molluscs, the putative Mdm2 
binding site is identical for both p53- and p63/73-like proteins. An 
Mdm2 homologue has also been cloned from the mussel, M. trossulus 
(GenBank Accession number ADF80285) and shown to bind p53 protein 
although no functional data are available (Muttray etal., 2010). While 
Mdm2 homologues have also been identified from the genomes of the 
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Figure 1.3 Structural analysis of M. arenaria p53 shows conservation of a putative 
Mdm2 binding site. Representation of a portion of the human Mdm2 protein (ribbon 
structure) bound to the minimal binding site of the human p53 protein (stick 
structure). Modified from Kussie et al. (1996) using RasMol software. Four 
residues critical for hydrophobic binding interactions (F19, L22, W^23 and L26, 
shown in green) are conserved between human p53 and clam p53. Two 
potentially important changes are the substitutions in M. arenaria p53 of glutamic 
acid and tyrosine, for S20 (red) and D21 (blue) in human p53, respectively 
(Holbrook et al ., 2009) (See Plate no.3 in the Color Plate Section). 


placozoan, Trichoplaxadhaerens and the northern deer tick, Ixodes scapularis, 
p53 in these species does not have a conserved Mdm2 binding site (Lane 
et al., 2010a,b) (Figs. 1.3 and 1.4). 

Given the conservation of key domains in human and clam p53 
proteins (Kelley eta!., 2001), the hypothesis that the two proteins func¬ 
tion in a similar manner was addressed experimentally (Holbrook etal., 
2009). Plasmids expressing either human p53 or M. arenaria p53 were 
introduced by transient transfection into the p53-null HI299 cell line 
derived from human liver. Functionality was assessed by monitoring the 
p53/mdm2 feedback loop and expression of p53-mediated downstream 
markers of growth arrest (e.g. CIP1) under non-stressed conditions. 
Expression of human p53 spontaneously induced expression of the 
cyclin/CDK binding protein CIP1 while clam p53 expression did not 
(Fig. 1.5). 

Expression of M. arenaria p53 transfected into the human p53-null 
HI299 cell line induced expression of human Mdm2 to approximately 
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Figure 1.4 Species tree comparing p53 family members in sequenced genomes for 
invertebrates and humans ( Homo sapiens). Comparisons include from left to right 
(a) percent similarity of amino acids in the TADs of invertebrate p53 genes compared 
to human p53 TAD, (b) presence or absence of Mdm2 binding site homologous with 
that of human p53, (c) number and nature of p53 family member genes present; 
(d) existence or absence of core sequence and (e) taxonomy. Human and M. arenaria 
p 53-like sequences were used to retrieve homologous gene models from the following 
genome projects located at the JGI Eukaryotic Genomics Web site (http://genome.jgi- 
psf.org/euk_curl.html for Capitella sp.; D. pulex; C. intestinalis; B.floridae; M. brevicollis; 
Trichoplax adhaerans and N. vectensis) and at WormBase (http://www.wormbase.org/ 
db/searches/blast_blat) for C. elegans; at Flybase (http://flybase.org/blast/) for D. 
melanogaster and at (http://www.spbase.org/SpBase/wwwblast/blast.php) for 5. 
purpuratus. Protein alignment was determined with ClustalX software. The species 
tree was constructed using PhyloWidget (Jordan and Piel, 2008) based on topology of 
Dunn et al. (2008). 


the same level as did human p53, suggesting that Mya p53 is also able to 
participate in an Mdm2 feedback loop. Conversely, the addition of human 
Mdm2 RNAi to H1299 cells induced the expression of human p53 and 
Mya p53 (Holbrook et al., 2009). The M. trossulus Mdm2 protein also 
binds Mytilus p53 as shown in two independent assays, the yeast two- 
hybrid assay and immunoprecipitation (Muttray etal., 2010). An obvious 
limitation of this in vitro assay is that it was done at 37 °C, much higher 
than ambient temperatures typically experienced by the clam, which may 
have a significant effect on protein folding and protein/protein interac¬ 
tions. Taken together, these data suggest that the M. arenaria and Mytilns 
p53 proteins share some functional similarity with human p53 as well as 
with other invertebrates, positioning the molluscs at a critical juncture in 
evolution of this gene family. 
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pcDNA-Human p53 pcDNA-M. arenaria p53 Vector alone 




Figure 1.5 Human p53 elicited expression of CIPlin p53-null H1299 cells, but M. 
arenaria p53 did not. Representative Western blot and summary densitometry data 
show that human p53 expressed in HI299 cells (black bars) elicited a robust CIP1 
response that remains elevated over the 48-h observation period. Neither M. arenaria 
p53 (grey bars) nor vector alone (pcDNA3.1/His, Invitrogen; white bars) altered 
CIP1 expression. There are no significant differences in CIP1 expression over time 
within any treatment group (P > 0.05; N = 3). Values were pooled for between- 
treatment comparisons; different letters denote significant differences; M, marker 
lane (after Holbrook et al ., 2009). 



Data Mining the Molluscan Genome for p53 Upstream 
Regulatory and Downstream Transcriptional Targets 


The study by Nedelcu and Tan (2007) compares the DNA binding 
domains of invertebrate p53 family members, but does not consider the 
TADs or p53-related regulatory or transcriptional targets for these species. 
Since the TAD of soft shell clams and mussels is 67—69% conserved with 
that of humans it is possible that clam p53- and p63/73-like proteins may 
share transcriptional activating functions with corresponding human pro¬ 
teins (Fig. 1.4). Considerable experimental evidence supports this 
hypothesis for the soft shell clam and will be discussed below (Kelley 
et al., 2001; Walker et al., 2006; Bottger ef al., 2008; Walker and Bottger, 
2008; Holbrook etal., 2009). 

Three genome sequencing projects have been completed for molluscs 
(for the archeogastropod Lottia gigantea; http://genome.jgi-psf.org/ 
Lotgil/Lotgil.home.html; for the sea hare, Aplysia californica: http:// 
genome.ucsc.edu/cgi-bin/hgGateway?hgsid 

+ 155184841&clade=other&org=Sea+hare&db=0 and very recently for 
the oyster, C.gigas, by the Oyster Genome Sequence Map Project, China, 
unpublished). Undoubtedly, completion of bivalve genomes and partic¬ 
ularly functional studies of p53-related proteins will be of great advantage 
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in developing a more thorough understanding of the biology the p53 
family of genes in cancer and stress biology. In this chapter, our analysis of 
genes in the limpet (Lottia gigantea ) genome that are either upstream 
regulators or are transcriptional target genes for molluscan p53 provides 
important information critical for designing experiments to determine 
the extent of functional similarity between bivalve and human p53- 
related proteins in their respective normal or cancerous cytoplasmic and 
nuclear environments (Fig. 1.6). 

Upstream regulatory oncogene and stress-induced genes important in 
human biology that are also present in the molluscan genome (Fig. 1.6) 
include Rb and E2F1 and ATM, ATR, ms and Chk2. Certain genes are 
missing (e.g. ARF, alternate reading frame), although HES1 and HEY1 


Human p53 pathways 


Oncogene—driven 
Rb Eg FI 

ARF & Hesl and Hey 1 

^^>|Mdm2 


Cellular scenescence 

CDKN1A SFN 
TP313 CDC2C 


DNA damage-driven 

ATM/ATR 

Mortalin 
(over-expressed 
in neuroblastoma and 
Puma/Noxa colon carconoma) 



DDIT4, TRIM22. GADD45, DDB2, BREC1. PCNA, TOPOII 

DNA editing and repair 


Proposed molluscan p53 pathways 


Oncogene-driven 


DNA damage-driven 

ATM/ATR 


No ARF , Hesl and Hey 1 

**'" ^|Mdm2 


Cellular scenescence 

CDKN1A SFN 
TP313 CDC2C 



Mortalin 
(Over-expressed 
in clam hemocyte 
cancer) 


DDIT4. TRIM22, GADD45. DDB2, BREC1, PCNA, TOPOII 

DNA editing and repair 



Figure 1.6 Comparison of genes present in H. sapiens and genomic DNA for L. 
gigantea p53 upstream regulatory and downstream transcriptionally regulated cellular 
senescence, DNA editing and repair and apoptosis pathways. Lottia genome searches 
were performed using BLAST for human sequences for p53 regulatory and 
transcriptional target genes from the following genome projects located at the JGI 
Eukaryotic Genomics Web site (http://genome.jgi-psf.org/euk_curl.html for 
Capitella sp.; D. pulex ; C. intestinalis; B. jloridae; M. hrevicollis; T. adhaerans and N. 
vectensis) and at WormBase (http://www.wormbase.org/db/searches/blast_blat) 
for C. elegans ; at Flybase (http://flybase.org/blast/) for D. melanogaster and at 
(http://www.spbase.org/SpBase/wwwblast/blast.php) for S. purpuratus. Protein 
alignment was determined with ClustalX software. 
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that behave in an analogous manner by blocking Mdm2 are highly con¬ 
served (Huang et al., 2004). The ubiquitin ligases, Mdm2, COP1 and 
Pirh2 are all present in the molluscan genome as is the ubiquitin-specific 
peptidase HAUSP. Unrelated p53 binding proteins such as mortalin, 
PARC, cullin 7 and Bcl-xl also occur. 

Genes are also found in the molluscan genome that are transcriptional 
targets for p53 in humans and are associated with outcomes in cellular 
senescence (e.g. CDKN1A, SFN, TP53I3, CDC25C), DNA editing and 
repair (e.g. DDB2, DDIT4, GADD45(Y, TRIM22, ribonucleotide 
reductase) and apoptosis (e.g. FAS, Bcl-2, Bcl-xl, Mcl-1, BAX, BAK, 
procaspases, 3, 7 and 9, APAF1, ICAD and CAD). Some genes important 
for mammals in regulation (e.g. CIP1) or apoptosis (e.g. all BH3-only 
proteins) are absent from the molluscan genome (Fig. 1.6). The inability 
of clam p53 to up-regulate CIP1 in the H1299 cell line described above 
(Holbrook etal., 2009) might be related to the absence of CIP1 from the 
molluscan genome (as is also the case for all non-chordate genomes). 
Clam p53 protein expressed in human cells might not recognize the 
promoter for the human CIP1 gene and therefore cannot initiate its 
expression. 



Bivalve Molluscs as Models for Animal Cancer 
and Stress Biology 


5.1. Haemocyte cancer in the soft shell clam and the 
involvement of the mitochondrial hspzo protein 
(mortalin) in cytoplasmic sequestration of P53 

The soft shell clam, M. arenaria, naturally develops a diffuse cancer of its 
haemocytes at high incidence (1—20% of populations) (Barber, 1996; 
Walker et al., 2009). This fatal disease is found in clams at multiple 
sites along the coasts of all states between Maryland and Maine as 
well as in eastern Canada (Bottger et al., in preparation a). 
Cancerous clam haemolymph contains 100% cancerous haemocytes 
(~5 X 10 8 cells ml ) with an average diameter of 7—10 |xm and a mean 
nuclear to cytoplasmic volume ratio of 1:1. These nearly round, mitotic 
haemocytes have a monotonous appearance, lack large pseudopodia, 
attach only loosely with thin cytoplasmic projections to plastic or glass 
and are not motile nor phagocytic (Fig. 1.7a). They contain membrane 
bound vacuoles (5—15/cell) that stain intensely for lipid with oil red O 
(Fig- 1.7b). 

Protocols for maintaining cancerous haemocytes from the soft shell 
clam in vitro and for amplifying them in vivo have been developed and are 
described in detail by Walker etal. (2009). These were the first successful 
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Figure 1.7 Clam haemocyte cancer: (a) interference image of freshly collected 
cancerous clam haemocytes; notice dumb-bell shaped haemocyte involved in 
cytokinesis near the centre; (b) Cytospin of CCH stained with oil red O; notice 
mitotic division figure on the left of the image; (c) Quantum-dot immunocytochemical 
localization in cytoplasmically sequestered p53 in an individual CCH; (d) quantum-dot 
immunocytochemical localization in cytoplasmic mortalin in the same CCH; (e) dual 
quantum-dot immunocytochemical localization of cytoplasmically sequestered p53 
and mortalin in the same CCH—yellow indicates co-localization of the two proteins 
(scale bar in (a) and (b) = 20 pirn and in (c—e) = 10 pim) and (f) co- 
immunoprecipitation of clam p53 and mortalin in normal (NCH) and CCH. Std, 
protein standard; Cl and C2, first and second elute from CCH lysates; N1 and N2, 
first and second elude from NCH lysates; NA, first elude from gel loaded with NCH 
and CCH lysate mix (negative control); Q, first elude from quenched gel loaded with 
NCH and CCH lysate mix (negative control); Ccyt, cytoplasmic protein extract from 
CCH and Cnuc, nuclear protein extract from CCH. (See Plate no.4 in the Color Plate 
Section) 


methods at maintaining malignant cells from any invertebrate in vitro. 
They provide ready experimental access to simultaneous in vitro and in 
vivo versions of this naturally occurring cancer model from the same 
organism for use by both biomedical and environmental researchers 
(see Appendix A of this review for a significantly updated version of the 
clam culture medium originally presented in Walker etal., 2009). 

Cancerous clam haemocytes express a highly conserved homologue 
for human wild-type p53 protein (GenBank Accession number 
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AF253323) that is rendered non-functional through sequestration in 
the cytoplasm by mortalin protein (the mitochondrial Hsp70), when 
the latter is over-expressed (Fig. 1.7c—e). Expression of clam mortalin is 
1634 times higher in CCFI than in normal clam haemocytes (NCH) as 
determined by quantitative PCR (Bottger et al., 2008). Cytoplasmic 
sequestration does not occur in NCH (Kelley et al., 2001; Walker et al., 
2006). A similar phenotype involving cytoplasmic sequestration of wild- 
type p53 has been observed in an unrelated group of human cancers, 
including undifferentiated neuroblastoma and colorectal adenocarcinoma 
(Moll et al., 1992, 1995, 1996; Kelley etal., 2001; Dundas et al., 2005; 
Walker et al, 2006, 2009; Walker and Bottger, 2008). 

These were the first structural and functional data for p53 and p63/73 
mRNAs and gene products in a naturally occurring, non-mammalian 
disease model for human cancers and also for any bivalve mollusc 
(Kelley etal, 2001). Based on the substantial homology of the TAD in 
these proteins (65—69% identity), it was suggested that clam p53- and 
p63/73-like proteins might both share transcriptional functions with 
human proteins. Subsequent studies have indicated that as in humans, 
both transcriptional and non-transcriptional p53-dependent mechanisms 
are involved in soft shell clam haemocyte cancer (Kelley etal., 2001; 
Walker et al., 2006; Bottger et al., 2008; Walker and Bottger, 2008; 
Holbrook et al., 2009). These will be discussed in detail in the following 
sections. 

5.2. Transcriptional induction of apoptosis by soft 

shell clam P53 

Genotoxic stress induced by treatment of cancerous clam haemocytes 
with etoposide promotes denovo transcription of wild-type p53 (602-fold 
increase over untreated CCH), translocation ofp53 protein to the nucleus 
(Fig. 1.8a and b), followed by DNA damage (Fig. 1.8b and c) and apo¬ 
ptosis (Bottger etal., 2008). The most attractive interpretation of these 
data is that treatment with etoposide elevates p53 levels in the cytoplasm 
of CCH and that de novo p53 protein overwhelms mortalin tethering, 
resulting in nuclear translocation of some p53 and apoptosis. After treat¬ 
ment with etoposide, transcriptional targets of p53 would be active in 
evaluating and repairing DNA, or failing that induction of the apoptosis 
of CCH that we have observed. 

5.3. Non-transcriptional induction of apoptosis 

by soft shell clam P53 

Treatment of CCH with the nuclear pore blocker wheat germ agglutinin 
(WGA) demonstrated that clam p53 is found only in the cytoplasm 
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Figure 1.8 Transcriptional initiation of apoptosis by soft shell clam p53. 
Distribution of p53 before and after treatment of CCH with etoposide; (a) and 
(b) show the distribution of p53 using an M. arenaria polyclonal antibody 
(Kelley et al 2001). (a) Zero time control for CCH in which clam p53 protein is 
sequestered in the cytoplasm with a perinuclear distribution (notice the clear nuclei, 
n), scale bar =10 |xm; (b) cancerous clam haemocytes from the same individual as 
Fig. 1.7a in which clam p53 protein has re-localized to the nucleus following 
treatment with 0.075 mM etoposide for 18 h, scale bar =10 |xm; (c) comet assay 
showing no evidence of DNA damage in control normal clam haemocytes after 18 h 
following sham injection with carrier minus etoposide; (d) comet assay showing 
extensive DNA damage in leukaemic clam haemocytes from the same individual 
as in part (a) following treatment with 0.5 mm etoposide for 18 h. 


(Fig. 1.9a), with some mitochondrial localization (Fig. 1.9d). When 
treated with MKT-077 (a chemotherapeutic rhodocyanin dye) only (no 
pre-treatment with WGA), clam p53 was initially found in the cytoplasm, 
but was also translocated into the nucleus at 6 h (Fig. 1.9b) with very little 
accumulation at the mitochondria (Fig. 1.9e). Pre-treatment of CCFd 
with WGA followed by post-treatment with MKT-077 demonstrated 
that clam p53 is found only in the cytoplasm (Fig. 1.9c), with increased 
mitochondrial localization of p53 over the results from treatment with 
MKT-077 alone (Fig. 1.9f). Treatment with MKT-077 induced apopto¬ 
sis, as detected by the TUNEL assay (results not shown). The most 
attractive interpretation of our results is that MKT-077 competes with 
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Non-transcriptional initiation 

of apoptosis 

by soft 

shell 

clam p53. 


Distribution of p53 in cytoplasmic, nuclear and mitochondrial protein extracts 
from CCH using an M. arenaria polyclonal antibody (Kelley et al., 2001) following 
treatment with MKT-077 for 6 h. Treatment with MKT-077 occurred both with and 
without WGA pre-treatment. Pre-treatment of CCH with WGA to block nuclear 
import employed FITC-labelled WGA (Sigma) transfected into CCH using the 
Chariot protein delivery system (Active Motif). WGA was allowed to couple with 
the Chariot compound for 30 min at room temperature to form complexes that were 
incubated for 1 h with CCH suspended in culture medium at concentrations of 
6 X 10 5 CCH/100 |Jil medium. Cytoplasmic, nuclear and mitochondrial protein 
extracts from CCH were treated as follows: for nuclear/cytoplasmic extracts— 
(a) WGA only, (b) 3.5 |xM MKT-077 only, (c) pre-treatment with WGA 
followed post-treatment with 3.5 |xM MKT-077; and for mitochondrial extracts— 
(d) WGA only, (e) with 3.5 p,M MKT-077 only and (f) with 3.5 |xM MKT-077 
following pre-treatment with WGA. Treatment with WGA blocks nuclear access for 
p53 and results in its translocation to the mitochondrion where apoptosis of CCH 
follows. Wfren normal clams are treated identically, p53 is transported to the 
mitochondria, but apoptosis does not result. Std, protein standard; nuc, nuclear 
protein; cyt, cytoplasmic protein. 


p53 for the mortalin p53 binding site. As a result unbound p53 is 
available to be transported to the nucleus or mitochondria. Since 
nuclear access is blocked, p53 is only transported to the mitochondria 
where it binds to and inactivates antiapoptotic Bcl-2 proteins in the 
mitochondrial membrane. Pro-apoptotic Bax proteins from the cyto¬ 
plasm are then recruited to the mitochondria and participate in mito¬ 
chondrial catastrophy. 

5.4. Haemocyte cancer in other bivalves 

In some locations, the blue and bay mussels, M. edulis and M. trossulus, also 
develop a cancer-like, proliferative disease of their haemocytes at high 
prevalence that is morphologically similar to that described above for M. 
arenaria (Bower, 1989; Barber, 2004). While mussel haemocyte cancer 
has received attention as a tool for environmental monitoring (Muttray 
et al., 2008; Pantzartzi etal., 2010), studies on the molecular basis for the 
disease have only recently begun (Muttray etal., 2005). The elucidation of 
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the structure of mussel p53 and p63/73 genes (Muttray etai., 2007) has 
provided the basis for future research into the function of their gene 
products in mussel haemocyte cancer. M. arenaria and Mytilus spp. p53 
and p63/73 share 69 and 57% sequence similarity at the protein level, 
respectively, and it is possible that the genes and gene products of the two 
species also share some of the functional properties that contribute to 
haemocyte cancer. Thus far, expression analysis has shown that p53 as 
well as an Mdm2-like gene product are expressed at higher levels in 
cancerous mussel haemocytes when compared to normal haemocytes 
(Muttray et al., 2010). The increase in p53 also results in an increase in 
transcription ofMdm2, similar to what has been reported for mammalian 
p53—Mdm2 feedback interactions. Indeed, Mytilus p53 and Mdm2 pro¬ 
teins interact, as has been described above. In addition, the DNp63/73 
isoform was expressed at higher levels in mussel haemocyte cancer. This is 
similar to some human cancers where it has been suggested as a biomarker 
of poor prognosis (Muttray etai, 2008). Because p53 was expressed at 
higher levels in cancerous haemocytes that escape apoptosis, the hypoth¬ 
esis that p53 may have acquired inactivating mutations was tested. 
Although specific single nucleotide polymorphisms (SNPs) were detected 
in cancerous mussel haemocytes, they were silent and thought unlikely 
to affect protein function (Vassilenko et al., 2010). However, variations 
in mRNA secondary structure and hence rates of transcription, and a 
potential change in splicing patterns may affect protein levels. 


5.5. Germinomas in bivalves 

Bivalve germinomas (gonadal cancer; germ cell cancers) were first 
described over four decades ago in hard shell clams (Mercenaria spp.) 
(Yevich and Barry, 1969). Early stages of the disease exhibit clusters of 
neoplastic cells that appear to originate basally in gonadal follicles 
(Rhodes and Van Beneden, 1996). As the disease progresses, density of 
the stromal extracellular matrix increases, accompanied by infiltrating 
phagocytic cells. In late stages, neoplastic cells fill the gonadal lumen, 
eventually obliterating the follicular structure (Fig. 1.10a—d). Unlike hae¬ 
mocyte cancer, the germinoma does not appear to be fatal. The drastic 
morphological changes and the apparent loss of germ cells suggest that 
there could be a significant impact on reproductive output at both the 
individual and population levels (Barber, 1996). 

Investigation of the molecular basis of the gonadal cancer in M. are¬ 
naria has focused on the effects of dioxin, aryl hydrocarbon (AhR) 
and the herbicide 2,4-D and molecular mechanisms that might mediate 
their toxicity (Butler etai., 2001, 2004; Van Beneden, 2005). Aberrant 
regulation of p53 gene has also been investigated. The level of p53 
protein under feedback control by Mdm2, is normally elevated in 
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Figure 1.10 Comparison of non-tumour and tumour-bearing gonadal tissue in 
female M. arenaria. (a) Typical gross appearance of mature gonad; note large, 
creamy appearance (arrow); (b) clam with advanced germinoma; note brown 
mottling (arrow); (c) histological section showing normal ovarian follicle 
containing mature eggs (H&E stain), scale bar =10 jxm; (d) histological section 
showing typical small, undifferentiated cells of advanced germinoma (H&E stain), 
scale bar = 10 pan. (For interpretation of the references to colour in this figure 
legend, the reader is referred to the Web version of this chapter.) 


response to cellular stress. In human cervical cancers, abnormal degra¬ 
dation of p53 is mediated by the E3 ubiquitin ligase (Huibregtse etal., 
1991). Since altered levels of the clam E3 homologue mRNA were 
detected earlier in laboratory exposure studies; the E3 cDNA sequence 
was determined (GenBank Accession number AF154109) (Kelley and 
Van Beneden, 2000). Protein expression levels of the clam homologues 
of E3 and p53 were examined in tumour- and non-tumour-bearing 
clams by Western blot analysis (Kelley and Van Beneden, 2000; 
Olberding etal., 2004). E3 protein expression was found to be signifi¬ 
cantly higher in tumour-bearing females compared to non-germinoma- 
bearing females (Fig. 1.11). Conversely, the level of p53 protein was 
found to be significantly lower in germinoma-bearing females com¬ 
pared to non-germinoma-bearing individuals (Fig. 1.11). These find¬ 
ings supported the hypothesis that unscheduled E3-mediated degrada¬ 
tion ofp53 releases normal cell cycle arrest at Gj, allowing abnormal cell 
proliferation and tumour formation. However, the direct interaction of 
p53 and E3 was not confirmed by protein/protein binding studies 
(Olberding et al., 2004). 
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Figure l.ll Clam germinoma tissue exhibited significantly higher levels of E3 
protein concomitant with significantly lower levels of p53 protein as compared to 
normal gonadal tissues. Summary graphs represent densitometry scans from the 
Western blots normalized to actin. Non-tumorigenic samples consisted of both 
male and female tissue; there were no statistically significant gender differences for 
either MaE3 or Map53 protein levels. Germinoma proteins were isolated from 
females; no male germinoma tissue was available for analysis. * significant 
difference; P < 0.05; error bars represent SD) (after Olberding et al ., 2004). 


5.6. Environmental stressors and potential links to bivalve 
haemocyte and germinoma cancers 

Linkages between environmental toxins and cancer are often difficult to 
demonstrate for lack of an appropriate model system (Van Beneden et al., 
1997; Van Beneden, 2005). Bivalve molluscs, however, are excellent 
candidates for such models. First, they are long-lived and sessile and filter 
feeding exposes them to both suspended and soluble compounds. 
Production of reactive oxygen species (ROS) is also a common response 
to environmental stressors in many marine invertebrates including M. 
arenaria (Abele et al., 2002; Lesser, 2006) and Mytilus spp. 
(Winston et al., 1996). All other parameters being equal, an increase 
in antioxidant capacity is often associated with higher exposures to 
xenobiotics (Frenzilli etal., 2004). 

5.6.1. Contaminants: field studies 

Based on limited, highly site-specific data, several environmental con¬ 
taminants have been circumstantially linked to haemocyte cancer in M. 
arenaria (Brown, 1981). Reports show significant correlation with oil 
(Yevich and Barry, 1969; Barry and Yevich, 1975; Brown etal., 1977; 
Lowe and More, 1978; Yevich and Barszcz, 1978; Walker etal., 1981), 
polyaromatic hydrocarbons (PAHs) (Lowe and More, 1978), polychlori¬ 
nated biphenyls (PCBs) (Reinisch et al., 1984; Bottger et al., in 
preparation b), industrial and municipal waste (Walker et al., 1981; 
Reinisch etal., 1984; St-Jean etal., 2005), pesticides (Kagley, 2003) and 
heavy metal pollution (Walker etal., 1981; Kagley, 2003) (Fig. 1.12). 
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Figure 1.12 Prevalence of haemocyte cancer in M.arenaria correlated with different 
contaminants. Results were compiled and averaged for oil (Yevich and Barry, 1969; 
Barry and Yevich, 1975; Brown et al., 1977; Lowe and More, 1978; Yevich and 
Barszcz, 1978; Walker et al., 1981), PAHs and PCBs combined (Bottger et al., in 
preparation b; Reinisch et al., 1984), industrial and municipal waste (Walker et al., 
1981; Reinisch et al., 1984; St-Jean et al., 2005) and heavy metal pollution (Walker 
et al., 1981; Kagley, 2003). Data reflect frequencies for haemocyte cancer and 
different contaminants averaged from literature values for a number of sites. 


Germinomas in the soft-shell clam, M. arenaria, were first observed 
following the March 1971 oil spill at Long Cove, ME (Dow and Hurst, 
1975). In the aftermath of the spill, it was speculated that they might be 
linked to exposure to petroleum hydrocarbons. Since both haemocyte 
and germinoma cancers have also been detected in unpolluted sites, it is 
not clear if any of these factors or any group of them is responsible for 
onset or progression of these diseases. 

5.6.2. Contaminants: laboratory studies 

Despite the overwhelming number of reports linking carcinogenic che¬ 
micals in the environment and clam haemocyte cancer, it has been diffi¬ 
cult to show clear cause—effect relationships in the laboratory (Harsbarger, 
1967). Haemocyte cancer did develop in Uniopictorum exposed to PAHs/ 
PCBs (Khudolei and Sirenko, 1978) but did not following the same 
treatment in M. edulis (Rasmussen, 1986), M. arenaria (MacCallum etal., 
2003) and Macoma baltica (Tay, 2003). 

Laboratory exposure studies under the direction of George Gardner at 
the EPA Laboratory (Narragansett, RI) were conducted in an attempt to 
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mimic the 1971 Maine oil spill. Adult M. arenaria were exposed to sedi¬ 
ments from Long Cove, ME, contaminated with petroleum hydrocar¬ 
bons; exposures, however, failed to induce germinomas. A companion 
study (Gardner et ah, 1991), however, reported that oysters ( Crassostrea 
virginica ) exposed to Black Rock Harbor sediments heavily contaminated 
with PAHs and PCBs developed a variety of neoplastic lesions, predom¬ 
inantly kidney lesions (9.3%) and germinomas (1.3%). The potential link 
of germinoma and environmental contaminants has been further inves¬ 
tigated in a series of controlled laboratory exposure studies (Rhodes and 
Van Beneden, 1996; Butler etal., 2004). Van Beneden and Gardner tested 
the hypothesis that either dioxin or herbicides may induce the germinoma 
in M. arenaria (Rhodes and Van Beneden, 1996; Butler etal., 2004). Clams 
were exposed to either (a) 3 H-TCDD (2,3,7,8-tetrachlorodibenzo-p- 
dioxin) in the presence or absence of the initiating carcinogen diethylni- 
trosamine (DEN) or (b) 2,4-D (Kelley and Van Beneden, 2000). 

Further investigation revealed that laboratory exposure to dioxin eli¬ 
cited early changes (<2 weeks post-exposure) in gene expression patterns 
that included induction of biosynthetic pathways, signal transduction 
cascades and ubiquitin-mediated proteolytic pathways (Rhodes and Van 
Beneden, 1996, 1997; Kelley and Van Beneden, 2000). Using differential 
display PCR (a differential mRNA expression screen, technology which 
preceded today’s microarrays), Rhodes and Van Beneden reported dif¬ 
ferentially expressed mRNA sequences in TCDD-exposed animals com¬ 
pared to unexposed controls (Rhodes and Van Beneden, 1997; Kelley and 
Van Beneden, 2000). So far, laboratory studies to directly link contami¬ 
nants with gonadal tumour induction have been largely unsuccessful. In 
summary, laboratory exposure to potential carcinogenic chemicals does 
not reproduce observations made in the field. 

5.6.3. Viral involvement 

Experimental evidence also exists for viral induction of cancer in clam 
haemocytes (Oprandy, 1981; Oprandy and Chang, 1983). A virus similar 
to a B-type retrovirus was isolated from M. arenaria and CCH were 
encountered following injection of this virus from CCH into healthy 
clams (Oprandy and Chang, 1983). ‘Virus-like particles’ have been iden¬ 
tified by electron microscopy (Elston etal., 1992) and by reverse transcrip¬ 
tase activity from CCH, but are absent from NCH, suggesting retroviral 
involvement in the disease (House et al., 1998; AboElkhair etal., 2009). 
Haemocyte cancer detected in M. arenaria in the field has been attributed 
to viral induction in isolated hotspots (Farley etal., 1986, 1991) and to 
species specificity of the disease (Kent etal., 1991; Muttray etal., 2009). In 
the laboratory, clam haemocyte cancer is most effectively transferred to 
normal clams after inoculation of normal clams with lysed CCH in 
diseased clam haemolymph, although haemolymph alone will also initiate 
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the disease in normal clams (Walker et al., 2009). Similar results have been 
obtained for Cerastoderma edule, following inoculation of normal cockles 
with haemolymph from diseased cockles. These results certainly suggest 
that a virus may be involved in the aetiology of bivalve haemocyte cancers 
(Collins and Mulcahy, 2003). However, irrefutable proof of the involve¬ 
ment of either environmental contaminants or viral agents or both in the 
aetiology of these diseases remains unavailable. 

These results point out that the aetiology of these diseases is complex 
and likely involves multiple causes including environmental, viral and 
genetic factors. Infectious and environmental aetiologies are not mutually 
exclusive and future research should focus on probing interactions 
between different disease factors. The use of haemocyte cell culture 
(Appendix A) and efforts to identify the source of haemocyte stem cells 
may be crucial in determining the aetiology of this disease. 



. Concluding Remarks 


Marine bivalves are ecologically and commercially important organ¬ 
isms that provide excellent models for both biomedical and environmen¬ 
tal researchers. For biomedical research, studies of naturally occurring 
cancers analogous to human diseases should indicate molecular mechan¬ 
isms that are held in common. As we have shown above, clam haemocyte 
cancer provides excellent in vitro and in vivo models for transcriptional 
and non-transcriptional outcomes reflecting those seen in similar human 
cancers under stress and with similar malfunctions in p53 functionality. 

Bivalves are simultaneously good models for the effects of toxic 
compounds, since they have low monooxygenase activity and thus, com¬ 
pounds slowly accumulate and tissue burdens can serve as integrative 
indices of bioavailable contaminants (Van Beneden, 2005). Ultimately, 
the naturally occurring, widely accessible and inexpensive-to-maintain 
clam haemocyte and germinoma cancer models may help to identify 
genes that are susceptible to environmental stress. Understanding the 
molecular basis of bivalve cancer and stress biology should significantly 
increase the usefulness of bivalve models as a sentinel species for sophis¬ 
ticated ecotoxicological studies and as environmental indicator species. 


6 .i. Bivalves as models for human cancer and stress biology 

Analysing the effects of genotoxic compounds on malignancy are cur¬ 
rently limited to highly regulated vertebrate cancer models where 
tumours are induced (e.g. the mouse). The CCH and germinoma models 
are more similar to outbreeding, human clinical populations than are 
those generated from inbred mouse strains intentionally exposed to 
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known tumour viruses (Kovalchuk etal., 2000). Fly and worm models for 
human cancers involving mortalin do not exist principally because adult 
somatic cells in these organisms do not divide, so an in vivo version of the 
model is unavailable. Also, existing fly and worm models have highly 
derived homologues for p53 only (p63 and p73 are lacking) (Brodsky 
et al., 2000; Jin et at., 2000; Derry etal., 2001; Schumacher et al., 2001). 
Additionally, very few members of the Bcl-2 family of proteins exist in 
either fly or worm models (Igaki and Miura, 2004; Brumby and 
Richardson, 2005; Lu and Abrams, 2006) and these genes are far less 
conserved than are those found in molluscs like the soft shell clam. 

Haemocyte cancer in clams is an example of a transplantable tumour 
that is not hampered by two of the characteristics that limit the applica¬ 
bility of mouse transplantable tumour models (Ostrand-Rosenberg, 
2004; Walker et al., 2009). For example, histocompatibility is not an issue 
since clams have only innate immunity (Medzhitov and Biron, 2003) and 
tumour cell inoculation is not subcutaneous as it is in mouse models, but 
occurs in an anatomically appropriate site for the tumour (the clam 
circulatory system) (Walker et al., 2009). Cancerous clam haemocytes 
provide excellent in vivo and in vitro models for human and other animal 
cancers displaying cytoplasmic sequestration of the p53 tumour suppres¬ 
sor by mortalin tethering. Clam haemocyte cancer is the only non¬ 
mammalian model sufficiently characterized to provide both in vitro and 
in vivo versions of the model from the same individuals for simultaneous 
use in studies of cytoplasmic sequestration and inactivation ofp53 proteins 
and their subsequent reactivation by treatments with genotoxic and non- 
genotoxic agents. 

6.2. In vitro and in vivo models for bivalve tumour cells 

The limited lifespan of human cells in culture has lead to the use of 
treatments that allow cells to achieve more passages. The lifespan 
of human cell lines can be extended by transfection with viral genes that 
produce products that sequester tumour suppressor molecules like p53 
and Rb. In these cultures some cells will occasionally acquire mutations 
that make them immortal and even tumorigenic (Masters, 2000). The 
intriguing feature of CCH is that although they do not necessarily appear 
to be immortalized, they are already tumorigenic. While CCFF main¬ 
tained in vitro are healthy and mitotic, refinements of our culture medium 
may be necessary to increase growth rates and to develop an immortalized 
cell line from CCH. However, if a virus is involved in this disease, the 
lifespan of CCH may already be extended via suppression of p53 tran¬ 
scription. Viral modification of the functions of CCH might be directly or 
indirectly involved in suppressing the activities of wild-type p53 in com¬ 
plexes with over-expressed mortalin that sequester p53 in the cytoplasm. 
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If this is the case, CCH may already be an immortalized cell line as defined 
for human cells in suspension culture (Masters, 2000; Walker etal., 2009). 

6.3. Use of bivalves as bioindicator species —understanding 
the signalling pathways involved 

Several recent studies of Mytilus sp. and other molluscs (e.g. the abalone, 
Haliotis tuberculata) have emphasized using changes in mRNA and protein 
levels for p53, heat-shock family member proteins (especially Hsp70 and 
Hsp90) and other cell stress-related gene products (e.g. SOD, superoxide 
dismutase and MT, metallothionein) as biosensors for environmental 
monitoring (Dondero etal., 2006a,b; Grattarola etal., 2006; Kourtidis 
et al., 2006a,b; Farcy et al., 2007, 2008; Banni et al., 2009; Pantzartzi 
et al., 2009, 2010). For example, many heat-shock family proteins are 
inducibly expressed in response to environmental stressors and function¬ 
ally integrate proliferative, survival and cell death signalling events (Beere, 
2001). In the case of clam haemocyte cancer, where the mitochondrial 
heat-shock protein mortalin is over-expressed, a response to stressors in 
the environment might imply a direct line from the environment to a 
well-studied molecular cancer mechanism that immortalizes molluscan cells 
(Kelley etal., 2001; Walker et al., 2006, 2009). Since mortalin is not heat- 
shock inducible, its over-expression must rely upon a different stressor. 

Environmental monitoring of appropriate gene expression profiles is 
potentially very useful to emphasize changes in the immediate environ¬ 
ment of bivalves and other members of marine communities. As a result, 
these issues can be addressed locally. However, the impact of potentially 
harmful environmental factors at the molecular level may be very com¬ 
plex. A recent paper by Pantzartzi et al. (2009) suggests that two cyto¬ 
plasmic Hsp90 genes from Mytilusgalloprovincialis contain promoters with 
p53 binding sites, suggesting the possibility that expression of these genes 
is transcriptionally regulated by p53. It has been shown that both Hsp70 
and Hsp90 participate with other proteins in multicomponent complexes 
that mediate translocation of cytoplasmic p53 to the nucleus using the 
molecular motor dynein (Galigniana et al., 2004; Giannakakou et al., 
2000, 2002). The response of bivalve cells to increases in Hsp90 resulting 
from stress-induced p53 transcription might drive changes in expression 
of mortalin, blocking apoptosis and immortalizing clam haemocytes. 
Even more intriguing is how stress can lead to mitochondrial transloca¬ 
tion of p53 and non-transcriptional induction of apoptosis in CCH. 
Certainly the dynein-dependent complex just described cannot alone 
be involved, since many mitochondria are distributed peripherally in 
the cytoplasm and are nowhere near the nucleus. While it is unknown 
even for human cancers what alternative multicomponent or other com¬ 
plexes involving mortalin and perhaps kinesin are involved in 
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mitochondrial translocation, recent studies have implicated both the 
DNAJ-like protein Tid 1 and loss of hyperubiquination in this process 
for several human colon carcinoma and breast cancer cell lines 
(Marchenko etal., 2007; Ahn etal., 2010; Kampinga and Craig, 2010). 
Despite our current understanding that changes in the expression of 
several stress-related genes can serve as excellent biosensors, the complex 
interactions in the normal signalling mechanism following a specific suite 
of environmental perturbations are virtually unknown. 
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Appendix A 


Walker’s Clam Haemocyte Culture Medium (Updated 
Formulation after Walker et al., 2009) 



Ingredient 

2 Liters 


1 Liter 


1. In Flask - stirring, add DI Water 

(ml) 


(ml) 

DI Water 

1000 

□ 

500 

2. Add Media 

(g) 


(g) 

Eagle's MEM (Sigma: M-3024) 

9.400 

□ 

4.700 

3. Add supplementary components 

(g) 


(g) 

CaCI2x2H20 

3.647 

— 

1.824 

KCI 

1.358 


0.679 

MgCI2x6H20 

8.714 


4.357 

MgS04x7H20 

6.316 


3.158 

NaCI 

48.520 


24.260 

NaHC03 

2.000 


1.000 

Hepes Buffer 

4.766 


2.383 

Glucose 

1.000 


0.500 

Glutamine 

0.292 

_ 

0.146 


(ml) 

□ 

(ml) 

Phenol Red (0.02% w/v) 

2.0 

1.0 

4. Add DI Water to bring the volumes up to: 

(ml) 


(ml) 

DI Water 

1740 

□ 

870 

5. pH to 7.2 (using 1M NaOH or HCI) | 


6. Add AntiBiotics and Nourishment 

(ml) 


(ml) 

Fetal Bovine Serum (FBS) [10% Final] 

200 


100 

Ampicillin (10 mg/ml) - 100 mg/L [Final] 

20 


10 

Penicillin/Streptomycin Mix (P4458-20ml)* 

40 


20 

^Penicillin (5000 U/ml) and Streptomycin (5 mg/ml) 




FINAL VOLUME: 

2000 


1000 


5. Vacuum-Filter using PES apparatus and then take the apparatus to the laminar 
flow hoods. Under sterile conditions, remove the top of the filter-apparatus, and 
cap/pour out aliquots into sterile containers. Store excess at 4 C (amino acids 

|should not freeze). 



_ 
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Abstract 

Stress regimes defined as the synchronous or sequential action of abiotic 
and biotic stresses determine the performance and distribution of species. 
The natural patterns of stress to which species are more or less well 
adapted have recently started to shift and alter under the influence of 
global change. This was the motivation to review our knowledge on the 
stress ecology of a benthic key player, the macroalgal genus Fucus. We first 
provide a comprehensive review of the genus as an ecological model includ¬ 
ing what is currently known about the major lineages of Fucus species with 
respect to hybridization, ecotypic differentiation and speciation; as well as 
life history, population structure and geographic distribution. We then 
review our current understanding of both extrinsic (abiotic/biotic) and 
intrinsic (genetic) stress(es) on Fucus species and how they interact with 
each other. 

It is concluded that (i) interactive stress effects appear to be equally 
distributed over additive, antagonistic and synergistic categories at the 
level of single experiments, but are predominantly additive when averaged 
over all studies in a meta-analysis of 41 experiments; (ii) juvenile and adult 
responses to stress frequently differ and (iii) several species or particular 
populations of Fucus may be relatively unaffected by climate change as a 
consequence of pre-adapted ecotypes that collectively express wide phys¬ 
iological tolerences. 

Future research on Fucus should (i) include additional species, (ii) 
include marginal populations as models for responses to environmental 
stress; (iii) assess a wider range of stress combinations, including their 
temporal fluctuations; (iv) better differentiate between stress sensitivity of 
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juvenile versus adult stages; (v) include a functional genomic component in 
order to better integrate Fucus’ ecological and evolutionary responses to 
stress regimes and (vi) utilize a multivariate modelling approach in order to 
develop and understand interaction networks. 


1. Introduction 

The shifting of environmental variables in the course of ongoing 
global climate change is expected to impact the performance, and ulti¬ 
mately the distribution, of numerous species in marine coastal systems 
(reviewed in Harley et al., 2006; Parmesan, 2006; IPCC Climate Change, 
2007). Many shallow coastal habitats, however, are defined by a physically 
demanding environment with steep abiotic gradients and drastic environ¬ 
mental fluctuations at small spatial and temporal scales. Marine organisms 
living in these intertidal or shallow subtidal habitats are regularly exposed 
to strong water motion and subjected to extreme fluctuations in temper¬ 
ature, pH, irradiance, salinity or nutrient availability, and the amplitude of 
these fluctuations far exceeds climate changes predicted for the coming 
decades (e.g. Thomsen and Melzner, 2010). Although organisms in these 
habitats cope with ambient abiotic stresses at least to the point of transient 
tolerance, they must also contend with stressful biotic interactions 
including competition, epibiosis, parasitism and herbivory, all of which 
have the potential to modulate the abiotic stresses (e.g. Wahl, 2008b). 
Understanding how organisms in harsh and fluctuating habitats cope with 
single and multiple stresses is essential to clarifying and evaluating the risks 
of global change. 

Foundation species in the intertidal/shallow subtidal regions of north¬ 
ern hemisphere temperate coasts frequently include members of the algal 
genus Fucus, which typically consists of three to four zoned species across 
the intertidal-shallow subtidal gradient (Fig. 2.1). Their ecology, physi¬ 
ology and genetics have been the subjects of intense research during 
the past four decades, although some species have been studied since 
the mid-1800s. Thus, results of these many studies provide insights into 
stress ecology and may form the basis for understanding how global 
climate change will affect northern intertidal/subtidal habitats. 

i.i. The concept of stress 

Stress is defined here as the impact of any set of abiotic and/or biotic 
factors that adversely affects individual ‘performance’ and ultimately 
impairs population growth rate through reduced individual survival, 
growth and/or reproduction (Grime, 1989; Vinebrooke et al., 2004). 
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(a) 



Figure 2.1 Zonation patterns characteristic of Fucus species along (a) intertidal and 
(b) atidal shores of the Baltic Sea. 

In this sense, stress is ubiquitous and widespread, particularly at the 
margins of species’ distributions (outside of which population growth 
remains negative) and at isolated locations within a distributional range. 

Stress remains a contentious concept despite its widely recognized 
importance. Whether or not a factor is considered stressful depends on 
the target organism; the intensity, duration and recurrence of the stress; 
and on the various interactions among stresses. As stresses typically co¬ 
occur, it is crucial to determine whether a set of stresses is likely to act 
additively, antagonistically or synergistically, which in turn, will deter¬ 
mine whether the combined impact will be minor or detrimental. 
Understanding interactions and feedbacks is one of the fundamental 
challenges for understanding ecological dynamics at a variety of scales 
(Green and Sadedin, 2005; Christensen et al., 2006; Gamfeldt and 
Kallstrom, 2007) and is a necessary prerequisite for any prediction of 
responses to changing stress regimes. Nevertheless, most studies of stress 
have focused on the physiological responses of organisms to an unfavour¬ 
able variable or set of variables (Crain et al ., 2008; Darling and Cote, 2008; 
Schiel, 2009) without linking those responses to the biotic interactions 
among the species being stressed or to the underlying phenotypes (Davis 
et al., 2005; Hughes et al., 2008). 

Stresses act as selective agents for stress tolerance. Adaptability of a 
population therefore is determined by the interaction between the extrin¬ 
sic components of stress—abiotic and biotic factors—and the intrinsic 
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component of stress (Bijlsma and Loeschcke, 2005) estimated as the 
standing genetic variation. If this evolutionary potential is compromised 
(i.e. in small, isolated populations or at range margins), then an extrinsic 
stress will be further magnified. It is now recognized that the temporal 
differential between ecological and genetic changes can be small or neg¬ 
ligible (Spielman et al, 2004). Thus, stress regimes can have a significant 
and nearly simultaneous impact in both ecological and evolutionary 
time scales. Understanding this coupling is the major challenge of eco¬ 
logical and evolutionary genomics (EEG) in natural communities (Tautz 
et al, 2010). 

1.2. Background studies and the constraints 
of monodisciplinary research 

General reviews of intertidal/shallow subtidal ecology include Benedetti- 
Cecchi (2000), Connell (1972), Gruner et al. (2008), Menge (1995); for 
intertidal physiological ecology, see Helmuth et al. (2002); for intertidal 
community genetics, see Schmidt et al. (2008) and for environmental 
change impact on marine ecosystems, see Harley (2006), Helmuth 
et al. (2006), Hillebrand etal. (2010), Korpinen etal. (2007a), Parmesan 
(2006), Schiel et al. (2004) and Thompson et al. (2002). Although fucoids 
and other macrophytic algae are sometimes included in these reviews, 
most reviews are overwhelmingly focused on sessile invertebrates (but see 
Davison and Pearson, 1996). 

Research on the genus Fucus has produced >1500 publications 
(excluding purely taxonomic contributions) over the past decades. 
Although the genus provides a unique system for examining individual, 
population, species and community level processes, these hierarchical 
levels have fostered persistent splits in research approaches (see review 
by Tomanek and Helmuth, 2002). For example, physiological ecologists 
are interested in testing thermal limits, photosynthetic capacity, desicca¬ 
tion tolerance and osmoregulation of species as the basis for upper- 
intertidal zonation or geographical range limits. Reviews of physiological 
stress and photobiology in algae sensulato include Bischof et al. (2002), 
Chapman (1995), Davison and Pearson (1996), Dethier and Williams 
(2009), Dethier et al. (2005), Schiel and Foster (2006), Wiencke et al. 
(1992, 2006). Community ecologists, however, focus on biotic interac¬ 
tions related to competition, epibiosis, predation and defence strategies 
against the former, as well as on phenology, reproduction, recruitment 
and life-history traits as the key factors shaping populations and their 
inter-specific interactions in a community. Specific reviews have been 
published on the ecology of brown seaweeds in general (Schiel and Foster, 
2006), on Laminariales and Fucales (Bartsch et al, 2008), and on the 
Fucales (Chapman, 1995). Last but not least, molecular ecologists 
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(who entered the field in the 1990s) focus on the role of genetic variation, 
fitness, gene flow and selection as the determinant processes leading to 
adaptation. Consequently, a large body of research has accumulated on the 
roles of various types of stress on Fucus spp. but relatively little integration 
or synthesis has emerged. 

Over the past 15 years physiological and community ecology studies 
have greatly expanded and become more rigorous, and when combined 
with advancing molecular approaches, have the potential to elucidate the 
complex pathways and networks of stress effects. Chapman (1995) 
focused his review of the genus Fucus on the years 1970—1993; our review 
focuses on subsequent studies. 

We begin our review with an overview of Fucus and its ecological 
importance. We then turn to a review of studies on abiotic stresses, both 
single and multiple, and their reported impacts at different ontogenetic 
stages of Fucus life history. Next we introduce biotic interactions, how they 
can be shifted by abiotic stresses, and how these shifts may increase or buffer 
stress effects. Woven throughout this part of the review we also provide the 
results of a meta-analysis of stress combination studies. We then turn to 
the genetic level and review what is known about genetic stress and the 
potential for adaptation in the new perspective of ecological and evolu¬ 
tionary genomics. In the final section, we discuss our conclusions and 
prospects, aiming to inspire the research agenda for the coming years. 


2. Methods 


In addition to a review of individual stress studies, a meta-analysis of 
multi-stress experiments was performed in an attempt to elucidate general 
effects of stress combinations on Fucus. 


2.1. Data search 

We searched for articles in Science Citation Index Expanded (ISI Web of 
Science) that conducted factorial experiments with at least two stresses 
using the following search parameters: (1) ts = (stress* or salinity or tem¬ 
perature or UV or irradiation or irradiance or desiccation or emersion or 
pollution or eutrophication or ‘nutrient enrichment’ or disturbance or 
osmoregulation or ‘wave exposure’ or ‘wave energy’ or ‘water motion’ 
or ‘wave action’ or ‘ice scour’ or competition or herbivor* or grazing or 
epibiot* or epiphyt* or fouling or pathogen* or parasit*), (2) ts = (inter¬ 
act* or antagonis* or synerg* or additiv* or non-additiv* or multiple or 
combined or factorial or experiment*), (3) ts = Fucus or fucoid* or 
fucales or fucaceae. The search yielded a total of 488 papers, but only 41 
experiments with fully crossed factorial experiments, adequate replication 
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and response to stress (measured as a quantitative performance variable) 
were chosen for our analysis (Table 2.1). 

2.2. Data acquisition 

Figures were scanned from the original articles and the data (means and 
standard deviations for the different stress treatments) extracted using 
GrabIT v. 1.7.2 XP add-in for Microsoft Excel (Datatrend Software, 
USA). Some papers included multiple experiments from which data 
could be extracted yielding a final meta-analysis data set with 41 entries. 
The full data set is available in spreadsheet format in Appendix SI. 

Note that for biological stresses (which often use an exclusion treat¬ 
ment to examine the effects of herbivory and competition), the exclusion 
treatment without the biological stress becomes the stress-free control in 
our meta-analysis. 

2.3. Groupings 

All groupings represent trade-offs, especially with respect to the number 
of studies available. Given the relatively small number of studies meeting 
our criteria and the large number of stresses tested, we opted for broad 
categories (Table 2.2). 

Data were divided into subsets that included two-way interactions 
between treatments (i.e. abiotic X abiotic, biotic X biotic, abiotic X bio¬ 
tic, abiotic X nutrient enrichment and biotic X nutrient enrichment). 
Three-way (or higher) interaction experiments were divided into separate 
two-way treatment combinations. In some cases, a single experiment 
could be divided and used in several analyses, thus the number of experi¬ 
ments used is larger than the number of articles listed in Table 2.1 and 
the number of comparisons used in any one grouping differed among 
tests depending on the combinations. For example, abiotic stress X 
nutrient enrichment and abiotic stress X biotic stress involve 15 and 
12 comparisons, respectively. The interaction between nutrient enrich¬ 
ment 1 X nutrient enrichment 2 was not analysed because only two 
experiments were found. 


2.4. Effect sizes 

Calculation of effect sizes followed Hedges et al. (1999). For each subset of 
the data we calculated the log response (effect) ratio (LRR), defined as the 
dimensionless ratio of the treatment over the control response value 
(effect = ln(stress treatment/control)). Thus, a LRR of 0.69 corresponds 
to a 100% increase in performance of the response variable, and a LRR of 
—0.69 corresponds to a 50% decrease. Each treatment and combination 




Table 2.1 Summary of studies exploring interactions including those articles used for the meta-analysis 




Interaction 

Focal species 

Performance 

Quality of the 

Method* 5 

Reference 

Used for 



trait 

interaction 3 



meta¬ 







analysis 

Abiotic—abiotic 







Nutrient 

F. serratus. 

Growth 

X: no interaction of 

M 

Kraufvelin 

Y 

enrichment X 

F. vesiculosus 


nutrient 


(2006) 


wave action 



enrichment and 







wave action 




UVR X salinity 

F. vesiculosus 

Photosynthesis 

An: low-salinity stress 
decreased the 

M 

Nygard and 
Ekelund 

Y 




sensitivity of F. 
vesiculosus towards 


(2006) 





UV 




UVR X 

F. gardneri 

Germination 

An: negative effects of 

L 

Hoffman etal. 

N 

temperature 


rate, germling 
cell number 

UVR on 
germination are 
reduced under 
high temperatures. 
Ad: negative effects 
of UVR and low 


(2003) 





temperature on cell 
division 




Light X 

F. vesiculosus 

Microbial 

An: low light 

L 

Wahl et al. 

N 

temperature 


fouling 

counteracts the 
microfouling— 
enhancing impact 
of warming 


(2010) 



Martin Wahl et al. 




Light X 
desiccation 

F. serratus 

Photoinhibition 

An: desiccation 
reduces the 
photoinhiting 
effects of high 
irradiation 

F 

Salinity X 
temperature 

F. vesiculosus 

Germination 

S: hyposalinity impact 
was enhanced by 
increasing 
temperature 

L 

Salinity X 
nutrients 

F. vesiculosus 

Growth 

Ad: low nutrients 
enhanced the effect 
of low salinity 

L 

Desiccation X 

wave action 

F. distichus 

Survival adults 

S c : desiccated thalli 
experienced a 
higher wave- 
induced mortality 

L, F 

Nitrate 
enrichment X 
phosphate 
enrichment 

F. vesiculosus 

Germination 

S: negative effect of 
phosphates 
enhanced by 
nitrate addition 

L 


F. vesiculosus 

Biomass 

production 

Ad: shading decreases 
biomass 

F 


Huppertz etal. N 
(1990) 


Maczassek N 

and Wahl 
(in prep.) 

Nygard and N 

Dring (2008) 

Haring etal. N 

( 2002 ) 


Bergstrom etal. Y 
(2003) 


Eriksson etal. Y 
(2006) 
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Interaction 

Focal species 

Performance 

Quality of the 

Method* 3 

Reference 

Used for 



trait 

interaction 3 



meta¬ 







analysis 

Nutrient 



production, 




enrichment X 



nutrients have no 




shade 



consistent effects 
but tend to vary 
among habitats of 
different 
complexity 




Parental 

F. vesiculosus 

Survival, 

S: parental exposure 

L 

Li and Brawley 

Y 

temperature X 


growth 

to increased 


(2004) 


embryonic 



temperature 




temperature 



enhanced heat- 
shock resistance of 
embryos 




Salinity X 

F. vesiculosu. 

Growth of 

Ad-An d : hypo 

L 

Munda and 

Y 

temperature 

F. distichus. 

germlings 

salinities and 


Kremer 



F. virusoides 


unfavourable 
temperature (high 
and low) affect 
germlings 
negatively and the 
combined effect is 


(1977) 



Martin Wahl et at. 





Salinity X F. vesiculosus, F. Survival of 

temperature spiralis embryos 


Biotic—biotic 

Competition X F. evanescens Growth rate 

grazing 

Competition X F. vesiculosus Colonization 

grazing 


success 


L Li and Brawley Y 

(2004) 


typically less than 
the sum of single 
stressor effects 

An: hypersalinity 
minimized 
negative heat- 
shock effects on 
survival 

Ad: competition by F 

red algal canopy 
and grazing impair 
growth rate 

Ad to (S) depending F 

on grazer 
community. 

Negative effect of 
competition tends 
to be accentuated 
by gastropod 
grazing 


Worm and Y 

Chapman 
(1998) 

Korpinen and Y 

Jormalainen 
(2008a) 
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Interaction 

Focal species 

Performance 

trait 

Epibiotism X 

F. vesiculosus 

Growth rate 

grazing 

Consumption X 

F. disticluis 

Endophytism 

parasitism 

Grazing X 
grazing 

F. vesiculosus 

Defence 

induction 


Competition X 
grazing 


F. vesiculosus, 
F. serratus 


Cover 


-1^ 

00 


Quality of the 
interaction 3 

Method* 5 

Reference 

Used fo 

meta¬ 

analysis 

Ad: epibiotism and 
grazing reduce 
growth 

M 

Jormalainen 
etal. (2008b) 

Y 

An: grazing by 
littorinids and 
gamarids reduces 
endophyte 
infection 

F 

Parker and 
Chapman 
(1994) 

N 

An: grazing by one 
consumer species 
reduces the impact 
of a second 
consumer species 

L 

Yun etal. (2010), 
Long etal. 
(2007) 

N 

Ad: in F. serratus, 
competition by 

F 

Jenkins etal. 
(1999) 

Y 


other canopy¬ 
forming algae 
decreased cover, 
limpet grazing had 
no effect. (S): in F. 
vesiculosus, grazing 
tended to amplify 
the negative effect 
of competition 


Martin Wahl et al. 





Competition X F. serratus, Cover, 

grazing F. vesiculosus recruitment 


Abiotic—biotic 

Warming X F. vesiculosus Cover 

competition 


Emersion X F. vesiculosus Fouling 

epibiotism resistance 


Ad: competition by F Cervin etal. N 

canopy algae and (2005) 

grazing reduce 
recruitment, 
sometimes to a 
degree that the 
interactive effect 
becomes An 
because there is no 
room for additive 
decrease as cover 
hits zero 

(Ad): heat-shock F Allison (2004) Y 

decreases cover 
development, no 
effect of inter¬ 
specific 
competition 

An c : emersion L Brock etal. N 

reduces fouling (2007) 

(by enhancing 
anti-fouling 
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Interaction 

Focal species 

Performance 

Quality of the 

Method* 3 

Reference 

Used for 



trait 

interaction 3 



meta¬ 







analysis 




metabolite 

concentration) 




Warming X 

F. vesiculosus 

Growth 

An c : warming reduces 

F 

Moore etal. 

N 

grazing 



grazing by driving 
a shift from a 
voracious to the 
less voracious 


(2007) 





consumer 




Wave exposure X 

F. vesiculosus 

Cover- 

rd 

o : exposure 

F 

Jonsson etal. 

Y 

grazing 



promotes negative 
effect of grazing 


(2006) 


Wave exposure X 

F. vesiculosus 

Vulnerability to 

An: exposure reduces 

F 

Nietsch (2009) 

N 

grazing 


grazing 

isopod grazing by 
toughening of 
thallus 




Increasing depth 

F. serratus 

Abundance 

An c : competition 

C 

Malm and 

N 

X competition 



with F. vesiculosus 


Kautsky 





relaxes in deeper 
depth 


(2003) 


Irradiance X 

F. vesiculosus 

Palatability 

An: light limitation 

L 

Weinberger 

N 

grazing 



decreases mannitol 
content and thus 
attractiveness of 


et al. (2011) 



Martin Wahl et al. 





Increasing 
depth e X 

F. vesiculosus 

Growth rate 

grazing 

Increasing 

F. vesiculosus 

Growth rate 

depth e X 
epibiotism 
Irradiance X 

F. vesiculosus 

Anti-grazer 

grazing 


defence 

Desiccation X 

F. distichus 

Growth, 

grazing 


reproduction. 

survival 

Emersion X 

F. distichus 

Growth, 

grazing 


reproduction 

Wave exposure X 

F. vesiculosus 

Grazing loss 


grazing X 
competition 


Fuchs to isopod 
grazers 


S: negative effect of 

F 

Jormalainen and 

N 

grazing increases 


Ramsay 


with deeper depth 


(2009) 


Ad: epibiotism and 

F 

Rohde etal. 

Y 

increasing depth 
reduce growth rate 


(2008) 


Ad: light limitation 

L 

Rohde etal. 

Y 

decreases anti¬ 
grazer defences 


(2004) 


Ad: desiccation and 

F, M 

Dethier etal. 

Y 

grazing impair 
growth and 
reproduction 


(2005) 


An c : emersion 

C 

Dethier and 

N 

reduces grazing 


Williams 


pressure (but this 
was dependent on 
season) 


(2009) 


S: grazing pressure 

M 

Engkvist etal. 

N 

shifts from 

F. serratus to 

F. vesiculosus under 
increased exposure 


(2004) 
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NJ 


Interaction 

Focal species 

Performance 

trait 

Quality of the 
interaction 3 

Method* 5 

Reference 

Used for 

meta¬ 

analysis 

Wave exposure X 
grazing 

F. vesiculosus, F. 
serratus 

Grazing loss 

(Ad): negative effect 
of grazing but no 
effect of exposure 
on overall grazing 
loss 

M 

Engkvist etal. 
(2004) 

N 

Nutrient 
enrichment X 
grazing 

F. vesiculosus 

Colonization 

success 

(S to An): depending 
on the grazer 
community 

F 

Korpinen and 
Jormalainen 
(2008a) 

N 

Nutrient 
enrichment X 
grazing 

F. vesiculosus 

Colonization 

success 

An d : grazers 
counteract 
negative (indirect) 
effects of nutrient 
enrichment 

F 

Worm et al. 
(2001) 

N 

Nutrient 

enrichment X 
epibiotism 

F. vesiculosus 

Growth rate 

S: nutrient 
enrichment 
enhances negative 
effect of epibiotism 

M 

Jormalainen 
etal. (2003) 

Y 

Nutrient 
enrichment X 
grazing 

F. vesiculosus 

Growth rate 

An: gastropod grazing 
counteracts the 
negative effect (due 
to increased 
epibiotism) of 

M 

Raberg and 
Kautsky 
(2008) 

Y 


Martin Wahl ef at. 





nutrient 

enrichment 


Nutrient 
enrichment X 
inter-specific 
competition 

F. vesiculosus 

Colonization 

success 

S: nutrient 
enrichment 
enhanced the 
negative effect of 
competition 

F 

Korpinen and 
Jormalainen 
(2008a) 

Y 

Nutrient 

enrichment X 
temperature X 
inter-specific 
competition 

F. serratus, 

F. evanescens 

Germling 

growth, 

germling 

survival 

S: warming and 
nutrient 
enrichment 
enhance 

competition effects 

L 

Steen (2004) 

N 

Nutrient 
enrichment X 
temperature X 
intra-specific 
competition 

F. serratus, 

F. evanescens 

Germling 

growth, 

germling 

survival 

Ad: temperature 
increased 
competition and 
reduced growth. S: 
in F. serratus, 

L 

Steen and 
Scrosati 
(2004) 

Y 


negative effect of 
competition on 


survival increased 
with warming and 

on 

U) 
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Interaction 

Focal species 

Performance 

Quality of the 

Method* 3 

Reference 

Used for 



trait 

interaction 3 



meta¬ 







analysis 




nutrient 







concentration 




Season X 

F. vesiculosus 

Germination 

S: eutrophication 

F 

Berger etal. 

N 

eutrophication X 



favours competing 


(2004) 


competition 



annual algae more 
in spring than in 







autumn 




Nutrient 

F. vesiculosus 

Colonization 

An: negative grazing 

F 

Korpinen etal. 

Y 

enrichment X 


success 

effect is dampened 


(2007b) 


grazing 



by negative effect 
of nutrient 







enrichment 




Nutrient 

F. vesiculosus 

Number of 

Ad: low light, 

F 

Eriksson etal. 

Y 

enrichment X 


recruits 

nutrient 


(2007) 


grazing X 



enrichment and 




shade 



grazing reduce 







recruitment 




Nutrient 

F. vesiculosus 

Palatability 

S c : nutrient 

L 

Henimi and 

N 

enrichment X 



enrichment 


Jormalainen 


grazing 



enhances grazing 


(2002) 



Martin Wahl et at. 





Nutrient 

F. uesiculosus 

Growth (apical 

enrichment X 


tip divisions) 

grazing 



Nutrient 
enrichment X 

F. serratus 

Palatability 

grazing 

Nutrient 

F. uesiculosus 

Biomass 

enrichment X 
competition X 
grazing 




(S): nutrient 

M 

Hemnii etal. 

Y 

enrichment tended 


(2005) 


to increase grazing 
losses 




S c : nutrient 

M 

Kraufvelin etal. 

N 

enrichment 
enhances grazing 


(2006) 


An: nutrients enhance 

F 

Worm et al. 

Y 

competition 
(when grazers are 
present) 


(2001) 



a Ad, additive (sum = A + B), no significant interaction; S, synergetic (sum > A + B), positive interaction; An, antagonistic (sum < A + B), negative interaction; Q, 
indicate a non-significant tendency (0.05 < P < 0.1) and X, no interaction. 
b L, lab experiment; M, flow-through outdoor mesocosm experiment; F, field experiment and C, correlative field data. 
c Interactive effect on algal performance interpreted from non-factorial designs where algal performance was not directly measured. 
d Interaction suggested, but statistical test lacking. 

e Light availability was the main abiotic factor (decreasing with increasing depth). 
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Table 2.2 Stress groupings used in the meta-analysis 


Abiotic stress 

Depth 


Desiccation 


Irradiance 


Salinity 


Shading 


Temperature 


Wave action 

Biotic stress 

Competition 


Herbivory 

Nutrient enrichment 

Nitrogen and phosphorous enrichment 


was compared separately against a ‘stress free’ control (Elser et al ., 2007). 
For each separate experiment we calculated the natural LRRs using the 
following equations: 


LRRmain effect stress A — (In JVLecUlAlBo) 

(In Mean A oBo) 

(1) 

LR-R-main effect stress B (In JVleailAOBl) 

(In Mean AO Bo) 

(2) 


L^R^Rcombined effect of stress A and B = (lnMean A iBi) - (lnMean A oBo), (3) 

where A and B refer to the different stresses in the two-way interaction. 
Al and B1 refer to the treatments with the stresses added, and AO and BO 
refer to the unstressed controls. In order to account for differences in 
replication and precision associated with individual studies, we calculated 
average effects of experimental treatments as weighted averages of the 
natural log response ratios (LRRs) obtained in Equations (1)—(3) (LRR*; 
Hedges et al, 1999). Statistical significance was tested by calculating the 
confidence interval around the weighted average (Hedges et al ., 1999). An 
effect was considered significant if the confidence interval did not cross the 
zero line. 

Given two stresses, an additive interaction is the sum of the two 
stresses. While the basic analysis provides an indication as to whether 
there is a synergistic or antagonistic trend, as compared to a simple 
additive combined effect, it does not explicitly test for an interaction. 
Therefore, we conducted a second analysis using the equations developed 
in Gruner et al. (2008). Note that in this analysis we used unweighted 
effect sizes. In this approach, an interaction effect is additive if the inter¬ 
action LRR does not differ from zero (i.e. the confidence interval does 
cross the zero line) and thereby generates a true statistical test of whether 
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the interaction effect differs from additivity. Given two stresses, synergistic 
interactions are indicated by a negative value and antagonistic interactions 
by positive interaction LRRs. 


3 . The Genus Fucus 

We describe stress ecology using Fucus because (1) the genus includes 
a set of ecologically important foundation species with a wide longitudi¬ 
nal and latitudinal distribution, (2) most Fucus species occur in stressful 
habitats of the intertidal or shallow subtidal and (3) some species (notably 
Fucus vesiculosus ) are already experiencing losses and gains in their geo¬ 
graphic distributions. 

3.1. Global distribution and diversification 

Fucus spp. commonly occur along protected and exposed rocky intertidal 
and shallow subtidal shores (Fig. 2.1), as well as in tidal marshes, over a 
wide latitudinal gradient within the North Pacific Ocean, North Atlantic 
Ocean and Baltic Sea (Liming, 1990). The genus is rare in the 
Mediterranean, where a single species is found only in the northern 
Adriatic Sea. Only one species, the cold-adapted Fucus distichus (see 
Coyer et al , 2006a for taxonomic notation used in this review) is found 
in both the Pacific and Atlantic oceans, ranging from Japan and Alaska in 
the North Pacific to Svalbard, northern Norway, the White Sea, Iceland 
and northern Scotland in the North Atlantic (with recent introductions to 
the North and Kattegat Sea and Bergen Harbor (Wikstrom et al, 2006; 
Sj0tun, personal communication). Other species are found only in the 
North Atlantic (the Pacific occurrence of Fucus spiralis likely is an intro¬ 
duction; Coyer et al, 2011a and references therein), where they range 
from northern Canada and northern Norway/White Sea to the Iberian 
peninsula, which is the southern limit for many but not all species. For 
example, while the distributions of Fucus serratus/high shore F. spiralis and 
F. vesiculosus end in northern and southern Portugal, respectively, the 
southern form of F. spiralis occurs in southern Portugal, the Canary 
Islands, the Azores and northern Morocco (Coyer et al, 2003; Billard 
et al, 2010; Zardi et al, unpubl.). 

In the North Atlantic, most Fucus species are characterized by hotspots 
of genetic and species diversity in southwestern Ireland and the Brittany 
peninsula of France or NW Iberia, all putative glacial refugia during the 
Last Glacial Maximum (Coyer et al, 2003; Hoarau et al, 2007; Coyer 
et al, 2011c; Neiva et al, 2010). On the other hand, genetic diversity is 
markedly reduced in populations of F. vesiculosus along the North 
American Atlantic coast a likely signature of postglacial recolonization 
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from Europe (Muhlin and Brawley, 2009) and in F. serratus at the edges of 
its southern distribution on the Iberian Peninsula (Coyer et al, 2003). 

Rapid ecotypic diversification is reflected in patterns of zonation along 
the shore, as well as morphological variability, hybridization and poly¬ 
ploidy. The ecotypic attributes have resulted in >150 described species, 
subspecies and forms (www.AlgaeBase.org, Guiry and Guiry, 2010), 
although recent molecular phylogenetic studies suggest the number of 
taxa to be <10 (Serrao et al, 1999a; Coyer et al, 2006a; F. Canovas, C. 
Mota, E. Serrao and G. Pearson, unpubl. data). 

The discovery of sibling and/or cryptic species by molecular methods 
is relevant to both ecological and evolutionary questions regarding stress, 
because smaller bounds need to be established on the distributions of some 
species. For example, recent work has demonstrated conclusively that 
Fucusradicans is not an ecotypic variant of F. vesiculosus (Bergstrom et al, 
2005), but is a separate species exclusively confined to the northeastern 
and northwestern Baltic and has evolved in the past 400—2000 years 
(Pereyra et al, 2009). Other recent studies have revealed that the high- 
and low-shore form of F. spiralis are reproductively isolated (Billard et al, 
2010) and allopatric F. spiralis found in southern Portugal and 
North Africa constitutes a new species (Coyer et al, 2011a; Zardi et al, 
unpubl.). 

Hybridization and introgression, which further complicate our 
understanding of the species’ distribution, are well documented between 
F. vesiculosus and F. spiralis (Wallace et al, 2004; Billard et al, 2005a,b, 
2010; Engel et al, 2005; Coyer et al, 2011a), as well as between Fucus cer- 
anoides and F. vesiculosus (Neiva et al, 2010), and F. serratus and F. distichus 
(Coyer et al, 2002a,b, 2007). Some hybrids (F. vesiculosus X F. spiralis) or 
polyploidy variants (F. vesiculosus, Wallace et al., 2004; Coyer et al., 2006b) 
commonly occur, while other hybrids are less fit ( F. serratus X F. distichus) 
(Coyer et al, 2007). 

In conclusion, a confusing taxonomy combined with complex pro¬ 
cesses of hybridization and rapid ecotypic differentiation make it advisable 
to genetically identify Fucus entities at the outset of any proposed study. 
With few exceptions (e.g. F. serratus, high shore F. spiralis, F. vesiculosus with 
vesicles), traditional (visual) taxonomy is unreliable. 


3.2. Life history and demography 

Fucoid life histories are monophasic with a diploid adult (van den Hoek 
et al, 1995). Haploid sperm and eggs are produced in numerous con- 
ceptacles within each of several large receptacles on apical tips of thallus 
branches. Reproduction is iteroparous, with both dioecious and her¬ 
maphroditic species. Eggs are fertilized close to (or on) the female and 
the diploid zygotes typically settle and mature within meters of the female 
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parent (Serrao et al., 1997; Pearson and Serrao, 2006; Schiel and Foster, 
2006, and references therein). Limited dispersal of gametes, together with 
synchronous production of eggs and sperm, increases the probability of 
selfing and inbreeding in hermaphroditic species (Coleman and Brawley, 
2005; Perrin et al, 2007) and possibly breeding among closely related 
parents in dioecious species. 

Fertilization is dependent on both water motion and temperature. 
Release of gametes typically occurs in F. vesiculosus at low turbulence 
and during daylight (Pearson and Serrao, 2006). Calm conditions 
are required to ensure fertilization success and facilitate recruitment of 
the non-dispersive zygotes within the adult habitat. Active photosynthesis 
is needed to detect turbulence (Serrao et al., 1996b); in the absence of 
turbulence, photosynthesis depletes dissolved inorganic carbon in the 
boundary layer covering the thallus and this depletion, in addition to a 
general semilunar periodicity, is used as one cue for synchronous gamete 
release (Pearson and Brawley, 1998; Pearson et al., 1998). Fertilization 
success typically approaches 100% (Berndt et al., 2002). Temporal patterns 
of reproduction may vary even at the local scale. For example, the repro¬ 
ductive period varies between 4 and 12 months among F. vesiculosus 
populations separated by 100 km or less along the German Baltic coast 
(K. Maczassek, personal communication). Along the Swedish South 
and East coast summer- and autumn-reproducing forms of F. vesiculosus 
are heterogeneously distributed with a tendency for a dominance of 
autumn breeders on the central mainland coast and a dominance of 
summer breeders along the coasts of islands (Oland and Gotland) 
(Berger et al, 2001). Autumn recruits are expected to suffer less from 
sedimentation and competing ephemeral green algae than spring recruits 
(Berger et al., 2004). Although such differentiation of reproductive phe¬ 
nology may be expected to create temporal isolation and promote genetic 
differentiation, no genetic differentiation (using five microsatellites) was 
observed between the summer and autumn reproducing forms 
(Tatarenkov et al. , 2007). 

Mortality in fucoids is typically very high during the microscopic 
stages (Schiel and Foster, 2006). Estimates of survival from an egg to a 
benthic settler vary from 1.5 to 10%, and from a settler to a visible 
germling from 5 to 12% (Chapman, 1995). As thalli grow, mortality 
decreases but remains temporally variable (Wright et al, 2004). 
Egg production in the well-studied F. distichus is about 1 X 10 5 — 
2 X 10 6 eggs m -2 mo 1 (Ang, 1991). Of these, only a few hundred 
per m~“ (i.e. less than 0.015%) succeed in the transition to a visible 
juvenile (Fig. 2.2). Fucus spp. typically reach an age of 5—10 years 
(Chapman, 1995). Due to high mortality of early life-history stages, there 
is often a paucity ofjuveniles in the mature stands, rendering populations 
vulnerable to catastrophic mortality (Dudgeon and Petraitis, 2005). 
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Space competition 
Light competition 
Low-Light stress 
High-Light stress 
UV stress 
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Wave exposure 
Ice disturbance 
Epibiosis 
Grazing 


Figure 2.2 (a) Typical survivorship curve of Fucus species (expressed as surviving 

individuals N m “) through sequential developmental stages and (b) an indication of 
the relative sensitivities of the various stages to various stresses. The survivorship 
curve is an approximation based on data from F. disticlms (Ang, 1991). 
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However, two recruitment strategies help counter the high early mortal¬ 
ity and contribute to population replenishment. First, the microscopic 
stages may exhibit delayed growth of several months, which creates a 
propagule bank under the canopy (Schiel and Foster, 2006). A second 
strategy involves fragmentation (Tatarenkov et al. , 2005) and regeneration 
(Schiel and Foster, 2006), which allows population re-growth after cat¬ 
astrophic mortality. Asexual reproduction by reattachment of fragments 
has been documented at high frequency only in F. radicans (Tatarenkov 
et al, 2005). Regeneration of fronds from attached adult holdfasts is 
common in F. vesicidosus (Kiirikki and Ruuskanen, 1996) and allows 
population recovery after mechanical damage, that is, by ice scour. 

Dispersal is highly variable. As distinct rocky intertidal habitats are 
often separated by soft bottom areas unsuitable for Fucus attachment, 
reproduction takes place mostly among the local individuals (Coyer 
et al., 2003). Direct estimates of average dispersal distances of fucoid 
zygotes are only a few meters, rarely exceeding 10 m (Dudgeon and 
Petraitis, 2001). Indirect estimates based on spatial auto-correlation and 
pairwise F sx (0) analysis suggested a panmictic unit (area of random 
reproduction) of 0.5—2 km (see Section 5) (Coyer et al, 2003). Long¬ 
distance dispersal and gene flow are possible by rafting ofintact individuals 
or detached pieces of receptacle-bearing thalli (Viejo and Arrontes, 1992; 
Thiel and Gutow, 2005). 


3.3. Latitudinal and zonational shifts 

Fucus spp. have high photosynthetic rates, surpassing phytoplankton spe¬ 
cies on a per-unit, bottom-area basis, and provide settlement substratum, 
food and shelter for a large number of microbial, algal and animal species 
(e.g. Kautsky et al., 1992; Middelboe et al., 2006; Rohde et al, 2008). 
Consequently, population, as well as community dynamics of Fucus spp. 
will affect the well-being of numerous community members (Altieri 
et al, 2007). 

In the atidal Baltic Sea, F. vesicidosus has declined substantially in large 
areas over the past 50 years (e.g. Berger et al, 2004; e.g. Vogt and 
Schramm, 1991, and references therein). The frequently observed upward 
shift of the lower depth limit has been attributed to direct and indirect 
effects of eutrophication such as decreased light availability and increased 
competition and sedimentation (Kautsky et al, 1986; Eriksson et al, 
1998; Berger et al, 2004; Korpinen et al, 2007b). However, changes in 
top-down regulation through a cascading effect of fish abundance on 
mesograzer abundance may contribute to local shifts by impacting com¬ 
peting filamentous algae (Eriksson et al, 1998) or the recruitment success 
of Fucus spp. (Korpinen et al, 2007a). Ongoing eutrophication in the 
Baltic changes the composition of F. vesicidosus associated invertebrate 
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community (Korpinen et al, 2010) and the decline ofE vesiculosus popu¬ 
lations is expected to entail even more substantial changes (Wikstrom and 
Kautsky, 2007). Salinity sets the distributional edge for F. vesiculosus in the 
northern Baltic Sea, which is expected to shift southwards with the 
climate change associated increase in precipitation in the northern catch¬ 
ment area (The BACC Author Team, 2008). 

A more complex situation exists along the rocky intertidal shores of 
the North Atlantic, as at least 2—3 species of Fucus commonly co-exist 
along with the large fucoid Ascophyllum nodosum and small, high intertidal 
Pelvetiacanaliculata. Latitudinal range shifts in seaweeds associated with sea 
surface temperature increases have been documented along the 
Portuguese coast (Lima et al, 2007). There is also evidence that the 
southernmost populations of F. serratus are maladapted and more vulner¬ 
able to stresses than more northerly core populations (Pearson et al, 

2009) . Although declines in overall abundance of fucoids (mostly A. 
nodosum) have been reported from a number of European areas (Davies 
et al., 2007), the specific changes in coverage of Fucus spp. have not been 
quantified. Reasons for fucoid decline have ranged from increased limpet 
grazing to increased numbers/intensity of storms. Similarly, the relative 
abundances of F. vesiculosus, F. distichus and F. spiralis are changing in the 
Canadian Maritimes, where F. vesiculosus is gradually increasing its cover 
within A. nodosum beds independent of harvesting regimes (Ugarte et al, 

2010) . The effect of climate change on range shifts in Fucus spp. specif¬ 
ically and along the rocky intertidal more generally is unknown 
(Thompson et al, 2002). Nevertheless, growing evidence suggests 
that the distributional ranges of Fucus spp. are in a process of reorganiza¬ 
tion as a consequence of species introductions and stresses associated 
with climate change, eutrophication and, more indirectly, overfishing 
(e.g. Ugarte et al, 2010). 

3.4. Genomic introductions and invasions 

Species introductions are also a major factor in shaping geographical 
distributions and community dynamics, including competitive displace¬ 
ment. For example, F. serratus was introduced to both Nova Scotia and 
Iceland from Europe within the last 100—200 years and in each area has 
expanded at a rate of 0.2—0.6 km yr~ (Coyer et al, 2006c; Brawley et al, 

2009) . In another example, Fucus evanescens (-F. distichus) was introduced 
to the Oslo fjord in the mid-1890s then expanded south into western 
Sweden (by 1933) and the southwestern Baltic (by 1992), where it expe¬ 
rienced less herbivore pressure both compared to native Fucus species and 
to populations in the native range (Wikstrom et al, 2006; Forslund et al, 

2010) . While the genus Fucus is considered a powerful invader itself 
(Williams and Smith, 2007), it may be impacted by invaders into its home 
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range. Thus, the introduced red alga Gracilaria vermiculophylla impacts 
the native Fucus vesiculosus doubly, by direct competition and by favour¬ 
ing Fucus grazers ( Idotea baltica) to which Gracilaria itself is relatively 
insensitive (Weinberger et al. , 2008). 

Another important aspect of species introductions is the notion of 
genomic invasions (e.g. Mallet, 2005), which can be inter-specific (i.e. 
congeneric contact between sibling species of Fucus) and/or intra-specific 
(i.e. conspecific contact among different biogeographic populations of a 
single species). In a natural example of genomic invasion, a process of 
‘genetic surfing' involving asymmetric introgression and exclusive north¬ 
ward spread of alien ( F. vesiculosus) organellar genomes following postgla¬ 
cial range expansion in F. ceranoides was recently discovered (Neiva et al, 
2010). In this case, the background nuclear genome of F. ceranoides was 
observed to retain remarkable integrity despite organellar introgression. 
Congeneric introductions of Fucus species create hybrid zones (as 
described in Section 3.1), thus complicating species identification and 
potentially resulting in long-term competitive displacement of parentals 
through evolution of new lineages via multiple generations of reproduc¬ 
tive FI hybrids, as well as by exchange of genes between parental species 
(introgression) via backcrossing with hybrids. In extreme cases, extensive 
backcrossing and interbreeding among hybrids could lead to complete 
homogenization of parental genomes or a hybrid swarm (extinction 
through hybridization) (see Coyer et al, 2007). 

Perhaps even more worrisome are conspecific introductions that go 
unnoticed. A new species of Fucus might be easily recognized as not 
part of the local marine flora but a foreign population of a resident 
Fucus species would not be. This is almost certainly the case in 
Nova Scotia (Brawley et al., 2009) where F. serratus was conclusively 
identified as an introduction from Europe and was of interest because it 
was not part of the Canadian Maritime flora. However, it is highly 
probable that European F. vesiculosus is also a part of the Canadian F. vesi¬ 
culosus flora although this has not been confirmed. The significance of 
this from the stress perspective lies in the fact that conspecific hybridiza¬ 
tion potentially leads to outbreeding depression and lower fitness of an 
individual or population, as has been shown in intertidal copepods 
(Edmands et al, 2005). Thus, increased extrinsic stresses might have even 
greater effects on such populations. Alternatively, it is also possible that 
intra-populational hybridization might simultaneously release beneficial 
variation and faster development (Edmands, 2008). Doubtlessly, human- 
mediated transport will increasingly bring hitherto separated species 
and populations into contact, as will climate change. For example, 
as the Arctic Ocean opens on a more permanent basis, contact between 
Pacific and Atlantic biotas will accelerate (Vermeij and Roopnarine, 
2008). 




64 


Martin Wahl et at. 



The Stressful Environment 


While marginal habitats may be dominated by single stresses (e.g. 
low salinity in the eastern Baltic, low light with increasing depth), envi¬ 
ronmental change typically involves the simultaneous shift of many stres¬ 
ses (Harley et ah, 2006; IPCC Climate Change, 2007; Darling and Cote, 
2008; Przeslawski et al, 2008). It is therefore important to understand the 
nature and influence of single stresses as well as their many combinations. 


4.1. Abiotic stress 

4.1.1. Single stresses 

Single stresses along environmental gradients can influence the perfor¬ 
mance and distribution of Fucus spp. although acclimation and adaptation 
to marginal and stressful habitats is possible to some degree (Kawecki, 
2008). At the largest scale, latitudinal gradients (thousands of kilometres) 
are governed primarily by sea surface temperatures and light regimes 
(Liming, 1990). At regional scales (tens of kilometres) there are transitions 
between fully marine and low salinity conditions in fjords and estuaries, as 
well as eutrophication gradients. Finally, the local scale (several metres) 
includes (1) the intertidal, where desiccation, high light and excessively 
high and low temperatures (air and water) create sharp and fluctuating 
gradients over a distance of a few metres and (2) the subtidal, where light 
and turbulence can be rapidly attenuated by depth, the former being 
further modulated by plankton blooms, turbidity or epibiosis. Below 
we provide an overview of the major abiotic stresses and how their effects 
may vary with Fucus life-history stage. 

4 .1.1.1. Irradiance Light is likely the most variable abiotic component 
of inter- and subtidal shores (Schubert et al, 2001) and is an obvious 
essential resource for all photo-autotrophic organisms. Nevertheless, 
both, an excess and deficit of light are significant stresses. Low-light 
stress occurs when irradiance is below the light compensation point of 
a photosynthetic species and compromises carbon accumulation (Lehvo 
et al, 2001; Middelboe et al, 2006). In the northern Baltic Sea, for 
instance, the lower distribution limits of F. vesiculosus correlate strongly 
with light penetration through the water; and in exposed areas (low 
sedimentation) F. vesiculosus individuals can survive down to 9—10 m 
depth (Kautsky et al, 1986; Eriksson and Bergstrom, 2005). When 
F. vesiculosus of the western Baltic, where the lower distributional limit 
is about 3 m, is transplanted to greater depths, some acclimation to the 
resulting low-light stress occurred (Rohde et al, 2008). However, the 
observed increase in chlorophyll-a concentration (by a factor 1.4 
following transplantation from 1 to 6 m depth) does not compensate for 
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the simultaneous decrease in irradiance (by a factor 7.3). With respect to 
seasonal irradiance regimes, F. vesiculosus is able to acclimate by shifting its 
compensation point between 35 p,mol m 2 s 1 (summer) and 
8 fjimol m “s 1 (winter) (Middelboe et al, 2006). Similarly, in the 
northeastern Pacific intertidal F. distichus becomes saturated at 
1000 ixmol m~ 2 s _1 in summer and 100 p,mol m~ 2 s in winter (Dethier 
and Williams, 2009). In the western Baltic the daily dose of irradiation 
(24 h average) falls below the compensation point of F. vesiculosus between 
October and March at 1 m depth and between August and April at 3 m 
depth (Wahl, unpublished) but—following the law of reciprocity—a few 
hours of high light per day may suffice for an equilibrated energy budget. 
Fucus species, like many brown algae, are able to store surplus solar 
energy in the form of mannitol (Lehvo et al., 2001), which sustains 
individuals through extended periods of sub-compensation. However, 
the observation that F. vesiculosus can survive the dark winter period in the 
western Baltic at 1 m depth but rarely occurs at 3 m depth or deeper 
(Rohde et al., 2008) may indicate limitations in the energy storage. 
Though it is likely that low-light stress affects all developmental stages 
equally, smaller individuals are often subject to additional shading by 
larger neighbours as a result of their understory position. In the western 
Baltic Sea, the development ofjuvenile F. vesiculosus over a growing season 
was limited by shade from the adult canopy in a field experiment 
manipulating both light and canopy cover (Eriksson et al. , 2006). 

Light exceeding the saturation level for photosynthesis may produce 
reversible photoinhibition as a protection mechanism or, at more extreme 
intensities, irreversible photodamage (Huppertz et al., 1990). 
Photodamage only occurs when the cellular mechanisms of photo-pro¬ 
tection are exceeded. Particularly severe stress may be caused by high- 
energy short wavelengths (UVA and B), which have negligible photo- 
morphogenic effects, but potentially large photo-destructive effects 
(Hanelt et al., 1997). Fucus species are often exposed to high levels of 
UV radiation, especially for intertidal populations during emersion. 
Repair mechanisms or shield compounds can accumulate in tissue to 
reduce UV stress, but may bear some metabolic cost and, thus, reduce 
physiological performance. While intertidal populations of Fucus exhibit 
regular photoinhibition, mostly in response to afternoon intensities, 
photodamage and reduction in physiological performance seem rare 
(Huppertz et al., 1990; Hanelt et al., 1997; Michler et al., 2002). 

The high mitotic activities during germination may be particularly 
susceptible to UV radiation (Wiencke et al., 2000), although some adap¬ 
tation to increased UV in particular habitats (e.g. UV absorbing 
compounds) can be expected. For example, embryos and juveniles of 
F. distichus in the upper intertidal on Spitsbergen, and F. serratus in the 
mid-intertidal on Helgoland, exhibited high sensitivity to UVA and 
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particularly UVB radiation at natural levels, while upper-intertidal F. 
spiralis from Helgoland were insensitive (Schoenwaelder et al, 2003). 
UV sensitivity in F. serratus was highest at the zygote stage and decreased 
non-linearly with age (Altamirano et al. , 2003). Phenolic phloroglucinol 
(a naturally occurring monomer of brown algal phlorotannins) tends to be 
more concentrated in younger tissue and was found to provide UV 
protection for developing embryos (Schoenwaelder et al., 2003). Early 
life stages of F. spiralis are less sensitive to ultra-violet radiation (UVR) than 
those of F. serratus presumably due to the higher density of phlorotannin- 
containing physodes in the former (Schoenwaelder et al, 2003). 

4 .i.i. 2 . Temperature Temperature stress is typically caused by the 
long-wavelength range of solar radiation (720 nm—300 pan), which is 
photosynthetically inactive, but thermally effective. The severity of 
stress caused by warming depends not only on the absolute temperature 
but also on the duration of exposure to it in air (emersion) and water 
(submersion), and the developmental stage considered. In intertidal areas, 
both high- and low-temperature stresses are generally more extreme 
during emersion. While germination success of F. vesiculosus is reduced 
by 90% at 25 °C compared to 15 °C in Kiel Bight (Germany), adults 
during emersion can survive transitory heating of up to 45 °C (Maczassek 
and Wahl, unpublished). Moderately warm temperatures, as in many 
other poikilotherm species, actually increased growth in A. nodosum 
(Keser et al, 2005). Severe high-temperature stress begins when the rate 
of protein denaturation cannot be neutralized by mechanisms such as 
increased production of chaperons (heat-shock proteins) or acceleration 
of protein biosynthesis (Csikasz-Nagy and Soyer, 2008). However, stress 
can also occur at sub-denaturating temperatures because impaired 
biochemical pathways may accumulate harmful intermediates (Davison 
and Pearson, 1996) and/or because repair mechanisms may be costly 
(Weidner and Ziemens, 1975). Elevated temperatures also negatively 
influence membrane properties that are crucial for gas exchange and 
nutrient uptake (Maheswari et al, 1999). 

Low-temperature stress reduces metabolic rates, but can be partially 
compensated by an increase in the concentration of key enzymes 
(reviewed in Middelboe et al, 2006). Nonetheless, growth may be 
severely reduced in winter (but see Lehvo et al, 2001; Dethier and 
Williams, 2009) and biomass losses to grazing, ice scouring or wave action 
cannot be regained before spring growth. As fucoid tissue can recover 
from extended freezing (Davison et al, 1989; Pearson and Davison, 1993; 
Pearson et al, 2000), the primary detriment of ice scouring is physical 
detachment of the holdfast. In some species, holdfasts remain intact 
following ice scour and may sprout new fronds in the subsequent growth 
season (see Section). In the northern Baltic Sea, subtidal F. vesiculosus 
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have two distinct strategies to avoid ice damage (Kiirikki and Ruuskanen, 
1996). In exposed areas, pack ice can reach several meters depth and an 
effective regeneration from holdfasts enables fast recovery of the vegeta¬ 
tion in the shallow zone. In sheltered areas, ice forms on the surface and 
here the large growing individuals decrease their buoyancy during winter, 
which causes them to lie flat on the bottom and thereby prevent them 
from freezing into the ice cover. 

4.1.1.3. Desiccation Extended exposure to air, especially in summer, 
can lead to variable degrees of desiccation for intertidal species of Fucus. 
Severe desiccation and the associated osmotic stress, especially when 
combined with high temperatures and light, can increase mortality 
(Pearson et al, 2009). For example, increasing duration of daily air- 
exposure for Fucusgardneri (F. distichus complex) had a negative effect on 
the growth rate of the adult thalli, as well as a cumulative negative effect on 
survival of all life-history stages (Wright et al, 2004). 

Juveniles are most susceptible to emersion and associated stresses 
(Schoenwaelder et al. , 2003; Henry and Van Alstyne, 2004). For example, 
fucoid germlings lose their photosynthetic ability within a few hours of air 
exposure, although they can recover quickly after re-immersion (Lamote 
et al, 2007). High mortality may result from desiccation, high light and 
high temperatures (Dudgeon and Petraitis, 2001; Wright et al., 2004). 
The smaller life stages, however, may gain a refuge from emersion stress 
under a canopy of adults (Lamote et al., 2007). For example, germling 
survival under an intertidal canopy approaches 100%, whereas survival is 
near zero on adjacent bare rock or on exposed habitats (Brawley and 
Johnson, 1991). 

4 . 1 . 1 . 4 . Pollution Pollution is defined as the accumulation of toxic 
compounds and/or abnormally high levels of nutrients in a local 
environment and/or tissue. Toxic compounds are an important marine 
stress because of bioaccumulation (Volterra and Conti, 2000). Fucus spp. 
readily adsorb and accumulate heavy metals (e.g. Mata et al., 2008) with 
concomitant negative effects on germination and growth (Andersson and 
Kautsky, 1996; Brooks et al, 2008; Nygard andDring, 2008). In addition, 
a biotic toxin (nodularin) produced by cyanobacteria is biosorbed by F. 
vesiculosus and induces oxidative stress (Pflugmacher et al, 2007). 
Increased nutrient loads mostly affect macroalgae indirectly through 
increased shading (pre-empting light) (e.g. Vogt and Schramm, 1991) 
and organic sedimentation (pre-empting space for recruitment; 
Eriksson and Johansson, 2003). Nutrient enrichment may enhance 
surface-fouling ranging from biofilms (Wahl et al, 2010) to ephemeral 
macroalgae (Korpinen et al , 2007b) and epifauna with diverse effects on 
the host (see Section 4.2). Nutrient enrichment also alters the selective 
environment favouring fast growing species (Berger et al., 2004) or by 
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changing the palatabilty of the thallus to grazers (Hemmi and 
Jormalainen, 2002). Direct effects arise through high concentrations of 
nitrate, which delayed attachment and decreased germination success in 
F. vesiculosus zygotes (Bergstrom et al., 2003). 

4.1.1.5. Osmotic stress Salinity stress mostly affects populations in the 
middle and uppermost intertidal. Rainfall increases the water potential 
(*/f) of the environment, whereas evaporation and desiccation obviously 
decreases the potential. Sensitivity to low-salinity stress varies enormously 
among Fucus species, populations and (presumably) life-history stages 
(e.g. Pearson et al., 2000). For example, F. vesiculosus and F. radicans have 
adapted to the low salinity of the Baltic Sea (Serrao et al., 1996a), but 
fertilization success in F. serratus decreases substantially with strongly 
reduced salinity: a fertilization success of 87% at 9 psu declines to 5% at 
6 psu (Malm et al, 2001). In a western Baltic population of F. vesiculosus, 
however, fertilization success was not affected by salinity between 7 and 
17 psu (Maczassek and Wahl, unpublished). Low salinity decreases 
swimming performance and fertilization ability of fucoid sperm (Serrao 
et al, 1996a) and increases the rate of polyspermy (Serrao et al., 1999b). 
The present range of F. vesiculosus in the Baltic Sea corresponds to the 
osmotic tolerance of its gametes (Serrao et al., 1996a), illustrating that 
environmental factors affecting performance of a particular life-history 
stage may be important determinants of distributional ranges and 
colonizing abilities. 

4.1.1.6. Water motion Water motion is an important consideration for 
shallow water organisms; as with light, too much or too little is stressful to 
an individual. In the case of Fucus and other large macrophyte 
genera, increasing flow speeds and/or wave action increases mechanical 
stress through increased drag, which must be counterbalanced by 
structurally resistant thalli and holdfasts (Haring et al., 2002). A negative 
correlation exists between increasing water motion and thallus size, as 
evidenced by extensive piles of beach wrack following high wave action. 
In F. vesiculosus and F. spiralis, a flow speed of 7—8ms _1 completely 
dislodges individuals larger than ~10 cm (Jonsson et al., 2006). 
Consequently, exposed sites have smaller-sized individuals than 
sheltered sites, either due to high mortality (especially of the larger 
individuals) or size reduction by thallus tattering of the larger 
individuals (Blanchette, 1997). The risk of dislodgement is more 
pronounced in areas where Fucus attaches to pebbles or mussel shells in 
lieu of bedrock. 

One strategy to deal with wave-induced forces is thallus toughening. 
In both the North Sea and the Baltic Sea, thalli from exposed F. vesiculosus 
were 30% more resistant to tear and breakage as compared with conspe- 
cifics from more sheltered sites (Nietsch, 2009). The formation of 
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bladders on F. vesiculosus seems to be a disadvantage under conditions of 
high turbulence and is negatively correlated with wave exposure 
(Burrows and Lodge, 1951). 

Fucoids enhance fertilization success by releasing gametes only under 
hydrodynamic conditions that optimize sperm—egg encounters (Pearson 
and Brawley, 1996; Serrao et al, 1996b; Pearson and Serrao, 2006). Water 
motion can interact with timing of gamete release during tidal cycles in 
species-specific ways, affecting gamete and zygote dispersal (Ladah et al, 
2008). Low wave exposure is also critical to successful attachment and 
survival of zygotes (Vadas et al ., 1990), and the rate offucoid recruitment is 
negatively correlated with wind speed (Lamote and Johnson, 2008). 
Consequently, water motion strongly influences local distribution pat¬ 
terns of Fucus spp. (Ladah et al., 2003, 2008). 

Conversely, calm conditions leading to a thicker boundary layer result 
in a slower absorption of nutrients and CO 2 , and a faster accumulation of 
0 2 and exudates at the thallus/water interface (J0rgensen and Revsbech, 
1985). The effects of depletion/enrichment, however, maybe compen¬ 
sated for by morphological adaptations, enhancing surface area and thallus 
rugosity that in turn reduces the boundary layer by enhanced microtur¬ 
bulence (Steen, 2003; Stewart and Embrey, 2003). Furthermore, calm 
conditions may facilitate sedimentation on algal thalli (Umar et al., 1998), 
resulting in shading and oxygen deficiency (if the sediment is rich in 
organic particles) (Kautsky et al, 1986; Duggins et al, 1990). High sed¬ 
imentation rates on rocky substrata also impede attachment of Fucus 
embryos and decrease the survival and growth of juveniles through scour 
and burial (Eriksson et al, 1998; Chapman and Fletcher, 2002; Schiel 
et al, 2006). In field and laboratory bio-assays, increasing amounts of 
deposited matter have been shown to decrease germling survival of F. 
vesiculosus by >60% (Berger et al, 2003), while experimental removal of 
sediment in the field increased recruitment twofold (Eriksson and 
Johansson, 2003). Thus, sedimentation regimes partly define recruitment 
opportunities of fucoid propagules, and strongly influence distributional 
patterns at the regional and local scale (Eriksson et al, 1998; Berger et al, 
2003; Eriksson and Johansson, 2005). 

4.1.2. Simultaneous abiotic stresses 

Stresses vary strongly in space and time. Spatial variation is usually linked 
to habitat properties such as depth, exposure or geographical position. 
Temporal variations may be rhythmic (diurnal, seasonal, decadal) or 
stochastic (rainfall, storms). Likewise stress sensitivity of an individual will 
vary with genotype, life-history stage or physiological phase (Fig. 2.2), 
which will determine the impact of a given stress at various scales (see 
Section 3.2). Temporal variations tend to decrease with increasing depth 
(e.g. Wahl et al, 2010) and lower latitudes. When different stresses 
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Figure 2.3 Hypothetical effects of single environmental factors on performance. 
Condition sets A and B exemplify conditions in different habitats or seasons where 
the combined effects of stressors determine a species’ performance. Environmental 
shift typically realizes as changes in the levels of several stressors, that is, a change in 
stress regime. If single stressor effects alone are considered, the stressor level with a 
minimum performance within a set of conditions determines performance. With 
additive stressor effects, performance is determined by the sum of negative stressor 
effects. With interactive stressor effects, performance cannot be predicted without 
knowledge on the quality of the interaction. 


fluctuate at different frequencies and/or amplitudes, the variability of their 
combined action may be complex (Fig. 2.3). 

Responses to combined abiotic stresses included antagonistic (5), 
synergistic (4) and additive (5) effects (Table 2.1). These studies collec¬ 
tively highlight the role of local conditions and local adaptation to stresses, 
not only in determining tolerance of single stresses but also in modulating 
the combined stress effects. For example, low nutrients and low salinity 
together decreased the performance of F. vesiculosus in an additive manner 
(Nygard and Dring, 2008). Antagonistic interactions were found in stres¬ 
ses including temperature, salinity, desiccation and UV, but the type of 
interaction was sometimes equivocal. In one case, low-salinity stress 
decreased the sensitivity of F. vesiculosus to UV in a population originating 
from the Atlantic (Nygard and Ekelund, 2006), but in a population from 
the Baltic Sea, the effects of salinity and UVR were additive. In another 
example, the combined effect of UV and temperature was either antag¬ 
onistic or additive depending on the response variable (Hoffman et al., 
2003). The microfouling—enhancing effect of warming in one study was 
buffered by light reduction (Wahl et al., 2010). Similarly, emersion and 
consequent partial desiccation antagonistically reduced the damage 
inflicted by high irradiation (Huppertz et al., 1990). Antagonistic inter¬ 
actions were also described for the combination of hyposalinity with 
unfavourably high or low temperature for various Fucus spp. (Munda 
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and Kremer, 1977), as well as of hypersalinity with stressful heat-shock for 
F. spiralis and F. vesiculosus (Li and Brawley, 2004). 

As an example of a synergistic stress effect, Maczassek (unpublished) 
found that experimental reductions of temperature (15—5 °C) and salinity 
(17—10 psu) lowered germination rates by 15 and 46%, respectively, when 
acting in isolation; but when the two stresses were present simultaneously, 
the germination rate was lowered by 90% (the additive expectation being 
61%). Other synergistic interactions were found between desiccation and 
wave action, between additions of different nutrients and between cold- 
temperature stress at parental stage and heat-shock at embryonic stage 
(Table 2.1). The last case provides an example of temporally separated 
synergistic stresses, where a stress in one ontogenetic stage (parent) ampli¬ 
fied the stress response in another stage (embryo) (Li and Brawley, 2004). 

4.1.3. Meta-analysis of abiotic stresses 

Results of the meta-analysis indicated that a single abiotic stress decreased 
algal performance of Fucus spp. by 40—50%, on average, as compared with 
the control (Fig. 2.4a and b). However, depending on the type of exper¬ 
iment, the average stress effect was significant only in experiments con¬ 
ducted in recruit/germling stages (lab experiments only, Fig. 2.4a) or only 
with adults (mainly field experiments, Fig. 2.4b). 

Averaged over all combinations, the effect of co-occurring abiotic 
stresses was highly variable but on average additive (Fig. 2.4a, Fig. Sla). 
The same was true with the combined effect of abiotic stress and nutrient 
enrichment (Fig. 2.4b, Fig. Sib). It is worth noting that nutrient enrich¬ 
ment alone had neither negative nor positive average effects (Fig. 2.4b). 
However, behind the average combined effects there was considerable 
variability of interactive effects within individual experiments (Fig. Sla 
and b) as noted earlier. This variability illustrates how a second abiotic 
stress may modulate the impact of a first abiotic stress and how much 
the outcome depends on the nature of stresses combined and the 
environmental setting. 


4.2. Biotic stress 
4.2.1. Single biotic stresses 

The biotic stresses of grazing, fouling and competition vary enormously 
over space and time. For example, the number of grazer and epibiont 
species decreases from the fully marine North Atlantic Ocean into the 
brackish Baltic Sea (e.g. Kautsky et al., 1986; Ojaveer et al, 2010) as does 
genetic diversity (Johannesson and Andre, 2006). Although the various 
biotic stresses interact in myriad ways in nature, understanding the effect 
of single biotic stresses is a necessary first step to understanding the higher 
order interactions. 
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Figure 2.4 Effects of single and combined stresses based on a meta-analysis of 41 
experiments extracted from articles in Table 2.1. Effect size (LRR) is plotted against a 
single stress or stress category (a—f) and their combined effects (see Table 2.2) 
following the method of Hedges et al. (1999). The number of experiments 
compared is given above each bar. The error bar depicts the 95% confidence 
interval of the mean effect ratio. When the confidence interval does not include 
zero, the stress is considered significant. Note that for different meta-analyses (a—f) 
different sets of experiments were used, which accounts for differences in the 
magnitude of a particular result. See Section 2 for details of the assumptions and 
calculations using the method of Hedges et al. (1999). See Figure SI for results based 
on the method of Gruner et al. (2008). 
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Figure 2.4 Continued 


4 . 2 . 1 . 1 . Competition and epibiotism The strong intertidal zonation 
patterns typically exhibited by Fuchs species (Fig. 2.1), together with 
other macroalgae and sessile animals, clearly reveals the importance of 
competition and competitive displacement among these species 
(reviewed in Karez and Chapman, 1998; e.g. Schonbeck and Norton, 
1978; Schonbeck and Norton, 1980). Ultimately, the ability of Fucus to 
persist and successfully compete intra- or inter-specifically is determined 
by its growth rates (Airoldi, 1998; Johnson et al., 1998), canopy height 
(Dayton et al, 1984), reproductive output (Malm and Kautsky, 2003), 
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production of allelopathic compounds (Gross, 2003) and tolerance 
towards biotic and abiotic stresses (e.g. Worm and Chapman, 1998; 
Pearson et al, 2009; Schiel, 2009). 

Competition for light and space occurs among Fucits individuals (con- 
specifics and congeners), as well as with epibionts attached to their thalli 
surfaces (Creed et al., 1996; Karez and Chapman, 1998; Karez, 2003a,b). 
In general, competition varies with ontogenetic stage. For example, 
periphyton and ephemeral algae may decrease germling survival 
(Lubchenco, 1983; Lotze et al., 2000), whereas microalgae and sessile 
invertebrates decrease growth, interfere with nutrient absorption and 
increase drag on the thallus (Brock et al., 2007; Jormalainen et al., 
2008b; Rohde et al., 2008; reviewed in Wahl, 1996; see also Wahl and 
Mark, 1999; Wahl et al., 2010). 

Epibionts on Fuats thalli also include bacterial biofilms (microfoulers), 
a largely unknown interaction that may be important, particularly as a 
host-specific assemblage of microbes influences many interactions of the 
host with its environment (Dobretsov, 2008; Lachnit et al., 2009). 
Furthermore, microfouling pressure is a conditioning phase in the fouling 
process and varies substantially with temperature, light and depth. For 
example, in the western Baltic fouling pressure (substratum occupation 
per unit time) by macroscopic species is 30 times lower in February as 
compared to July, and the number of simultaneously settling species is 10 
times lower in winter than in summer (Thomsen and Melzner, 2010; 
Wahl et al, 2010). By colonizing the functional interface between the host 
alga and the environment, dense epibiosis will clearly modulate most 
interactions of the host (Wahl, 2008a). 

4.2.1.2. Herbivory and parasitism Herbivory on Fttats, primarily by 
isopods and, to a lesser extent, by amphipods and snails, ranges from mild 
to severe and may differentially impact Fttats from the germling stage 
through the adult stage (Pennings et al., 2000; Long et al, 2007; Dethier 
and Williams, 2009). High grazing rates by isopods (mainly the genus 
Idotea) can eliminate local Fttats stands (i.e. Engkvist et al., 2000), while 
grazing by herbivorous fishes is less common in the distributional range of 
Fuats (Connell and Slatyer, 1977; Floeter et al, 2005). 

In the species-poor northern Baltic Sea, grazing pressure by isopods 
varies seasonally (e.g. Kotta et al, 2006): low in winter (when isopods are 
inactive) and early summer (when isopod densities are low); high in late 
summer, which is characterized by high densities of juvenile isopods 
(Engkvist et al, 2000; Korpinen et al, 2010). Grazing pressure also varies 
spatio-temporally with exposure (wave action) (e.g. Jonsson et al, 2006) 
and the patchy distribution of consumers (e.g. Korpinen and Jormalainen, 
2008b). Isopod grazing pressure may also vary with fishing pressure on 
isopod predators (e.g. Bostrom and Mattila, 2005; Jormalainen and 
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Ramsay, 2009). Simultaneous grazing by several consumer species may 
interact. Thus, defences developed by F. vesiculosus against one herbivo¬ 
rous species may reduce grazing by another species (Long et al ., 2007; Yun 
et al, 2010). 

Very little is known about parasitism or infections in Fucus. Coles 
(1958) identified parasitic nematodes that caused galls on F. vesiculosus 
and F. serratus along the Cornwall coast and Zuccaro et al. (2008) recently 
detected a number of fungal species associated with F. serratus. Unlike A. 
nodosum, which harbours an obligate endophytic fungus (Xu et al, 
2008, and references therein), few, if any mycobionts are associated with 
species of Fucus (but see Parker and Chapman, 1994 for impact by endo¬ 
phytes) . However, the potential for increased biotic interactions involving 
invasive parasites or pathogens is on the rise in many marine systems 
(e.g. Torchin et al., 2002). 

4.2.2. Simultaneous biotic stresses 

Because Fucus exists within a multi-species community, it will simulta¬ 
neously or sequentially interact with a multitude of micro- and macro- 
foulers, competitors (annual or perennial macroalgae, sessile animals) and 
consumers (crustaceans, molluscs), as well as an unknown number of 
pathogens and parasites. Consequently, the co-occurrence of different 
biotic stresses is expected to be as common as the co-occurrence of abiotic 
stresses. However, only few studies have investigated the interactive effects 
of two or more biotic stresses (Table 2.1). 

Several studies illustrate the complexity of interactive biotic effects. 
The competitive impact of turf algae (Rhodophyta) on F. vesiculosus is 
alleviated by limpet grazing, but not the similarly negative effect of larger 
canopy algae (Jenkins et al, 1999). Similarly, limpet grazing and compe¬ 
tition with turf algae had an additive effect on Fucus recruitment, but at 
the same time limpet grazing inhibited recruitment of A. nodosum that has 
the potential to suppress Fucus by canopy effects (Cervin et al, 2005). In 
another study, epibiotism and grazing had an additive effect on the growth 
of F. vesiculosus, but the effects were functionally interrelated in the sense 
that fouling attracted grazers (Jormalainen etal. 2008b). Finally, in F. disti- 
chus, grazing by littorinids and gammarids has been reported to reduce 
the rate of fatal endophyte infections (Parker and Chapman, 1994). Thus, 
if grazers feed on both the host and the epibiota, epibionts may attract 
grazers with consequent ‘shared doom’ for both, epibionts and host. On 
the other hand, if the epibionts are resistant to grazers, the host may gain 
associational resistance. 

4.2.3. Results of the meta-analysis of biotic stresses 

Experiments with two biotic stresses revealed additive (5), synergistic (1) 
and antagonistic (2) interactions (Table 2.1). The meta-analysis of single 
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biotic stresses (competition, herbivory and epibiotism) showed a mean 
10—50% decrease in performance relative to controls (Fig. 2.4c and d) and 
no differences in biotic stress responses among life-history stages. With all 
life-history stages combined, both the single and combined effects of 
biotic stresses were statistically significant, which was not the case when 
the data were stratified by life-history stages (Fig. 2.4c). The meta-analysis 
indicated that the combined effect of two biotic stresses, in general, was 
additive (Fig. 2.4c, Fig. Sic). Still, the high variation in interaction effect 
sizes reflects the contingency of multi-stress responses in which interactive 
effects are more complex and depend on the host spectrum and feeding 
preferences of the grazers (Wahl and Hay, 1995; Karez et al., 2000). 


4.3. Protection against abiotic and biotic stresses 

Fucus spp. have evolved several defence mechanisms for a variety of 
abiotic and biotic stresses. One example is development of an elastic, but 
tough thallus that resists mechanical damage (e.g. Nietsch, 2009), whereas 
wound healing and re-sprouting from holdfasts should favour fast popula¬ 
tion recovery when damage has occurred (e.g. Malm et al., 2001; 
Tatarenkov et al. , 2005). High concentrations of phlorotannins may offer 
protection from UV radiation (Pavia et al., 1997; Swanson and Druehl, 
2002) or deter herbivory (e.g. Hay and Steinberg, 1992; Jormalainen and 
Honkanen, 2008). Metabolic repair mechanisms of UV damage exist, but 
may be jeopardized under multiple stresses (e.g. Hoffman et al., 2003). 
Palatability of Fucus species to grazers is affected by a number of lipophilic 
and polar metabolites; in particular, it is augmented by high concentrations 
of mannitol and proteinic amino acids, as well as by high water content, and 
reduced by high concentrations of phlorotannins and high toughness of the 
thallus (Tuomi et al., 1989; Hemmi and Jormalainen, 2004; Nietsch, 2009; 
Jormalainen et al., 2011; Weinberger et al., 2011). Although the reduction 
of nutritive value (e.g. fewer storage compounds such as mannitol) may save 
energy and discourage grazing (Viejo and Arrontes, 1992; Worm et al, 
1999; Engkvist et al, 2004), it may also jeopardize the chances to survive 
seasons of low energy supply (Lehvo et al., 2001). 

Chemical defences are present in many terrestrial and marine species, 
but their production, storage and ultimate detoxification may well impose 
fitness costs on the producer (reviewed in Stamp, 2003; Pohnert, 2004). 
Costs of chemical defence traits have been documented in macroalgae 
(Dworjanyn et al., 2006; Jormalainen and Ramsay, 2009), but because the 
relevance of costs varies with resource availability (Davison and Pearson, 
1996; Dethier et al, 2005), and may vary with respect to the measured 
fitness component, costs have not been detected in all studies of chem¬ 
ically defended algal species (Rohde and Wahl, 2008; Pansch et al, 2009; 
Appelhans et al., 2010). In addition, net defence costs are only expected to 
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arise under conditions in which natural enemies do not impose significant 
performance costs and thus when defences are superfluous. 

Nevertheless, at least four Fucus species employ chemical defences 
against fouling and herbivory. For example, F. vesiculosus alters the com¬ 
position of epibiotic fouling by producing a combination of pro- and anti¬ 
fouling compounds (e.g. Lachnit et al, 2010), and macrofoulers such as 
barnacles are deterred by secondary metabolites of F. vesiculosus (Brock 
et al., 2007; Jormalainen et al., 2008b; Rohde and Wahl, 2008) and F. 
evanescens (-F. distichus, Wikstrom and Pavia, 2004). Seasonal fluctuations 
in anti-microfouling chemical activity have been reported in F. vesiculosus 
(Wahl et al, 2010), but whether this is resource-driven (availability of 
energy) or demand-driven (fouling pressure) and whether similar rhythms 
exist in anti-macrofouling defences has not yet been investigated. Induced 
chemical defences against herbivores have been detected in F. distichus (van 
Alstyne, 1989), F. vesiculosus (Rohde et al., 2004), F. serratus (Rohde and 
Wahl, 2008) and F. evanescens (Rohde and Wahl, 2008). While some 
chemical defences are produced or deployed at variable intervals after 
the onset of herbivory (‘induced’), they may also be reduced before 
herbivory is terminated (Weinberger et al., 2011). 


4.4. Modulation of stresses 

When two negatively interacting species exhibit unequal sensitivities 
towards a given abiotic stress, the least sensitive species may indirectly 
benefit, as the more sensitive species is more strongly inhibited. 
Furthermore, when abiotic stress jeopardizes the production of defence 
chemicals, the overall stress impact will be magnified. Thus, when the 
costs of defence and of stress compensation mechanisms must be paid from 
the same limited resource, defences may decline, subsequently increasing 
biotic stress between two species. Accordingly, acute abiotic stress should 
most impact inducible defences that are produced de novo, rather than as a 
by-product of some primary metabolic pathway. 

The meta-analysis revealed that, on average, interactions between 
abiotic and biotic stresses were additive (Fig. 2.4d and e), although there 
was a wide scatter of variable effects among separate experiments, varying 
from synergistic to highly antagonistic ones (Fig. Sid and e). Below, we 
focus on studies examining how abiotic stress modulates various biotic 
interactions among Fucus species, their grazers, competitors and epibiota. 
Although several studies did not utilize a factorial design, they do provide 
indirect evidence about the possible type of stress interactions. 

4.4.1. Irradiance effects on biotic interactions 

Sensitivity to high and low light stress contributes to the patterns of 
intertidal zonation exhibited by Fucus species (Karez and Chapman, 
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1998). Although epibiont shading may help to protect Fucus thalli against 
excessive irradiation (as shown for other organisms by Przeslawski et al., 
2008), the effect is often transient, as epibiota rapidly increase and thus 
exacerbate light limitation for the host (Oswald et al, 1984). Low light 
may weaken anti-microfouling defences in F. vesiculosus and uncontrolled 
bacterial fouling, in turn, can reduce available light by more than 90% 
within a few weeks (Wahl et al, 2010). Low light and fouling had additive 
negative effects on growth (Rohde et al, 2008) and both macro- and 
microepibionts may raise the light compensation depth by 2 m in the 
western Baltic (Rohde et al, 2008). Susceptibility to light stress (and 
competitive interactions) also is influenced by emersion time. For exam¬ 
ple, F. serratus is more susceptible to high light stress during emersion than 
F. vesiculosus (Malm and Kautsky, 2003) and this can alter local distribution 
patterns of the species. 

Low light stress also affects the susceptibility of Fucus spp. to herbivory, 
although there is no indication that irradiation affects the behaviour of 
herbivores directly (e.g. Eriksson et al, 2006). For example, light limita¬ 
tion can reduce the production of chemical defences in F. vesiculosus thus 
enhancing consumption by the isopod I. baltica (Weinberger et al, 2011). 
At the same time, reduced growth rates under light limitation (e.g. Rohde 
et al., 2008) increases the effect of herbivory, as grazed tissue is more slowly 
replaced; an effect that may be countered by decreased mannitol produc¬ 
tion in low light, which subsequently renders the alga less attractive to 
isopod grazers (Weinberger et al, 2011). A synergistic effect of grazing 
and low light stress has been demonstrated, as grazing losses in F. vesiculosus 
increase under low light, possibly by decreasing resistance to herbivory 
(fewer defence chemicals) (Jormalainen and Ramsay, 2009). 

Deleterious effects of UVR on the defence capacity have been shown 
in some macroalgae (Cronin and Hay, 1996; Pavia et al, 1997), but not in 
others (Macaya et al, 2005). While enhanced UVR exposure in various 
Fucus spp. increases phlorotannin concentration (Pavia and Toth, 2000; 
Schoenwaelder et al, 2003; Henry and Van Alstyne, 2004), it is unclear 
whether increased levels constitute an anti-herbivore defence or protec¬ 
tion from UVR (reviewed in Amsler and Fairhead, 2006). Defence 
against epibionts is another consideration and in F. vesiculosus, is strongest 
at moderate irradiance (Wahl et al, 2010) or levels typical of its normal 
depth distribution (1—2 m) in the western Baltic Sea. 

4.4.2. Temperature effects on biotic interactions 

Increasing temperatures (>15 °C) enhance biotic stress in Baltic Sea Fucus 
by (1) increasing micro- and macrofouling rates (Wahl et al, 2010); 
(2) increasing grazing rates of I. baltica in summer (but not in the other 
seasons) (M. Zimmer, personal communication) and (3) decreasing levels 
of defence chemicals (under both high and fluctuating temperatures) 
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(Weinberger et al. , 2011). In contrast, anti-fouling defence chemicals in F. 
vesiculosus were unaffected by a 2-week exposure to temperatures between 
8 and 23 °C (Wahl et al., 2010). Thus, temperature and grazing stress can 
generate synergistic interactions, whereas temperature and fouling stress 
are additive. 

Freezing temperatures can shift competitive interactions among Fucus 
spp. Laboratory experiments revealed that F. vesiculosus had a higher tol¬ 
erance to freezing than its congener F. serratus: aerial exposure to —15 °C 
for 1 h killed F. serratus, while F. vesiculosus survived (Malm and Kautsky, 
2003). Thus, for this species pair freezing and competitive stress is likely to 
be antagonistic for F. vesiculosus (direct physiological stress partially com¬ 
pensated for by reduced competition) and negatively synergistic for 
F. serratus. 

4.4.3. Salinity and emersion effects on biotic interactions 

Emersion may decrease salinity by removing seawater and/or exposing 
intertidal organisms to rain and/or freshwater runoff or increase salinity by 
evaporation. The combination of diverse factors associated with emersion 
will determine the strength or even the sign of stress during emersion. In 
F. vesiculosus, for example, emerged individuals experienced a transient 
refuge from aquatic herbivory and epibiont settlement, as well as an 
accumulation of anti-fouling defence chemicals at the thallus surfaces, 
which produced a more effective defence after re-immersion (Brock 
et al., 2007). Thus, emersion and fouling may act as antagonistic stresses. 

4.4.4. Water motion effects on biotic interactions 

Canopy sweeping or whiplash removes or tatters early post-settlement 
stages (Vadas et al., 1992) and illustrates how wave action, which is not a 
strong stress for small individuals or recruits, may be lethally amplified by 
the presence oflarger (con)specific individuals. Thus, whiplash represents 
a synergistic stress for juveniles and may shift competitive interactions 
between F. serratus and F. vesiculosus in favour of F. serratus because F. serratus 
juveniles are less sensitive to whiplash (Vadas et al., 1992; Malm and 
Kautsky, 2003). However, because F. vesiculosus can regenerate from hold¬ 
fasts, the slight competitive advantage of F. serratus recruits ultimately may 
be negligible (Kiirikki and Ruuskanen, 1996; Malm and Kautsky, 2003; 
Malm and Isaeus, 2005). 

Whiplash also reduces the abundance of competing filamentous epi¬ 
phytic and understory macroalgae (Kiirikki, 1996), as well as the effect of 
swimming grazers (either intertidally or subtidally) (Davis et al., 2003), 
thus creating an antagonistic interaction between wave action and both 
competition and herbivory. As isopods (I. baltica) experience more attach¬ 
ment difficulties on F. serratus than on F. vesiculosus due to differences in 
surface morphology, in high water motion grazing intensity is shifted to 
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F. vesiculosus (Engkvist et al, 2004). Thus, the effects of water motion and 
grazing in a mixed stand of F. serratus and F. vesiculosus are antagonistic for 
F. serratus and synergistic for F. vesiculosus. Though not tested, it is con¬ 
ceivable that where the two species compete, F. serratus may benefit from a 
combined wave action and grazing stress, as long as F. vesiculosus density 
remains at or above some threshold level. In the more exposed mid- to 
high-intertidal zone, non-swimming grazers such as limpets (Patella vul- 
gata) are dominant and their grazing preference for F. vesiculosus and F. 
spiralis over F. serratus further affects the competitive hierarchy by allowing 
F. serratus to move upshore (Jonsson et al, 2006). An additional consider¬ 
ation is that while the risk of thallus dislodgement is higher at wave- 
exposed sites, water motion may indirectly lead to reduced grazing by 
inducing increased toughness of the thallus (Nietsch, 2009). 

4.4.5. Eutrophication effects on biotic interactions 

Nutrient enrichment has been the most frequently studied abiotic mod¬ 
ifier of biotic interactions, with individual responses being additive (2), 
antagonistic (3) or synergistic (9) (Table 2.1). 

With respect to nutrient enrichment—herbivory interaction, grazers 
counteract the negative effect of high nutrient availability on colonization 
success of F. vesiculosus by reducing competition for space (Lotze et al., 
2000). In later stages, gastropod grazing counteracts nutrient stress by 
removing epibiota (Jormalainen et al, 2003; Raberg and Kautsky, 
2008). Field experiments on F. vesiculosus in the Baltic Sea (where abun¬ 
dant gastropod grazers regularly graze early recruits and nutrient enrich¬ 
ment severely restricts recruitment success) showed that negative effects of 
grazing and nutrient enrichment on recruitment success remained addi¬ 
tive, but varied depending on the community structure of other grazers 
(Korpinen and Jormalainen, 2008a) and the presence of canopy species 
(Eriksson et al ., 2007). On the other hand, synergistic effects of nutrient 
enrichment and herbivory may arise when nutrient enrichment 
improves the nutritional quality of Fucus thalli. Thus, high-quality Fucus 
tissue greatly benefits growth and reproduction of I. baltica (Hemmi 
and Jormalainen, 2002) and the amphipod, Gammarus locusta (Kraufvelin 
et al, 2006), which concomitantly leads to higher grazing pressure. 
Palatability of F. vesiculosus to littorinid grazers also increases as a result 
of increased nutritive value (increased N: Yates and Peckol, 1993; 
increased mannitol: F. Weinberger, personal communication). Finally, 
negative effects of grazing on F. vesiculosus may be stronger under increased 
nutrient conditions (Hemmi et al., 2005). 

With respect to nutrient enrichment and competition, interactions 
were mainly synergistic (Table 2.1), suggesting that nutrient enrichment 
benefits competitors/epibiota of Fucus more than Fucus itself 
(Jormalainen et al, 2003, 2008b). As expected, growing F. vesiculosus in 
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a nutrient enriched environment without epiphytes positively affected 
growth (Jormalainen et al, 2003). Thus, two effects of epibiosis have a 
synergistic impact on their Fucus host by shading and inhibiting 
nutrient uptake (through insulation or a ‘smothering’ effect). 
Furthermore, nutrient enrichment may completely reorganize compet¬ 
itive hierarchies as eutrophication typically favours ephemeral, fast 
growing algae (Berger et al, 2004) at the cost of slow growing perennials 
(Korpinen et al, 2010). 

The addition of an abiotic stress to two biotic stresses may render the 
biotic—biotic interactive effects more complex and non-additive. For 
example, under ambient levels of nutrients, grazers had a large positive 
and competitors a slightly negative effect on growth of F. vesiculosus 
recruits, with an additive or a minor synergistic interaction (Worm 
et al, 2001). Under nutrient enrichment, however, the combined effect 
turned highly antagonistic as the increased competition cancelled out the 
positive grazing effect. 

4.4.6. Biotic modulation of abiotic stress 

Biotic interactions may, in turn, modulate the impact of abiotic stresses. 
Thus, as shown for non-Fucus spp., epibionts may protect their host from 
UVR radiation and desiccation stress (Penhale, 1977; Przeslawski et al, 
2008). Flowever, they may also render the fouled host more vulnerable to 
drag (Witman and Suchanek, 1984; Wahl, 1997). Furthermore, canopy¬ 
forming algae, which might act as inter- or intra-specific competitors, 
may at the same time shelter small Fucus from abiotic stress during emer¬ 
sion (Brawley and Johnson, 1991; Lamote et al, 2007). 

Results of the meta- analysis showed that the average combined effect of 
herbivory or competition (or both) with nutrient enrichment did not differ 
from additivity (Fig. 2.4e and f). However, the nutrient enrichment— 
herbivory interaction was less consistent than the nutrient enrichment- 
competition interaction, which was also found in the individual studies. 


5 . Genetic Levels of Stress Response 

On a temporal scale, responses to stress range from fast acclimation 
by phenotypic plasticity to slower adaptation by stress-driven selection 
(i.e. adaptive evolution or ecotypic differentiation). The efficacy of both 
relates to an individual’s or a population’s genetic composition. 
Traditionally, physiologists and ecologists have been more interested in 
immediate responses, while population geneticists have been more inter¬ 
ested in changes over multiple generations (Allendorf and Luikart, 2007). 
However, it is now realized that the temporal differential between 
ecological and genetic changes can be small or negligible and most species 
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are not driven to local extinction before genetic factors play a role 
(Spielman et al, 2004). In short, dynamic selective regimes in the envi¬ 
ronment can have a significant and nearly simultaneous impact on both 
ecological and evolutionary processes, which subsequently shape genetic 
structure and performance of populations. The ability to distinguish these 
effects is an emerging challenge. 

5.1. Sensitivity versus genetic diversity of a population: 
stress from the evolutionary perspective 

The extrinsic component of stress (abiotic and biotic factors in the envi¬ 
ronment) shapes the phenotype through natural selection and leads 
to changes in allele frequencies and new gene complexes, collectively 
resulting in adaptation (i.e. the integration of the phenotype with 
the genotype). But it is crucial to also recognize the intrinsic component 
of stress tolerance (Bijlsma and Loeschcke, 2005), estimated as the 
standing genetic variation or diversity (evolutionary potential) of a 
population. 

Genetic diversity stems from the relationship between the effective 
population size (roughly the number of individuals that successfully 
contribute to the next generation) and the processes of gene flow, 
selection, genetic drift and mutation. Effective population size (N e ) is a 
crucial parameter for assessing the adaptive potential of a population. 
Large sexually reproducing populations with high N c have more genetic 
diversity and collectively experience fewer negative ramifications from 
the stochastic effects of genetic drift (random gain and loss of alleles 
through time). Moreover, at least some individuals within large popula¬ 
tions are likely to have genotypes that are already able to cope (pre¬ 
adapted) with a new selective pressure such as a higher temperature. In 
contrast, populations with low N c —often arising from habitat fragmen¬ 
tation or isolation—tend to have less genetic diversity, are highly subject 
to the negative effects of genetic drift, and are collectively more likely to 
be affected by stress due to their low adaptive potential (Harley et al., 
2006). As pre-adapted genotypes in these populations would be few or 
absent, it follows that they also are at higher risk of local extinction. 
Estimating iV e is extremely difficult. Direct estimates require extensive 
tagging studies over long time periods, whereas indirect estimates 
are based on temporal changes in allelic frequencies. Recent work on 
F. serratus using microsatellites suggested that despite large census popula¬ 
tion sizes lV e was very small with a concomitant loss of fitness (Coyer et al ., 
2008). If loss of genetic variation is combined with a deteriorating 
environment, then a negatively synergistic effect is likely to arise in 
which a population’s susceptibility and tolerance to extrinsic stress is 
magnified. 
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5.2. Indirect detection of adaptation to stress on relevant 
spatial scales 

Neutral loci are not under selection and as they are affected only by gene 
flow and genetic drift, cannot be used to directly measure adaptation. 
Neutral loci do, however, provide necessary basic information for studies 
of adaptation by resolving genetic population structure, measuring gene 
flow and providing indirect evidence for adaptation when correlated with 
environmental parameters. 

Genetic population structure can identify spatial scales over which 
selection (e.g. adaptation to stress at a population level) may act and has 
been assessed at one or more spatial scales for F. serratus (Coyer et al. , 2003; 
Hoarau et al, 2007), F. vesiculosus (Engel et al., 2005; Tatarenkov et al, 

2007) , F. distichus (Coyer et al, 2011c) and F. spiralis (Engel et al, 2005; 
Coyer et al, 2011a) using microsatellite (Coyer et al, 2002a,b,c, 2009; 
Engel et al, 2003, 2005) and mitochondrial loci (mt DNA) (Hoarau et al, 
2007; Engel et al, 2008). Population structure may be homogenized by 
gene flow as high gene flow dampens local adaptation to a stress (Kawecki, 

2008) . However, estimates of gene flow for F. serratus based on population 
differentiation at the regional scale (1—100 s km) suggest that gene flow 
can be as little as 500 m to 2 km (Coyer et al, 2003), between 10 m and 
<1 km in F. vesiculosus from the Baltic Sea (Tatarenkov et al, 2007) and 
150 m to 20 km in Brittany (estimated from Billard et al , 2005a). In the 
Kattegat Sea, gene flow distances for F. serratus increased to 70 km (Coyer 
et al, 2003). Taken together, the aforementioned results illustrate how 
oceanographic currents and local coastal topologies affect dispersal and 
variation (see below). 

Changes in putatively neutral allele frequencies along various stress 
gradients in replicated geographic locations provide strong correlative 
evidence for adaptation. This, however, is not conclusive as long as the 
actual genes involved in the stress response remain unknown. Conversely, 
strong changes in phenotypic traits may be observed across a gradient that 
are not mirrored by neutral allele frequencies. In Fucus spp., both of these 
situations have been encountered. For example, no correlation was found 
in the neutral allele frequency of F. vesiculosus along a gradient in salinity 
from the North Sea (35 psu), to the Kattegat, and the SW Baltic Sea 
(12 psu) using neutral markers (Tatarenkov et al, 2005; Johannesson and 
Andre, 2006), despite strong phenotypic differences in emersion stress 
tolerances (F. vesiculosus is submerged in the atidal Baltic Sea) (Pearson 
et al,2000) and fertilization success (Serrao et al, 1996a,b). Deeper within 
the Baltic Sea (from Bornholm and then extending north into the 
Bothnian Gulf), however, a strong shift in allele frequencies is correlated 
with decreasing salinity and a transition to clonal reproduction of F. radi- 
cans (Tatarenkov et al , 2005). The salinity shift is even considered a strong 




8 4 


Martin Wahl et at. 


enough selective force to contribute to the recent and rapid formation of 
the newly described endemic species F. radicans within the past 400—2000 
years (Pereyra et al., 2009). In F. serratus, there is no change in neutral allele 
frequencies over salinity gradients either within or between the Kattegat 
and Baltic Seas (Coyer et al., 2003) or along steep salinity gradients in a 
Norwegian fjord (Coyer et al., 2011b). However, using a genome scan 
approach these authors found evidence for selection in six candidate loci 
(three derived from expressed sequence tag (EST) libraries and three from 
anonymous libraries; see Section 5.3) of which at least one was related to 
osmotic stress (Coyer et al., 2011b). Thus, a lack of correlation using 
neutral markers does not imply a lack of selection. 

At the local scale of a single patch or stand, populations may display 
variation in quantitative genetic traits as has been shown in F. vesiculosus 
from the Baltic Sea regarding the resistance to both epiphytes and grazers 
(Jormalainen and Honkanen, 2008; Jormalainen et al., 2008b; Koivikko 
et al., 2008; Jormalainen and Ramsay, 2009), as well as to concentrations 
of phlorotannins (Jormalainen and Honkanen, 2004). Similarly, strength 
and/or inducibility of anti-herbivory and anti-fouling defences may vary 
enormously among conspecifics of F. vesiculosus, F. serratus and F. evanescens 
(Rohde et al., 2004; Wahl et al., 2010; E Symanowski, unpublished). 
This variation indicates the potential for local adaptation to natural 
enemy regimes. 

5.3. Direct detection of adaptation to abiotic and biotic stress 

The interplay between environment and genes is the focus of the new field 
ofEEG (Ouborg and Vriezen, 2007; Tautz et al. , 2010). Until full genome 
sequences are available for the various Fucus spp. (whose genomes range in 
size from 800 to 1000 Mb, Kapraun, 2005), a more rapid and economical 
approach is to develop EST libraries, which represent the expressed genes 
of the individual under a given stress condition (compared to a recovery 
condition). EST libraries provide a rich source for gene discovery and the 
development of gene-linked markers that can be tested for evidence of 
selection using genome scans (Storz, 2005). Comparative gene expression 
studies based on EST libraries provide a snapshot of organismal physiology 
and a fine-scale picture of plastic and constitutive changes in cellular 
metabolism associated with acclimation and genetic adaptation to the 
stress under investigation. Thus, EST-linked markers and genomic scans 
can distinguish phenotypic plasticity (acclimation capacity) from ecotypic 
differentiation directly rather than correlatively. 

Salinity dines and heat-shock in Fucus spp. have recently been inves¬ 
tigated in an EEG framework. Species distributions from the North, 
Kattegat and Baltic Seas are governed by both salinity and (emersion) 
temperature gradients. Replicated clinal variation along such parallel 
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environmental gradients provides good evidence for natural selection in 
the wild. For example, four populations of F. serratus along a 12-km long 
salinity gradient (2.6—33.0 psu) in a northern Norwegian fjord were 
investigated with a genome scan approach using both anonymously and 
EST-linked microsatellite loci (Coyer et al., 2011b). Six genes revealed 
putative signatures of selection, of which one was annotated to a putative 
mannitol transporter (Coyer et al, 2011b). As mannitol is regulated in 
highly saline environments to control cell turgor (Iwamoto and Shiraiwa, 
2005, and references therein), significant selection has occurred for at least 
some part of the mannitol pathway over the 12-km scale. Disentangling 
the dual roles of mannitol in osmoregulation and energy storage will 
require manipulative mesocosm experiments, which have not yet been 
performed. 

Intertidal gradients present several physiological stresses and it is along 
these gradients that classical transplant experiments have been conducted 
with several species of Fucus (see Section 4.1 and Table 2.1). Using an 
EEG approach, Pearson et al. (2009) compared emersion stress at the 
distributional centre versus the distributional edges in two Fucus species. 
Southern edge populations of F. serratus (northern Portugal) were less 
resilient to desiccation and heat-shock than central (UK) populations 
and expression of heat-shock genes was higher at the same temperatures, 
suggesting greater cellular stress (Pearson et al, 2009). In contrast, F. vesi- 
culosus showed no such divergence in heat-shock response and little 
variation in gene expression. These results suggest that for F. serratus, edge 
populations are maladapted, less genetically diverse and display lower 
fitness. Thus, changing climate conditions may threaten small, fragmen¬ 
ted and/ or marginal populations because of inherently reduced fitness and 
lower adaptive capacity. 

Using comparative gene expression studies of several heat-shock pro¬ 
teins in naturally occurring, sympatric populations ofBaltic Sea F. radicans 
and closely related F. vesiculosus, Lago-Leston et al. (2010) were able to 
show that F. radicans was more sensitive to mild heat-shock than F. vesicu¬ 
losus. Also from earlier experiments, the tolerance of F. vesiculosus did 
not vary substantially in contrast to its physiological variability in desic¬ 
cation and freezing tolerance (Pearson et al., 2000). At present, the inter¬ 
pretation of these results remains general, but suggests that changes in 
transcriptional regulation, combined with local ecological dynamics, may 
significantly impact functional traits in marginal habitats, coincident with 
genome-wide reductions in genetic diversity due to genetic drift (Pearson 
et al, 2009). 

At the level of biotic stress (i.e. community genomics), comparative 
transcriptomics is still in its infancy with the exception of host—parasite 
and plant—herbivore studies conducted mainly on model systems (e.g. 
Walley and Dehesh, 2010). 
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In order to explore the genetic basis of inducible chemical defence 
against herbivores in F. vesiciilosus, a cDNA microarray has been developed 
from libraries based on grazed and non-grazed F. vesiciilosus (Weinberger 
et al., unpublished). F. vesiciilosus is capable of inducing a chemical defence 
when exposed to grazers (Rohde et al, 2004). Both chemical elicitators 
(oligoguluronate, methyl-jasmonate) and direct feeding by isopod grazers 
altered gene expression. 

While Fucus has not yet entered the post-genome era, the rapidity 
of advances in next generation sequencing technologies (Metzker, 2010) 
is paving the way for the types of multi-level experiments essential to 
understanding complex ecological interactions involving stress (Gilad 
et al, 2007). Importantly, the demands for multi-factorial experimental 
designs and sufficient replication are already accessible and the costs are 
coming within reach of individual investigator budgets. 



. Conclusions 


The simultaneous effects of different stresses and the dynamics of 
positive and/or negative feed-backs form the selective environment on 
the genetic potential to determine whether a given species can (1) tolerate 
a transient stress regime (plasticity, acclimation); (2) adapt to a permanent 
regime (ecotypic differentiation via selection) or (3) become locally extinct 
(unable to maintain positive population growth). As the combination of 
stresses is variable over space and time, the challenge for an individual is to 
cope with the spectrum of stresses occurring during its lifetime; the chal¬ 
lenge for the population is to adapt to shifting stress regimes within a few 
generations. Finally, in the context of climate change, abiotic stress-driven 
effects on biotic interactions are expected to shift stress impact from the 
population level to the community level by altering competitive hierarchies 
and feedbacks among trophic and other interactions. 

At this point, we have a reasonable understanding of the effects of 
individual stresses on Fucus spp. at the population level, but our knowl¬ 
edge of their higher order interactions remains rudimentary. 
Consequently, our ability to predict interactions and their effects on 
community maintenance/assembly declines precipitously as the number 
of stresses increases. The model in Fig. 2.5 is a hypothesis of the network 
of single stressor effects and their interdependencies that may generate 
various interactive effects as they are currently understood. However, 
rigorous testing of this network of system-level hypotheses will require 
a different approach that takes into account non-linear and stochastic 
dynamics (Section 7). 

Results from the meta-analysis (Fig. 2.4) show that although all com¬ 
bined effects were, on average, additive this was not necessarily the case for 
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Figure 2.5 Example of an interaction network of abiotic and biotic stresses on F. 
vesiculosus (based on Wahl et al., 2010; Weinberger et al ., 2011, Table 2.1). Black 
circles used as arrow heads indicate negative interactions, black diamonds indicate 
positive interactions. Dark grey shading indicates biotic variables and light grey 
shading indicates abiotic variables. Examples: shading weakens energy reserves (a) 
leading to weakened anti-fouling (AF) and anti-grazing (AG) defences, which 
enhances fouling (2) and grazing (2 r ) pressure. Both reduce (3, 3 r ) the 
photo synthetic area of the thallus (c), which amplifies (4) the energy shortage 
under low light conditions. Low light (1/) and temperature—stress (‘warming’, 5) 
reduce growth (b), which jeopardizes the alga’s ability to compensate for tissue lost 
to grazers, the activity of which is enhanced by temperature stress (5 7 ). Emersion 
reduces fouling and strengthens anti-fouling defences (7')- Eutrophication, in 
contrast, enhances fouling (S f ) and competition (8) by favouring filamentous algae 
relative to Fucus. 1-2-3-4 is one of several examples for a positive feedback loop, 
which may be amplified or dampened by other abiotic and biotic stresses. 

individual case studies. There, the interaction between two stresses was 
often synergistic or antagonistic contingent upon the environmental con¬ 
text. The statistical test of interaction effects (Fig. SI) suggested a 
similar result although a slight trend towards synergistic interactions was 
stronger and more evident given the range of variation in individual 
studies. We emphasize that the results from the meta-analysis must be 
interpreted with caution, as the number of published studies fulfilling all 
requirements (factorial design, adequate replication, controlled treat¬ 
ments) was small. 


7 . Prospects 

The many species of Fucus span 45° of latitude in the northern 
hemisphere with dominance in the North Atlantic. Collectively they 
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Table 2.3 Distribution of Fucus species (left) and regions with habitats (right) 
used in the meta-analysis 


F. distichus 

1 

Baltic Sea (subtidal) 

16 

F. evanescens 

2 

Eastern Atlantic (intertidal) 

15 

F. gardneri 

3 

Eastern Pacific (intertidal) 

1 

F. serratus 

9 

Western Atlantic (intertidal) 

1 

F. spiralis 

2 



F. vesiculosus 

26 




inhabit a broad suite of habitats: exposed outer coast to sheltered bays/ 
fjords; fully marine to nearly freshwater; rocky substrata to saltmarsh; high 
intertidal to shallow subtidal. 

Despite this diversity, however, experimental work on stress is 
highly skewed with >60% of the experiments having been conducted 
on F. vesiculosus in the Baltic Sea or eastern North Atlantic (Table 2.3) 
and >95% of all studies having been conducted in the mid-range of 
the biogeographic distribution of the genus. Clearly, a fuller understand¬ 
ing of how members of the genus have colonized and continue to 
successfully exploit northern hemisphere shores, requires greater 
emphasis on the less studied species, as well as greater emphasis on their 
performance across dines and along distributional edges. Indeed, 
studies at distributional edges are likely to offer insights on how Fucus 
(and other species) will respond to global climate change (Pearson 
et ah, 2009). 

7.1. Experiments and modelling 

Our call for research on all species of Fucus and over their respective clinal 
and distributional ranges also acknowledges that our growing capacity to 
integrate physiology, genetics and experimental ecology provides numer¬ 
ous avenues to understanding the role of stress. In methodological terms, 
this will require (1) the design of manipulative experiments in a factorial 
manner to address multivariate interactions, (2) the up-scaling of experi¬ 
ments by including several interacting groups and by allowing enough 
time for effects to spread from the individual to the community level, (3) 
comparison of stress impact on marginal versus core populations, (4) 
attention to genotypic—phenotypic links and (5) the use of multivariate 
modelling techniques to explore and analyse the network system. 
Factorial ecological experiments with appropriate replication for >3 
stress factors, however, are logistically difficult (if not impossible) in the 
field and are only slightly less daunting to establish and maintain in the 
laboratory (mesocosms). Thus, future studies should strategically acquire 
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data that are appropriate for complex analyses as it is precisely the mul¬ 
tivariate effects that need to be understood. 

Structural equation modelling (SEM) provides a framework for esti¬ 
mating causal effects through the study of path relations (reviewed in 
Grace, 2006; Grace et al, 2010). Because the focus of SEMs is on under¬ 
standing direct and indirect pathways, the approach is well suited for both 
building and testing hypotheses of multiple stress processes (e.g. Tonsor 
and Scheiner, 2007, Arabidopsis ). In a marine setting, SEM has been used 
to model indirect facilitation of the benthic communities of kelp forests 
(Arkema et al, 2009). Its application to the hypothetical interaction 
network proposed in Fig. 2.5 has not yet been attempted but would 
greatly contribute to our understanding of community stress. 

7.2. Next generation molecular ecology 

Ecologists as a group tend to shun molecular methods in their research 
programs (Johnson et al, 2009). Nevertheless, the capacity to examine 
interactions among stress reactions—almost all of which are complex, 
multi-gene-regulated traits—will open fundamentally new avenues in 
community genetics—genomics. At this stage, we still know very little 
about the genetic basis of osmoregulation, thermal tolerance, chemical 
defence and the boundaries between plastic and selective responses for any 
of these traits. 

As distinguishing phenotypic plasticity from ecotypic differentiation 
can be approached directly by examining gene expression/regulation via 
differential transcription and structural changes in the DNA sequences of 
key genes, it is increasingly possible to integrate our understanding of 
ecological stress with complex genetic responses (Tautz et al, 2010). 
Thus, integrating the genetic perspective into field studies of stress 
ecology will enable us to (1) more finely identify plasticity in stress responses 
and expression of stress-related genes; (2) evaluate the adaptive potential of 
populations and (3) investigate how selection driven by one stress affects 
sensitivity towards other stresses. For example, differential mortality in 
response to one stress will lead to the survival of a non-random (selected) 
sub-sample of the original genotypic diversity, which in turn will harbour a 
sub-sample of traits ancillary to tolerance for the selective stress (e.g. tol¬ 
erance to a subsequent stress, productivity, light harvesting, etc.). Traits of 
the genotypic sub-sample may be non-random if the traits are linked to the 
tolerance of the first stress, or random if they are not. 


7.3. How will climate change affect fucus? 

As a foundational taxon, the genus Fucus has a large influence on the 
associated community as a whole; likewise, the associated community 
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influences Fuchs. This duality highlights the importance of gaining a more 
integrative understanding of stress responses in species of Fucus and the 
communities they support. 

Species with wide ecological niches are more easily able to cope 
with change; likewise, the ongoing and rapid diversification of the 
genus creates multiple overlapping niches. This suggests that Fucus 
spp. from stressful and fluctuating habitats may be better equipped 
to tolerate (via pre-adaptation) the superimposition of increased 
fluctuations (e.g. Schneider, 2008) than species from habitats with 
more subtle fluctuations (e.g. the deep subtidal, or low latitude habitats). 
As the predicted changes in the North Atlantic region over the next one or 
two centuries (3—6 °C increase in temperature, a 0.5 units decline in pH, 
IPCC Climate Change, 2007) are small relative to the daily environmental 
fluctuations characterizing any temperate northern hemisphere intertidal 
habitat, the direct impact of rapid climate changes on Fucus may be small. 

In contrast, the indirect impacts could be large if Fucus spp. were to be 
competitively displaced (as, e.g. by de novo interactions with introduced 
competitors, consumers and/or parasites and pathogens). However, the 
niche breadth of several sympatrically occurring Fucus spp. taken together, 
provides further insurance in terms of holding space and by ensuring its 
own ‘invasive’ potential (Peterson, 2003). As more than a third of Fucus 
spp. have been characterized as successful invaders (Williams and Smith, 
2007), the genus has high expansion potential, be it by anthropogenic or 
natural means. Invasiveness, that is, the ability of a species to establish in 
new habitats, might be a proxy for its capacity to tolerate climate change. 
Reproductive success in a new environment is favoured by phenotypic 
plasticity and/or rapid ecotypic selection. Even if the effects of climate 
change are disproportionally strong on marginal habitats (Melzner et al., 
2010; Thomsen and Melzner, 2010), Fucus spp. certainly represent prom¬ 
ising candidates to colonize areas abandoned by more sensitive 
macroalgae. 

The scale and importance of indirect effects of stresses, however, are 
only beginning to be understood. As a complicated yet realistic hypo¬ 
thetical example, moderate warming may not affect Fucus directly, but 
could have a large indirect impact via a synergistic interaction involving (i) 
a slight increase in the metabolic rate of consumers (^increased grazing 
pressure); (ii) a slight decrease in the defence strength of Fucus leading to; 
(iii) slightly higher fouling pressure and, in turn leading to; (iv) a slight 
decrease in available solar energy, further reducing; (v) defence strength 
and (vi) overall fitness (Fig. 2.5). The key realization is that accumulation 
of small effects acting in a non-linear manner could suddenly move 
towards a tipping point and major regime shift (Scheffer and van Nes, 
2004; Christensen et al., 2006; Harley et al., 2006). 
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7.4. Open questions 

Stress ecology has gained a new dimension of relevance and urgency in 
our contemporary world. At the same time we are in a position to 
integrate disciplines and exploit new types of data and multivariate anal¬ 
yses to address the role of climate change at the species and community 
level. With these considerations in mind, the following questions are now 
tractable. 

• Are there consistent combinations of stress interactions that may pro¬ 
vide general predictions about community responses to climate change? 

• What suites of genes are associated with particular stress interaction 
networks affecting, for example, salinity or temperature tolerance, 
chemical defence and susceptibility to predation? 

• What is the potential of biotic interactions to modulate the effects of 
abiotic stress and community resources; and how does this feedback to 
the long-term resilience of Fucus communities? 

• Will marginal and core populations differ in their response to global 
change? 

• What specific insights about climate change effects can be gained from 
understanding the responses of Fucus -based communities as compared 
with other foundational species and their associated communities (e.g. 
seagrasses, mussels)? 


Appendix A 

Supplementary Materials 


Supplementary data associated with this chapter can be found, in the 
online version, at doi:10.1016/B978-0-12-385536-7.00002-9. 
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Abstract 

The physiological mechanisms that regulate adaptive plasticity of clonal 
organisms are key to their success in changing environments. Here, we 
review the mechanisms that regulate morphological plasticity of colonial 
hydrozoans. There is a heritable, genetic basis to colony form, but 
environmentally-induced plasticity and self-reinforcing developmental 
physiology explain much of total phenotypic variance. Morphological devel¬ 
opment of colonial hydrozoans emerges from interactions among (1) beha¬ 
viors which drive gastrovascular transport, (2) architecture of the 
gastrovascular system that determines hydrodynamic characteristics of vas¬ 
cular flow, and, (3) gene products that vary in response to physiological 
signals provided by gastrovascular transport. Several morphogenetic signal¬ 
ing mechanisms have been identified, including, reactive oxygen species and 
nutrient concentrations in the hydroplasm, and hydromechanical forces asso¬ 
ciated with gastrovascular transport. We present a conceptual model of the 
interacting forces that drive hydrozoan morphological development. Several 
avenues for future research are suggested by the synthesis of information 
from prior studies of hydrozoans. Elucidating the morphogenetic signaling 
pathways responsive to metabolites or hydromechanical forces and the epi¬ 
genetic effect of vascular architecture on colony form may give new insight 
into the self-maintenance of indeterminately growing and continuously devel¬ 
oping vascular systems. 



Introduction 


Clonal organisms grow and propagate by self-replication of genet¬ 
ically identical units (i.e. ramets) that acquire resources, provide structure 
or reproduce and can survive on their own if separated from one another 
(Jackson, 1985). Multicellular clonal (or colonial) organisms occur in 
many lineages of the Eukarya including several invertebrate phyla, fungi, 
benthic macroalgae and plants. Ecologically, they cover much of the 
biosphere and are particularly common in marine benthic communities 
(Boardman etal., 1973; Jackson, 1979, 1985; Buss and Blackstone, 1991). 
Clonal organisms may occur as solitary individuals or be physically con¬ 
nected and physiologically integrated, often by means of a shared vascular 
system (Jackson, 1985; Harvell and Grosberg, 1988; Buss and Blackstone, 
1991; Blackstone and Bridge, 2005). Here, we examine the scope for 
alternative patterning of those connections. 

Besides the phylogenetic diversity and ecological success of clonal 
lineages of eukaryotes, there are two other compelling reasons to study 
the morphological patterning of colonial organisms. First, longevity, size 
and shape of a clone can vary enormously such that a developmental 
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programme cannot pre-determine the timing and location of every mor¬ 
phogenetic event during the life span of the organism. This contrasts with 
the model systems of development (unitary, [aclonal] organisms; e.g. fruit 
flies, mice, worms) in which gene expression early in life results in form 
that thereafter grows. Second, different forms show different ecological 
traits and patterns of selection on whole colony morphology (Buss and 
Grosberg, 1990; Yund, 1991). Despite the lack of developmental genetic 
programmes to specify a finished form, coherent and repeatable colony 
morphologies are nevertheless generated from the physiological state of 
the organism interacting with the environment to regulate gene expres¬ 
sion. We review what is known about the development and maintenance 
of morphology in Hydrozoa (Phylum Cnidaria) with respect to herita- 
bility of colony form and the extent of phenotypic plasticity. We discuss 
existing evidence and present data supporting hypothesized physiological 
mechanisms of plasticity and constraint on plasticity imposed by colony 
architecture. 

Synthesizing the above, we present a conceptual model describing 
colony growth and development as a context-dependent dynamical sys¬ 
tem of physiological responses to environmental cues. The outcome of 
colony development from these responses is contingent upon the existing 
colony architecture as well as genotype. In this way, feedback mechanisms 
and environmental signals operating during development enable a self- 
maintaining system leading to one of multiple potential morphological 
states for a given genotype. Future directions for research are offered based 
on (1) unresolved questions and (2) the opportunity to use indetermi¬ 
nately growing species as models of internal vascular systems that develop 
continuously in an organism’s internal environment. 

1.1, Hydrozoans as models for studying physiological 
regulation of colony development 

Clonal organisms differ markedly from aclonal (or unitary) organisms 
with respect to their morphological construction, ecological traits and 
patterns of growth (Harper, 1977; Jackson, 1979, 1985; Buss, 1987; 
Harvell and Grosberg, 1988; Caswell, 2001; Winston, 2010). These 
fundamental differences in the morphology and ecology of clonal and 
aclonal organisms derive primarily from differences in their respective 
patterns of development. In contrast to unitary organisms in which devel¬ 
opment of the adult form occurs during a relatively narrow window of 
time early in ontogeny, development of clonal organisms continues with 
growth by means of active somatic stem cell lineages (Buss, 1987; 
Blackstone and Bridge, 2005). Variation in the relative growth rates of 
different modules determines their spacing and arrangement (i.e. form). 
Growth may be indeterminate and, consequently, morphogenetic 
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patterning occurs over a wider span of the lifetime. Continuous devel¬ 
opment can be adaptive in organisms with a sessile habit and great lon¬ 
gevity because they likely face environmental challenges they cannot 
escape during their life span. A developmental programme that responds 
to environmental cues (i.e. environmental signalling, sensu Blackstone and 
Bridge, 2005) produces morphologically plastic phenotypes that amelio¬ 
rate unfavourable conditions (Bradshaw, 1965; Thompson, 1991; Stearns, 
1992; Blackstone, 2003). 

Morphological plasticity is widespread among clonal organisms and its 
ecological significance clearly established (Foster, 1979, 1980; Palumbi, 
1984; Schemer and Goodnight, 1984; Harvell, 1986; Stearns, 1992; 
Scheiner, 1993a; Schmitt, 1993). In the past two decades, research on 
colonial animals has emphasized physiological processes that regulate 
morphogenesis as colonies develop and features that limit the extent of 
morphological plasticity. Colonial hydrozoans (Phylum Cnidaria), espe¬ 
cially those in the Family Hydractiniidae, are useful models of the phys¬ 
iological regulation of development. Hydrozoan colonies have relatively 
simple morphology, consisting of polyps (modules for feeding, transport¬ 
ing gastrovascular fluid, defence and reproduction) and stolons (network 
of vascular tubes; Fig. 3.1). The principal colony-wide physiological 
system is the gastrovascular system, fluid-filled canals that run throughout 
the stolonal network and join the gastric cavities of polyps (Berrill, 1949; 
Crowell, 1957;Mackie, 1973, 1986). The gastrovascular system performs 
most of the extracellular physiological functions of the colony including 
digestion, transport of food and metabolites and waste removal (Lenhoff, 
1968; Rees etal., 1970; Murdock, 1978; Gladfelter, 1983; Blackstone, 
1999, 2001, 2003). Hydrozoans are useful in experimental manipulation 



Figure J.l Hydrozoan colonies growing on glass slides showing principal elements 
of a colony. (A) H. symhiolongicarpits characterized by feeding polyps (gastrozooids), 
the ectodermally fused stolonal mat and peripheral stolons extending beyond the 
mat. (B) P. earned, which lacks a stolon mat, displays polyps atop sparsely branched 
stolons. White scale bars at the bottom of each photo are 1 mm. 
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of colony form, being effectively an exposed vascular system unencum¬ 
bered by organs and tissues. 

Variation in spacing and arrangement of polyps and stolons is partially 
heritable and representative of the range of morphologies observed 
among diverse clonal taxa (Buss and Blackstone, 1991; Dudgeon and 
Buss, 1996). Morphological variants have been shown to correlate with 
ecological traits (Buss, 1979; Jackson, 1979; Yund, 1987, 1991; Black- 
stone and Yund, 1989; Buss and Grosberg, 1990). In addition to heritable 
variation of morphology (McFadden et al ., 1984; Blackstone and Buss, 
1991), colonies can display remarkable plasticity in growth and form 
(Crowell, 1957; Buss and Blackstone, 1991; Blackstone and Buss, 1992, 
1993; Dudgeon and Buss, 1996). Knowledge of the environmental sig¬ 
nalling that underlies morphological plasticity may provide new insight 
into development not afforded by study of the gross morphology of 
unitary model organisms but directly analogous to the indeterminate 
development of vascular systems in both colonial and unitary organisms. 
Finally, including colonial organisms as model systems for studying devel¬ 
opment broadens perspective of the conceptual framework of develop¬ 
ment as a dynamical system of interaction between gene expression, the 
modular features of organisms and the environment. Such a perspective 
may complement emerging views of the modularity of development of all 
multicellular organisms (Schlosser and Wagner, 2004). 


2. Experimental Studies of Hydrozoan Colony Form 

2.i. Genetic basis of colony form 

The spacing and arrangement of polyps and stolons varies widely both 
among and within cnidarian species. At one end of the spectrum of 
morphological variation, colonies (called ‘runners’) show widely spaced 
polyps along long, sparsely branched stolonal connections. At the other 
end, colonies (called ‘sheets’) show closely packed polyps along short, 
densely branched and anastomosed stolons (Buss and Blackstone, 1991). A 
genetic basis to colony form is suggested by the fact that colony form is 
routinely used to distinguish among ecologically distinct morphotypes 
within species or, identify cnidarian species in situ (e.g. Knowlton etal., 
1992) and a clone-specific proportion between stolon length and polyp 
number (Braverman, 1963), such that replicate explants from the same 
colony produce similar outcomes of morphology (McFadden et al., 1984). 
This is remarkable given the nature of indeterminate growth and variable 
environments that can influence growth rates of colony parts as well as the 
size, shape or longevity of the colony. To what extent is cnidarian colony 
form under direct genetic control? 
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For most colonial invertebrates, the respective contributions of heri¬ 
table genetic variation and environmentally determined variation to phe¬ 
notypes are unknown. Morphological plasticity has been observed in 
laboratory and field studies of cnidarian species (Crowell, 1957; Dustan, 
1975; Foster, 1979, 1980; Graus and McIntyre, 1982; Brown etal., 1985; 
Willis, 1985; Sebens and Done, 1992; Bruno and Edmunds, 1997) and so, 
environmental variance (V E ) (see Table 3.1 for a summary of all abbrevia¬ 
tions) is expected to represent a large proportion of total phenotypic 
variance of colony form. 

Attempts to partition phenotypic variance of colony form into genetic 
and environmental components have been made for a few species of 
hydrozoans. The approach used most frequently is growing clonal repli¬ 
cates of several genotypes in a common garden and measuring clonal 
repeatability, r. This approach does not technically distinguish genetic 
and environmental components because shared environmental history 
among clonal replicates may overestimate the genetic component 
(Falconer, 1986; Lynch and Walsh, 1998). Nevertheless, clonal repeat¬ 
ability should represent an upper estimate of broad-sense heritability, 
(the sum of additive and non-additive genetic variance). Values of 


Table 3.1 Abbreviations for genetic effects and their descriptions 


Component 

description 

v P 

Total phenotypic variance of a trait; Vp = Vq + V E 

v G 

Measures the amount of phenotypic variation among averages 
of different genotypes; V G = V A + Fna 

V E 

Component of phenotypic variance attributable to 

environmental variation measured among individuals with 
the same genotype 

V A 

Additive genetic variance; the linear component of 
phenotypic change associated with allelic substitution 

I 7 na 

Non-additive genetic variance; the non-linear component of 
phenotypic change associated with allelic substitution 
caused by dominance and/or epistasis 

r 

Clonal repeatability; the ratio of the between-individual 
component to the total phenotypic variance, which 
includes both genetic and environmental components, the 
latter due to non-localized, permanent or historical effects; 
r=(V G +V Eg )/V p 

hi 

Broad-sense heritability; the fraction of phenotypic variance 
attributable to differences in genotype; hi = Vc/Vp 

h N 

Narrow-sense heritability; the fraction of phenotypic variance 
attributable to additive genetic variance; h ^ = Vj\/Vp 
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heritability range from 0 to 1. Those approaching 1 indicate a large 
genetic component, and those near 0 indicate a large environmental effect 
on morphology. 

Common-garden studies of clonal repeatability have been used to 
distinguish colony form among species and genotypes within species of 
Podocoryne carnea and four species of Hydractinia. In some cases, only 
statistics establishing differences in colony form among species or geno¬ 
types are given implying genetic sources of variation, but not estimates of 
repeatability (McFadden etal., 1984; Blackstone and Buss, 1991). For 
studies that present such estimates, h ^ from measures of clonal repeatabil¬ 
ity of different Hydractinia species are typically quite small for overall 
colony morphology (Table 3.2). Values of repeatability range from 0.17 
to 0.31 in Hydractinia polyclina, 0.21 in Hydractinia symbiolongicarpus and 
0.24 in Hydractinia [GM] (an unnamed Hydractinia species in the Gulf of 
Mexico; D. Ferrell, personal communication). These low values of 
repeatability, representing an upper estimate of lig, imply that most var¬ 
iation in colony form is environmentally determined. Yund (1991) 
observed values of clonal repeatability in H. polyclina of 0.48 and 0.58 
for stolonal mat and peripheral stolon production, respectively, somewhat 
greater than those measured for overall colony form. 

Blackstone and Buss (1991) provided an important addition to the 
studies that report values of clonal repeatability. In addition to common- 
garden experiments (although repeatability was not estimated), it is the 
only breeding study of Cnidaria that has attempted to quantify the her¬ 
itable genetic components of colony form. They mated two males to each 
of five females to analyse co-variation among both full- and half-siblings 
(average number of offspring per family = 24 after 50 days). In a nested 
analysis of variance, the among-male component measures /ig, whereas 
the among-female component, which estimates covariance among half¬ 
siblings, measures additive genetic variance in order to estimate narrow- 
sense heritability (h^). All of their analyses indicated a strong effect of the 
male parent and a weak effect of the female parent, which correspond to a 
large non-additive genetic variance component and a small additive 
genetic variance component. Because of the small number of families 
and offspring in the study they focused on qualitative interpretations of 
the analysis, rather than exact calculations of heritability. Blackstone and 
Buss (1991) inferred from these results that genetic control of colony form 
was manifested principally from dominance and epistasis. 

For comparisons with other estimates of from measures of clonal 
repeatability, we estimated broad- and narrow-sense heritabilities from 
the data presented in Blackstone and Buss (1991) (we assumed common 
environmental variance = residual environmental variance, as is sup¬ 
ported by the lack of position effects in their analyses). From their data, 
we estimated h ^ = 0.84, which was comprised largely of non-additive 




Table 3.2 Measures of clonal repeatability as upper estimates of broad-sense heritability in Hydractinia spp. 


Species 

Analysis 

Variable 

r 

/jg Study 

H. polyclina 

Repeatability 

Peripheral/mat 

0.31 

Yund (1991) 

H. polyclina 

Repeatability 

Stolon tips 

0.17 

D. Ferrell (personal communication) 

H. symbiolongicarpus 

Repeatability 

Stolon tips 

0.21 

D. Ferrell (personal communication) 

H.[GM] 

Repeatability 

Stolon tips 

0.24 

D. Ferrell (personal communication) 

H. symbiolongicarpus 

Half-sibling breeding 

Perimeter/area 172 


0.84 Blackstone and Buss (1991) 


Note: Other similar studies did not report values for r, or h g, and they could not be inferred from the data presented. 


Dudgeon and Kiibler 







The shape of things to come 


115 


genetic variance (dominance and epistasis, F NA = 0.69; V A — 0.15), and 
environmental variance comprised only 16% of phenotypic variance. The 
breeding study by Blackstone and Buss (1991) suggests a much greater 
genetic component underlying colony form of hydractiniid hydrozoans 
than do other studies measuring clonal repeatability. Unfortunately, clonal 
repeatability for neither species can be calculated from the data provided 
in Blackstone and Buss (1991) for direct comparison. 

The different interpretation of the importance of genetic versus envi¬ 
ronmental control of colony form has at least two potential explanations. 
Estimates of clonal repeatability were made on colonies freshly collected 
from the field, whereas colonies used in breeding experiments were cul¬ 
tured in the laboratory until reproductively mature (2 months for 
Hydmctinia ; Blackstone and Buss, 1991). The difference may also reflect 
different morphometric variables used in analyses among studies by differ¬ 
ent investigators. Morphological variation in clonal organisms is often 
described qualitatively in terms of the suites of traits characteristic of 
different forms (Buss, 1979; Jackson, 1979; Lovett-Doust, 1981; Harper, 
1985; Buss and Blackstone, 1991). Consequently, there are several potential 
descriptor variables that may effectively represent morphology to varying 
degrees, but may differ in the strengths of the genetic bases underlying the 
respective traits. There is also the possibility that the colonies used in various 
experiments sample different amounts of the genetic scope for plasticity in 
the population at large. Reviewing multiple studies, we can only conclude 
that there is strong evidence for interrelated genetic and environmental 
determinants of clonal morphology of hydrozoans. 

The correlation between ecological traits of interest (e.g. competitive 
ability, demographic parameters) and colony morphology in hydractiniid 
hydrozoans, specifically, and clonal organisms in general, emphasizes the 
genetic basis of form (Buss and Grosberg, 1990; Yund, 1991; Bruno and 
Edmunds, 1997). In particular, traits like competitive ability, for which 
there appears to be an underlying morphogenetic basis can only be 
efficaciously acted upon by selection if much of the genetic basis to colony 
form is in the form of additive genetic variance. The present evidence is 
inconclusive as to the extent to which colony form is heritable, but 
available data suggest that little of the variation in morphology is attrib¬ 
utable to additive genetic variance. 

2.2. Colony integration by means of the gastrovascular system 

In the simply constructed cnidarians, the physiological system regulating 
morphological development is the gastrovascular system, which carries out 
much of extracellular physiological function of the colony. The hypothe¬ 
sized role of the gastrovascular system in colony-level patterning means that 
we need to understand the dynamics of gastrovascular transport. 
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The gastrovascular system integrates polyps of a colony with respect 
to sharing of resources, coordinating behaviours and defensive 
responses. Several early studies recognized the importance of gastrovas¬ 
cular transport to colony growth (Crowell, 1957, 1974; Hale, 1960; 
Wyttenbach, 1968, 1973; Karlson and Marferin, 1984; Marferin, 
1985; Beloussov etai, 1989). It is not certain whether polyps remain 
locally autonomous, or whether they serve as organizing centres initi¬ 
ating behaviours (Mackie, 1986). The role of polyps may vary with 
taxon. Parrin etal. (2010) summarize the behavioural control of cnidar- 
ian gastrovascular dynamics. Briefly, muscular contractions of polyps 
outweigh ciliary action in driving gastrovascular flow in hydrozoans 
and the derived anthozoan group, the pennatulacean octocorals 
(Parker, 1920; Berrill, 1949; Crowell, 1957; Hale, 1960; Fulton, 
1963; Wyttenbach, 1968; Brafield, 1969; Schierwater et al, 1992; 
Blackstone, 1996; Wagner etal., 1998; Dudgeon etal., 1999; Anctil 
etal., 2005; Parrin etal., 2010; but see Roosen-Runge, 1967 for an 
exception among Hydrozoa). Studies of other octocorals and hexacorals 
indicate that gastrovascular transport is achieved by cilia-driven flow, 
which may be the ancestral state for the cnidarian gastrovascular system 
(Murdock, 1978; Gladfelter, 1983; Schlichter, 1991; Gateno etal, 1998; 
Parrin etal ., 2010). 

Behaviours of polyps and stolons driving the dynamics of gastrovascu¬ 
lar flow have been quantified in vivo in hydrozoans. Polyp behaviour 
can be characterized by the differing input—output relationships between 
the polyp and either the external (via the mouth) or the internal (via the 
polyp—stolon junction) environment, and delineated with respect to the 
timing of ingestion of food (Wagner et al., 1998; Dudgeon, 2001a). Polyps 
prior to feeding, or long after regurgitation show infrequent contractions, 
which occur without periodicity, and the volume of fluid exchanged with 
the stolon is extremely small (Dudgeon etal, 1999; Dudgeon, 2001b). 
Fed polyps initiate the first phase of oscillatory behaviour (expansion- 
contraction cycles of polyps with respect to their length and width, hence 
volume) within 5—15 min of ingestion. In the second phase of oscillatory 
behaviour, polyps behave semi-autonomously and display constant ampli¬ 
tude oscillations. Digestion occurs in the gastric cavity of the polyp during 
this stage, but small volumes may be exported to the common gastrovascular 
system. The final, third phase is initiated by the export of dense streams of 
particulates from the polyp into the stolon and terminated by regurgitation 
of undigested food, which triggers a return to pre-feeding behaviour. This 
phase of maximal transport rates occurs between 2 and 8 h after feeding 
(Schierwater etal, 1992), but tends to occur later and for longer durations 
in single polyp systems, or very small colonies (Dudgeon etal, 1999). 

Direct evidence of cause and effect for polyps moving gastrovascular 
fluid through stolons, first in one direction, then the other, would be 
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alternating oscillations of polyps and stolon lumens. In a simple colony 
with a single polyp and known length of stolon, time series of changes in 
polyp volume correspond with changes in lumen diameter of the stolon 
(Fig. 3.2A). Stolon diameter is at a minimum when polyp volume is 
maximal and stolon diameter increases as polyp volume declines, the 
former reaching its maximum diameter when polyp volume is minimal. 
Since fluid moves in both directions in a stolon (to and from a polyp), it 
must go from being at rest to a maximal velocity in one direction and then 
decline to 0 velocity as it switches to moving in the other direction. 
Flowspeed is greatest immediately upon expansion of the stolon lumen 
as polyp volume decreases from maximum to minimum volume 
(Fig. 3.2B, see also Van Winkle and Blackstone, 1997). Flowspeeds in 
the stolon approach 0 when polyp volume is at, both, maxima and 
minima. Flowspeeds in the stolon are greater when moving away from 
the polyp (during export by polyps) than when moving towards the polyp 
(during import by polyps). The volumetric rate of fluid transport in the 
stolon follows the dynamics of flowspeed (Fig. 3.2C). Volumetric rate is 
maximal as the stolon lumen rapidly opens and expands during export by 
polyps. As the lumen approaches its maximum diameter, volumetric rate 
drops quickly and approaches 0 as the stolon is maximally expanded. 
Volumetric rate again increases as it moves towards the polyp upon closing 
of the stolon lumen, but the rate of transport back towards polyps is less 
and the duration slightly longer than during export. 

The dynamics of polyp behaviour and gastrovascular transport change 
dramatically following ingestion and continue to change through regur¬ 
gitation. Wagner etal. (1998) identified alternative hypotheses to explain 
the dynamics observed and developed an ordinary differential equation 
model to test whether the hypothesis of nutrient generation (from diges¬ 
tion) and metabolism could produce model polyp behaviours iti silico 
qualitatively similar to observed polyp behaviours following ingestion 
of food and during digestion. Modelled polyp behaviour showed a 
start-up of oscillatory behaviour, increased polyp length and oscillation 
amplitude during digestion and a steady state of oscillations until nutrients 
were depleted (Wagner et al ., 1998). The model equations (differentiated 
with respect to time) represented (1) the quantity of food remaining in the 
polyp gastric cavity (a monotonically declining function until the next 
ingestion period), (2) the nutrients available to the polyp, which are 
generated by the digestion of food and depleted by metabolism and export 
from the polyp and (3) a function in which two state variables of the polyp, 
representing an underlying non-linear oscillator, are sensitive to the quan¬ 
tity of nutrients in the gastric cavity. Oscillations in values of these state 
variables are reflected in oscillatory behaviour by polyps. The alternative 
hypothesis that mechanical strain on polyp tissue caused by ingestion of 
food elicits polyp oscillations is compromised by observations that 
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Figure 3*2 Gastrovascular dynamics in P. carnea 3—5 h after feeding. Data represent 
behaviour of a single polyp connected to a stolon network with a volume capacity in 
excess of 5 nl. Dynamics over a ~630 s interval showing oscillations of polyp volume 
(solid circles and line) and lumen diameter (dashed line) of stolons (A), polyp volume 
(solid circles and line) and flowspeed (dashed line) of gastrovascular fluid through 
stolons (B) and lumen diameter of stolon (dashed line) and volumetric rate (solid 
circles and line) of gastrovascular fluid transport (C). Volume estimates of polyps 
were sampled every 8 s. Stolon dynamics (lumen diameter, flowspeed and 
volumetric rate) were sampled at 30 frames s 1 . Polyp contractions export fluid 
into stolons which open rapidly to near their maximum diameter as the fluid 
accelerates to maximum speed and then declines to 0 velocity before switching 
direction as the polyp begins to expand. 
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oscillations may not begin for up to 15 min after strain is initiated upon 
ingestion and lacks a mechanism to induce regurgitation by polyps. The 
hypothesis that nutrients emanating from food substrates elicit oscillatory 
behaviour by polyps insilico, including a rise in oscillation amplitude to a 
steady state as nutrient concentration plateaus and the eventual return to a 
quiescent state when nutrients drop below levels that elicit polyp oscilla¬ 
tions, is consistent with observed in vivo polyp behaviours. 

As expected, quantitative features of oscillating polyps were influenced 
by environmental conditions. Hydrozoan polyps oscillated with greater 
frequency with increasing temperature. At temperatures of 34 and 35 °C, 
near the upper thermal limit of H. symbiolongicarpus, polyp behaviours 
underwent a transition in their behaviour from stable limit cycles to 
intermittent and erratic changes in oscillation frequency, to a chaotic 
regime of oscillations in which oscillation frequencies varied widely for 
indefinite periods of time (Buss and Vaisnys, 1993). Chaotic polyp beha¬ 
viours were reversible by lowering temperature upon which stable limit 
cycle behaviours resumed. 

Quantitative models that characterize the dynamics of gastrovascular 
transport have time-varying behaviour of oscillations in polyp length or 
volume as output, not a morphogenetic pattern. Pumping behaviour of 
polyps is still not mechanistically linked to the expression of colony 
morphology. Conceptual models to explain colony-level morphogenesis 
and its plasticity (described below) invoke patterns and rates of gastrovas¬ 
cular transport, or substances transported within it, during the phase of 
maximal gastrovascular activity. 

2.3. Gastrovascular physiology and colony form 

Patterns of fluid transport through the gastrovascular system emerge from 
the collective behaviour of polyps interacting with one another and the 
environment. Thus, physiological states of the gastrovascular system can 
represent environmental sensors from which adaptive morphological 
plasticity develops. For instance, Blackstone (2003) pointed out that 
developmental pathways that are sensitive to environmental conditions 
and link those cues to timing the action of morphogenetic mechanisms 
(i.e. initiation of polyps and stolon branches relative to stolon elongation) 
will be favoured by selection. Among modular organisms displaying 
continuous development, genetic information relevant to colony-level 
patterning has been hypothesized to represent a series of rules iteratively 
applied in response to the dynamics of the colony’s physiological state 
(Ulam, 1962; Braverman and Schrandt, 1966; Buss, 2001; Kaandorp and 
Kiibler, 2001). Essentially the same idea of colony-level patterning was 
called environmental signalling by Blackstone and Bridge (2005), which 
they defined as developmental pathways influenced by cues that initiate 
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outside the organism and underlie phenotypic plasticity. Both views share 
the idea that physiological mechanisms that can sense and respond to 
environmental cues are inextricably linked to morphological develop¬ 
ment of colonies. 

All evidence from hydractiniid hydrozoans indicates that gastrovascu- 
lar transport mediates morphological development of colonies, with, 
potentially multiple signalling mechanisms at work. Colonies that exhibit 
higher volumetric rates of flow (via larger diameter stolons and/or faster 
flowspeeds) exhibit more rapid elongation of stolons and lower rates of 
stolon branching and polyp bud formation (Blackstone, 1996, 1998; 
Dudgeon and Buss, 1996). For instance, colonies of H. symbiolongicarpus 
that display a more reticulate vasculature compared to P. earned show lower 
rates of flow than P. carnea at all developmental stages except at, and prior 
to, mat formation (very early in development; Blackstone, 1996). 
Colonies of P. carnea derived from inbred lines and selected for morphol¬ 
ogies with long, sparsely branched stolons show greater volumetric rates 
of vascular transport than those with shorter branches and anastomosed 
stolons (Blackstone, 1998). 

Given the correlation between gastrovascular flow rate and growth 
rates of modules, simple principles of fluid mechanics, therefore, pre¬ 
dict that the relative sizes and arrangement of stolons (i.e. vascular 
architecture) should influence subsequent colony development. The 
runner-to-sheet morphological continuum spans a wide range of 
colony form in Hydractinia spp. The vascular architecture of runner 
colonies of Hydractinia spp. is characterized by a highly variable size- 
frequency distribution of stolon diameters (Dudgeon and Buss, 1996). 
In contrast, the size—frequency distribution of diameters in the stolonal 
network of sheet colonies shows little variation. Stolon diameter is 
significant because volume flux and other hydrodynamic characteristics 
of transport vary with the dimensions of the stolon through which the 
fluid flows. For instance, in accord with the Law of LaPlace (Vogel, 
1988) wider stolons disproportionately increase their diameter com¬ 
pared to narrower stolons when pressurized. Moreover, at the low 
Reynolds numbers relevant here (10 “—10 3 ; Dudgeon and Buss, 
1996), wide stolons transport a much greater volume of fluid than 
narrow stolons for a given pressure drop and stolon length (Vogel, 
1994). Consequently, runner colonies have an unequal distribution of 
flow to the growing periphery, whereas sheet colonies have equivalent 
volumes of flow delivered to the growing tips. These flow patterns 
reinforce their respective growth trajectories, as in a self-maintaining 
system. 

In addition to observing co-variation between rates of transport and 
colony form, experimental evidence shows that manipulating gastrovas¬ 
cular transport causes changes in colony form. Diminished flow rates in 
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stolons, either by treatments that shift mitochondrial redox state towards 
oxidation or excess feeding typically causes elevated rates of stolon branch 
and polyp bud formation resulting in sheet morph colonies (Blackstone 
and Buss, 1992, 1993). Moreover, in P. carnea a dose—response relationship 
between gastrovascular flow to peripheral stolon tips and colony devel¬ 
opment was observed; the greater the perturbation to gastrovascular 
transport, the more colony patterning was altered (Blackstone, 1997). 

There are two alternative, but not mutually exclusive, mechanisms by 
which gastrovascular transport may signal morphogenetic events. 
Endodermal cells lining the lumen of the stolon vasculature may respond 
to chemical substances (e.g. metabolites, diffusing morphogens, nutrients, 
etc.) transported in the flow, forces directly caused by gastrovascular flow 
or both. 

Multiple, independent lines of evidence from P. carnea indicate that 
reactive oxygen species (ROS) signal morphogenetic events. In fact, an 
increased oxidative state alone within colonies ofP. carnea was sufficient to 
stimulate polyp bud and stolon tip formation in the absence of diminished 
flow (Blackstone and Buss, 1993; Blackstone, 1999, 2001, 2003). 
Treatment with pharmacological oxidizing agents, like 2,4-dinitrophenol 
(DNP) partially uncoupled oxidative phosphorylation diminishing ATP 
production and vascular transport while shifting colony redox state 
towards oxidation (Blackstone and Buss, 1992, 1993; Blackstone, 
2009). Relative oxidation of colony physiological state (via DNP or other 
oxidizing agents with different associated effects, like carbonyl cyanide m- 
chlorophenylhydrazone [CCCP]) increased rates of stolon branching and 
polyp formation (i.e. sheet growth; Blackstone, 1999, 2001, 2003). 
Treatment with pharmacological reducing agents, like sodium azide, 
blocked the electron transport chain, preventing 0 2 uptake, thereby 
diminishing ATP production and vascular transport and shifted colony 
redox state towards reduction (Blackstone, 2001, 2003, 2009). A shift to a 
more reduced physiological state correlated with lower rates of polyp bud 
and stolon tip formation (i.e. runner growth) compared to controls 
(Blackstone, 2001, 2003, 2009). 

ROS may participate in multiple signalling pathways from elongating 
stolons and inhibiting branching to cell death (Blackstone et al ., 2005). 
Not only do stolons extend and sprout branches to varying degrees but 
they can regress. Treatment of hydroid colonies (P. carnea and Eireneviri- 
dnla) with nitric oxide donors caused stolon regression (Cherry Vogt 
et al ., 2008) as did treatment of peripheral stolon tips with the chemi¬ 
cal antioxidant, vitamin C (Blackstone et al ., 2005). 

Perturbations to the redox state of a colony generate signatures of ROS 
that emanate most strongly from the epitheliomuscular cells (EMCs) near 
the base of polyps (Blackstone, 2003; Blackstone etal ., 2005). This obser¬ 
vation is consistent with polyp oscillations driving the transport of 
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gastrovascular fluid (Wagner etal, 1998; Dudgeon etal, 1999; Dudgeon, 
2001a,b, Bumann and Buss 2008) and provides an ecological explanation 
linking adaptive morphological plasticity to environmentally heteroge¬ 
neous food supplies via ROS signalling associated with polyp behaviour 
(Blackstone, 2001, 2003; Blackstone and Bridge, 2005). Polyps actively 
pumping following ingestion become relatively oxidized and ROS ema¬ 
nating from mitochondrion-rich EMCs and dispersed in gastrovascular 
fluid appear to signal increased rates of stolon branching and polyp for¬ 
mation to exploit locally abundant food supplies (Blackstone, 2001). 

Environmental signalling may operate through pathways in addition to 
ROS signalling in some hydrozoans. This is suggested by tests of different 
hypotheses as well as results from dinitrophenol experiments using H. 
symbiolongicarpus (Blackstone and Buss, 1993). DNP treatment shifted 
redox state towards oxidation but colonies showed reduced rates of bud 
and tip formation (i.e. runner-like growth). Developmental plasticity is 
sensitive to colony nutritional status, diffusing morphogens and vascular 
architecture, which determines volume flux and hydromechanical forces 
associated with gastrovascular transport (Turing, 1952; Crowell, 1957; 
Wyttenbach, 1974; Meinhardt, 1976; Muller and Plickert, 1982; Plickert 
etal., 1987; Lange andMiiller, 1991; Blackstone andBuss, 1993; Dudgeon 
and Buss, 1996; Blackstone, 2001; Dudgeon, 2001c). 

With respect to nutritional status, Bumann and Buss (2008) elaborated 
upon earlier works by Crowell (1957) and Braverman (1974) and found 
that food supply and the distribution of fed and unfed polyps markedly 
affected colony morphology. They suggested that ROS may not be the 
cause of morphogenetic patterning, but may simply reflect epitheliomus- 
cular activity of large, fed polyps generating local circulation patterns 
among themselves that biases the distribution of nutrients and, hence, 
zones of growth in a colony. 

The different vascular architectures of runner and sheet colonies of H. 
symbiolongicarpus were surgically manipulated to alter hydrodynamic char¬ 
acteristics of gastrovascular transport and assay subsequent trajectories of 
colony form (Dudgeon and Buss, 1996). Reduced variation in flow rate 
among stolons of runner colonies resulted in subsequent sheet growth. 
Conversely, increased variation in flow among stolons within sheet col¬ 
onies resulted in subsequent runner growth. Even without disrupting the 
colony-wide volume flow rate through metabolic perturbations, colony 
development was malleable in response to changed gastrovascular flow 
patterns. Based on these results, Dudgeon and Buss (1996) suggested that 
hydromechanical forces of gastrovascular fluid flowing through colonies 
act as morphogenetic signals, as occurs in vertebrate vascular systems 
(Davies etal., 1986; Nollert etal., 1990; Ziegelstein etal., 1992; Chen 
et al., 1999; Yeh et at., 1999; Osawa et at., 2002; Shay-Salit et at., 
2002; Urbich et al., 2002; Resnick etal., 2003; Milkiewicz et al., 2007). 
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Buss (2001) examined whether hydromechanical forces were plau¬ 
sible mechanisms regulating polyp formation and stolon branching 
based on the results of Dudgeon and Buss (1996) and observations by 
Crowell (1957) and Braverman (1963). Crowell (1957) noted that inter¬ 
polyp distances were greater in Campanulariaflexuosa colonies with wider 
stolon radii and Braverman (1963) observed that polyp number to stolon 
length was a constant in a clone ofP. carnea. Buss (2001) established the 
relationship of spacing between a polyp and bud along a stolon of a given 
radius to polyp size in experimental colonies. This relationship was 
compared to a derived relationship that balanced opposing forces of 
the pressure of fluid over a cross-sectional area imparted by a polyp with 
the shear forces of gastrovascular fluid moving past itself and along the 
lumen wall of a stolon. The predicted relationship between polyp size 
and spacing between polyps normalized to stolon radius was observed 
and, therefore, is consistent with the hypothesis that shear stress signals 
polyp morphogenesis. 

Dudgeon and Buss established experimental replicates consisting of 
single polyps of P. camea in the middle of a linear stolon of defined length 
(Fig. 3.3A) and grew them in six different seawater treatments in which 
viscosity was manipulated by the addition of Dextran over the range of 
1.0 cp (unmanipulated control seawater) to 2.35 cp. Shear stress equals the 
product of dynamic viscosity and the velocity gradient across the stolon 
radius, so changing viscosity indirectly manipulates shear stress during 
gastrovascular transport. Growth rates of polyps and stolons were assayed 
in each treatment after 3 days and significant differences in morphological 
patterning were evident (Fig. 3.3B and C). An increase in seawater 
viscosity increased the rate of production of polyps per unit stolon length 
relative to controls (Fig. 3.4). In contrast, no relationship was evident 
between seawater viscosity and the number of new stolon tips (Fig. 3.5). 
These results suggest that shear forces along the lumen wall of stolons are 
signals transduced by endodermal cells that regulate polyp morphogene¬ 
sis. Shear stress appears not to be involved in signalling stolon branch 
formation. Buss (2001) suggested that circumferential wall stress may act 
as a mechanosensor signalling stolon branch formation. 

The limits of morphological plasticity are potentially broad such that 
runner phenotypes may be transformed into sheet phenotypes and vice 
versa. Developmental plasticity of colony morphology is regulated by 
multiple signalling mechanisms including ROS and nutrient concentra¬ 
tions associated with gastrovascular transport and the hydrodynamic 
forces of that transport. Physiological regulation of development by mul¬ 
tiple pathways also occurs in well-studied vascular systems of vertebrates 
and there are strong similarities between them. Both hydrozoan and 
vertebrate vascular systems are regulated by ROS and hydro-(hemo-) 
dynamic factors suggesting common mechanisms underlying plasticity. 
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Figure 3*3 Single polyp/single stolon systems of P. carnea used in viscosity 
manipulation experiments. (A) A representative replicate prior to manipulation— 
polyp positioned near the centre of a 2.7 mm length of stolon. (B) A control replicate 
after 72 h, note single new polyp bud emerging at lower right. (C) A replicate in 1.93- 
fold viscosity of seawater following 72 h treatment, note three new polyps formed close 
together to the left of the original polyp. 


2.4. Vascular architecture and constraints on plasticity 

The range of morphological plasticity of a colony may exceed the differ¬ 
ences among genotypes, or even among species, but that range itself is 
likely to vary among colonies. For instance, in H. symbiolongicarpus, rep¬ 
licate colonies of a runner genotype were transformed more easily into a 
sheet phenotype than replicates of a sheet genotype were transformed into 
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Figure 3*4 Production of new polyps relative to stolon growth in colonies grown 
under different seawater viscosities. Values represent mean =b s.e. of N = 5 replicates 
per treatment. The allocation to polyp growth increases relative to stolon growth 
with increased seawater viscosity up to 2.2 cp. 
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Figure 3-5 Production of new branch tips of stolons grown under different seawater 
viscosities. Values represent mean d= s.e. of N = 5 replicates per treatment. New 
stolon tips branch independently of the dynamic viscosity of seawater. 
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a runner phenotype (Dudgeon and Buss, 1996). There are two possible 
explanations for this asymmetry. Genetic differences between the clones 
used may have differed in their respective reaction norms to surgical 
perturbations of stolons. Alternatively, differences in their existing gas- 
trovascular architectures may have influenced their respective ranges of 
plasticity. The greater extent of laterally branching anastomosed stolons 
of sheet colonies should facilitate sharing of resources and coordination of 
behaviours within the colony (i.e. enhanced colony integration). At the 
low Reynold’s numbers that occur in the gastrovascular system ‘slug flow’ 
characterizes transport of fluid, so ‘mixing’ of nutrient-laden fluid with 
nutrient-poor fluid occurs only by insertion of these respective parcels at 
intersections of stolon branches (M. LaBarbera, personal communica¬ 
tion) . Enhanced colony integration may reduce disparities in local envi¬ 
ronmental signals at different points in the colony and buffer against 
different morphogenetic responses. 

We tested the hypothesis that a runner morphology displays greater 
plasticity in (1) the polyp/stolon investment and (2) growth among sto¬ 
lons within a colony in response to food limitation than does a sheet 
morphology. A single ‘sheet morph’ genotype of P. carnea was used to 
unconfound genotypic and morphological sources of plasticity by surgi¬ 
cally constructing experimental morphologies. Plasticities of different 
experimental morphologies were assessed by comparing patterns of 
growth of colonies before and after a transition from a ‘fed to repletion’ 
environment to a ‘food-limiting’ environment. In addition to the sheet 
morph treatment (unmanipulated controls), two other morphology treat¬ 
ments were established. The experimental treatment consisted of remov¬ 
ing lateral stolons that connected radially oriented stolons. A ‘sham- 
manipulated’ sheet served as a control for severing stolon connections, 
but lateral stolons were not removed enabling stolons to fuse. Colonies 
remained in a food replete environment for 4 days after initiating treat¬ 
ments and then were switched to a food-limiting environment in which 
only the same selected polyps received a single brine shrimp nauplius larva 
on feeding days for the remainder of the experiment. 

In the food replete environment, all replicates showed similar patterns of 
gastrovascular transport (Fig. 3.6) and growth among stolons within a 
colony (Fig. 3.7). In the food-limiting environment, stolons in sheet col¬ 
onies continued to grow and produce polyps at similar rates, but runner 
colonies did not. In these colonies, stolons with unfed polyps grew more 
slowly (Fig. 3.7), but increased production of polyps (relative to stolons) 
compared to stolons with fed polyps (Fig. 3.8). Moreover, differences in 
gastrovascular behaviour were consistent with those observed with respect 
to polyp and stolon growth. In sheet colonies, both stolons with fed polyps 
and those with starved polyps showed similar maximum relative diameters 
when oscillating (Fig. 3.9). In contrast, stolons of runner colonies with 
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Figure 3.6 Relative diameters (diameter of the stolon lumen divided by the total 
[external] diameter of the stolon) near stolon tips during a period of maximal 
gastrovascular transport after feeding. Measures were taken prior to beginning the 
experiment. Data points in the time series represent the mean (dzs.e.) of four stolon tips 
in each of five colonies for each treatment (N = 20). Gastrovascular behaviour shows 
no differences among morphology treatments prior to limiting colony food supply. 


starved polyps exhibited lower maximum relative diameters, hence, dimin¬ 
ished gastrovascular transport compared to stolons in the same colony with 
fed polyps (see also Bumann and Buss (2008) for the distribution of labelled 
food in colonies). These results indicate that colonies that vary in form, yet 
derive from the same clone, show different morphological plasticities in 
response to food supply. Moreover, the differences in plasticity resulted 
from differences in the transport of gastrovascular fluid imposed by the 
architecture of the gastrovascular system. 

Blackstone (2001) and Bumann and Buss (2008) have also conducted 
experiments with different objectives using P. carnea in which only 
selected polyps were fed with slightly different outcomes than reported 
above. Both prior studies observed that fed areas of colonies grew polyps 
and stolon tips more rapidly than unfed areas of colonies that continued 
linear extensions of stolons. Selected polyps in the study by Blackstone 
were fed to repletion but colonies in some treatments in the study by 
Bumann and Buss (2008) were food-limited. Like the latter study, in our 
study the selectively fed polyps were food-limited and some surgically 
manipulated to an experimental morphology. In all three studies, 
environmentally induced plasticity established patterns of vascular archi¬ 
tecture and subsequent expressions of plasticity. 

Vascular architecture is important in colony development in at least 
two ways. Vascular architecture defines the characteristics of gastrovascu¬ 
lar flow, including hydrodynamic forces on endodermal cells lining the 
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Days after initial surgery 

Figure 3*7 Specific growth rate of stolons having fed (open circles) or unfed (solid 
circles) polyps in (a) sheet colonies, (b) sham-manipulated colonies and (c) runner 
colonies over a 19-day time course after surgically manipulating colony form. 
Arrows on day 4 indicate the transition to food limitation. Data points represent 
means (zhs.e.) of N > 3 stolons of each feeding level per colony and M = 5 colonies 
per treatment. On day 19, data points that have the same letters are not significantly 
different. No differences in growth rate between stolons with fed versus unfed polyps 
among morphology treatments prior to onset of food limitation (arrow). After 15 
days of food limitation, growth differs between stolons with fed versus unfed polyps 
only on runner colonies. 
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Figure 3*8 Ratio of polyp to stolon growth in runner, sham-manipulated control 
and sheet colony treatments from days 2—19 of the experiment. Dark grey bars 
represent data for fed stolons and light grey bars represent data for unfed stolons. 
Columns represent means (zbs.e.) of N > 3 stolons of each feeding level per colony 
and M = 5 colonies in the sheet and runner treatments, but M = 3 in the sham- 
manipulated control treatment. Columns sharing the same letters are not 
significantly different. Polyp production differed between stolons with fed versus 
unfed polyps on runner colonies, but not on other morphs. 


vascular lumen. Vascular architecture indirectly determines the rate of 
delivery of morphogenetic substances distributed to endodermal cells of 
the lumen wall. Architectural effects are known to affect floral develop¬ 
ment in the modular, vascular plants (Diggle, 1995, 2003). We continue 
to investigate the relative strengths of architectural, genetic and environ¬ 
mental factors in morphogenesis of colonial invertebrates. 

Waddington (1942a) coined the term ‘epigenetics’ to describe inher¬ 
ited changes in phenotype caused by mechanisms other than changes in 
the underlying DNA sequence. This definition remains the broadest 
conceptualization, but typically epigenetics implies DNA methylation, 
histone modification, RNA actions or other regulatory mechanisms that 
transmit a non-genetic cellular memory to daughter cells, and in some 
cases transgenerationally, to specify a phenotype (Bonasio et al, 2010; 
Riddinhough and Zahn, 2010). Those are specific epigenetic mechan¬ 
isms, but the essence of the basic definition includes mechanisms at other 
levels of development. As colonial invertebrates continue to develop as 
they grow, opportunities for variation to be transmitted to differentiating 
cells also occur at larger scales. Existing structures influence the develop¬ 
ment of new structures. The structural state of vascular architecture acts as 
an epigenetic factor influencing the range of morphological plasticity of 
colonies. Lateral and circumferential vascular connections buffer against 
disparities in gastrovascular transport in different parts of a colony, damp¬ 
ening morphological plasticity of colonies. In response to changes in food 






Figure 3*9 Relative diameters (lumen diameter of the stolon divided by the total 
[external] diameter of the stolon) near tips of stolons having fed (open circles) or 
unfed (solid circles) polyps in (A) sheet colonies, (B) sham-manipulated colonies and 
(C) runner colonies, during a period of maximal gastrovascular transport after 
feeding. Measures were taken on day 19 of the experiment. Data points in the 
time series represent the mean (dzs.e.) of N > 3 stolon tips of each feeding level in 
each of M = 5 colonies per treatment. Treatments sharing the same letters indicate 
that their relative maximum diameters are not significantly different. Gastrovascular 
behaviour in stolons with unfed polyps differed from those with fed polyps, but only 
among experimental runner morphs. 
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supply, the reaction norm of runner morphs is broader than that of sheet 
morphs, and the epigenetic factor controlling the norm of reaction is the 
dynamic pattern of vascular connectivity in colonies. 


3. Colony Development as a Dynamical System 

Buss and Blackstone (1991) applied Waddington’s (1942b) view of 
development as a dynamical system to colonial organisms. Their experi¬ 
ments, with a single clone of Hydractinia echinata, provided evidence that 
colonies regulate rates of polyp and stolon production, but morphology was 
not canalized by selection. They suggested that correlations between form 
and function were largely the outcome of a self-organizing process 
(se«i)/Bonabeau et al, 1997) of dynamical behaviours of a developing col¬ 
ony. Genetic programmes of colonial cnidarians do not specify a stereo¬ 
typical morphology, nevertheless repeatable and identifiable morphs are 
formed. Morphogenetic patterning occurs in response to physiological 
signals established by the interactions between collective polyp behaviour 
and the existing vascular architecture (Blackstone and Buss, 1992; Dudgeon 
and Buss, 1996; Wagner et al, 1998; Blackstone, 1999, 2001, 2003; 
Dudgeon etal, 1999, Bumann and Buss 2008; data presented above). The 
collective behaviour of a colony emerges from simple behaviours of polyps 
interacting with one another and their environment (Wagner etal, 1998; 
Dudgeon etal, 1999; Dudgeon, 2001a,b). With each new polyp or stolon 
branch added to the existing vasculature the physiological state may change. 

Colonies develop, therefore, by genotype- and context-dependent spa- 
tiotemporal patterns of gene expression dynamically responding to the 
physiological state (environment) and architecture of the gastrovascular 
system (Fig. 3.10). A population of cnidarian (or other invertebrate) col¬ 
onies may include a variety of genotypic classes that vary in phenotypic 
expression in response to ontogenetic stage and environment. Some gen¬ 
otypes may be canalized early, some late, some not at all. There may be 
latitude in phenotypic values early in ontogeny that is canalized to a narrow 
outcome later in ontogeny (Fig. 3.10A), but Buss and Blackstone did not 
observe such a pattern in the genotype they studied. Plastic responses of 
colonies may vary from linear to threshold responses. Fig. 3.10B, C and E 
represent genotypes in which plasticity declines during ontogeny. 
Genotypes represented by Fig. 3.10B show a non-linear, but monotonic 
change in phenotype across environments early in ontogeny and a less 
responsive, linear change of phenotype across environments late in ontog¬ 
eny. Genotypes represented by Fig. 3.10C show a strong threshold response 
of phenotype to environmental change early, but a less plastic, linear 
response late in ontogeny. Possibly, genotypes may have threshold plastic 
responses to environmental change throughout ontogeny (Fig. 3.10D). If 
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A Canalization of phenotype 



B Monotonlcally plastic phenotype 




D Threshold plasticity throughout 



E Multiple stable points late in ontogeny 



Environment 


Figure 3* 10 Phenotypic outcomes for different genotypic classes with respect to 
environmental conditions and ontogenetic stage (i.e. architectural inertia). Surfaces 
represent genotypes and show a range of hypothetical interactive effects with 
environment and existing architecture on phenotypes. (A) Canalized phenotype, (B) 
monotonically plastic phenotype with greater sensitivity to environmental variation at 
early stages, (C) threshold response of plasticity in early, but not late ontogeny, (D) 
threshold response of plasticity throughout ontogeny, (E) multiple stable phenotypic 
points late in ontogeny, but a threshold response early. These five patterns represent a 
sample of genotypic classes of effect that may exist in populations. 


effects of vascular architecture increase later in ontogeny, then a plausible 
phenotypic response late in ontogeny may be one of discrete, alternative 
stable phenotypes of a genotype over some range of environmental con¬ 
ditions (Fig. 3.10E, see also Zeeman, 1976). Colonies may be very plastic 
early in ontogeny, but late in ontogeny they are not. Instead, colonies 
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become canalized to one of a few discrete, stable outcomes. Which of the 
multiple possible phenotypes prevails is contingent upon interactions 
between the existing phenotype and environment as the colony grows. 

The existing vascular architecture for a given ontogenetic stage is as 
important as the environmental conditions in determining plasticity, 
because architecture creates conditions for self-maintenance. We predict 
that canalization of a phenotype is much greater at later ontogenetic stages 
given the increasing inertia of architectural effects at large colony sizes. 

Modular organisms offer an illuminating contrast to unitary model 
systems (e.g. fruit flies, worms, mice) for study of development 
(Blackstone and Bridge, 2005). Both show tightly regulated spatiotem- 
poral patterns of gene expression in the formation of their respective 
structures or modules. Body plans of modular organisms may vary greatly 
in the number, spacing and arrangement of modules relative to one 
another. At the level of patterning modules, environmental signalling 
plays an expanded role in morphogenesis (Blackstone and Bridge, 2005). 

Developmental genetic programmes of unitary and modular organ¬ 
isms can be arrayed along a spectrum of variation in the extent to which 
environmental signalling influences morphological pattern. At one 
extreme, gene expression in unitary organisms in response to the envi¬ 
ronment is largely the internal environment created by prior gene expres¬ 
sion in the target and neighbouring cells with less opportunity for external 
environmental stimuli to influence morphological patterning. 
Phenotypic effects are more constrained to variation in size or shape of 
parts, with only rare variation in the number of parts, or body axes 
(Blackstone and Bridge, 2005). Toward the other extreme, gene expres¬ 
sion in clonal organisms is highly responsive to external environmental 
signals. One can interpret the environmental signalling axis as represent¬ 
ing the extent to which genes and their products directly drive morpho¬ 
logical development. With increased environmental signalling, the role of 
genes changes to one of establishing a physiological state that is thereafter 
maintained by epigenetic feedbacks (or ‘local rules’; Ulam, 1962; 
Braverman and Schrandt, 1966; Buss, 2001; Kaandorp and Kiibler, 
2001) and modified in response to environmental signalling. 

Such dynamic regulation of development and resulting adaptive plas¬ 
ticity may enable maintenance of greater genetic variation in natural popu¬ 
lations of clonal organisms. In Hydractinia spp., distributions of morpholog¬ 
ical phenotypes in natural populations remain variable in the face of 
intraspecific competition in which a clear dominance hierarchy of runners 
to sheets is known (Buss and Grosberg, 1990; Yund, 1991). A number of 
other processes, especially environmental variability, may act to maintain 
this variation through their effect on the dynamics of development. 
Selection against competitively inferior genotypes may be weakened, or 
masked, if competitive ability co-varies with morphological plasticity. 
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Future Directions and Conclusions 


Much has been learned in the past two decades about the physiological 
regulation of developmental plasticity in hydrozoans, which may serve as a 
model for many clonal organisms. These studies have reinvigorated old 
questions and stimulated new ones. We identify six avenues of research to 
enhance understanding of relationships between the development and 
ecology of colonial cnidarians, specifically, and clonal organisms, generally. 


4.1, The heritability of colony form 

Evolutionary ecology studies of colonial cnidarians clearly documented a 
genetic basis of colony form (Foster, 1979; McFadden et al. 1984; Buss 
and Grosberg, 1990; Blackstone and Buss, 1991; Yund, 1991; Knowlton 
etal., 1992; Bruno and Edmunds, 1997). In light of data emphasizing 
morphological plasticity, disparities among estimates from different stud¬ 
ies of broad-sense heritability and the lack of sufficient data to estimate 
narrow-sense heritability, this question needs greater resolution. It would 
be useful to have good estimates of the components of phenotypic 
variance for the different morphometric characters that capture the qual¬ 
itative descriptions of colony form, in the context of competitive relation¬ 
ships among morphs. Estimates of additive genetic variance (Fa) in 
populations could test the hypothesis that selection from competition 
strongly influences morphological phenotypes. Narrow-sense heritability 
and the response to selection from competition can be measured only 
from estimates of F A . Additional estimates of contributions of non¬ 
additive genetic, and environmental, variances to phenotypic variance, 
from breeding studies, would resolve disparities among studies estimating 
broad-sense heritability and the scope for plasticity. 


4.2. The heritability of morphological plasticity 

We ask whether morphological plasticity observed among colonial inver¬ 
tebrates is, itself, heritable. In other words, is there positive covariation 
between the plastic responses of parents and their offspring? Genotype by 
environment interactions (G*E interactions) contribute a significant frac¬ 
tion to phenotypic variance within populations of colonial cnidarians 
(Bruno and Edmunds, 1997), but it has yet to be quantified relative to 
other components of variance. The extent to which G*E interactions 
represent heritable variation would contribute to understanding quanti¬ 
tative genetics, as well as the developmental plasticity of clonal organisms. 
Whether plasticity is a heritable trait was hotly debated by theorists during 
the 1990s (Via and Lande, 1985; Schlichting, 1986; Stearns, 1989; 
Schemer, 1993b; Schlichting and Pigliucci, 1993; Via, 1993; Via etal., 
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1995; Pigliucci 1996) and colonial invertebrates appear to be ideal systems 
for empirical studies addressing the question because clonal variation in 
plasticity is known (Bruno and Edmunds, 1997) and cnidarians are ame¬ 
nable to the requisite experiments. 

4.3. Architectural versus genotypic effects 

A related issue is the relative contribution of architectural versus geno¬ 
typic variation in setting the limits of phenotypic plasticity. Experiments 
using a single clone suggested that vascular architecture influenced plastic 
responses of phenotypes. Experiments that cross different vascular geom¬ 
etries with multiple genotypes are needed to assess the importance of 
architectural effects more generally. The likely existence of genetic var¬ 
iation for phenotypic plasticity (whether or not it is heritable) raises two 
further questions. Are effects of vascular architecture consistent across 
genotypes? What are the relative effect sizes of plasticity associated with 
different architectures compared to those among genotypes? 

Variation in competitive ability among morphotypes could be studied in 
the context of architectural versus genotypic effects. Competitive ability is 
linked to morphology, but it is unknown if competitive ability changes as 
the morphology of a colony is altered. If competitively superior morphs 
become inferior following a change in morphology, then competitive 
ability appears more tied to form than genotype. Alternatively, a genotype 
may be competitively superior regardless of its morphological plasticity. 

4.4. Physiological mechanisms underlying morphological 
plasticity 

Blackstone demonstrated that ROS-signalling is a mechanism by which 
environments are sensed through polyp behaviour triggering adaptive 
morphogenetic plasticity. Redundancy of signalling pathways for plastic¬ 
ity would be advantageous in complex dynamic environments. Where 
mitochondria perceive anoxia, the presence of oxygen in the cell should 
prevent signalling of hypoxia-induced compensatory pathways (genes 
encoding hypoxia-inducing factor 1 [ HIF-1] have been identified 
in Hydra magnipapillata; D. Bridge, personal communication). Under 
natural conditions redox state and p0 2 are correlated (Lesser, 2006). 
It is unknown whether colonies in nature operate multiple pathways of 
oxygen-dependent regulation. Hydractiniid hydrozoans, like several 
Hydractinia species and P. carnea, live on shells occupied by hermit crabs 
in sandy to muddy marine habitats characterized by variable p0 2 , redox 
potential in sediments and temperature (Fenchel and Reidl, 1970). The 
behaviours of their host hermit crab subject colonies to fluctuations in 
p0 2 ranging from hyperoxia to extreme hypoxia within minutes as the 
shell inhabited by a colony is buried in reduced sediments (Vasquez, 
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2009). In addition to ROS signalling, it is plausible that colonies operate 
other signalling pathways to sense environmental p0 2 directly. 

Morphological plasticity of colonies in response to manipulated seawa¬ 
ter viscosity suggested shear stress signalling in colony-level morphogene¬ 
sis. Additional experiments are needed to test that possibility. Previous 
experiments manipulating the dynamic viscosity of seawater used minimal 
architectures (single polyp, single stolon). Shear stress could have increased, 
decreased or remained the same depending on possible changes of transport 
rates in the gastrovascular system in response to viscosity manipulations. 
Controlled viscosity manipulations with concurrent measures of transport 
rates could determine the effect of manipulation on vascular shear stress and 
consequent morphogenetic effects. Simple observations of correlations 
between flow rate and the fraction of cells undergoing mitosis (i.e. mitotic 
index) in unmanipulated colonies suggests a direct link between gastrovas¬ 
cular flow rate and stolon elongation (Fig. 3.11). 


4.5. Empirical tests of nutrient-driven polyp behaviour 

The model presented by Wagner et al. (1998) was based on behaviour of a 
single polyp fed a defined quantity of food. Manipulating nutrient content 
could test predictions of this model. Experimental diets can be given to 
colonies to test specific levels of nutrients that elicit polyp oscillatory 
behaviour (Dudgeon etal., 2009). If model predictions are consistent with 
observed polyp behaviours, then scaling up the model to two- and three- 
polyp cases arranged in various geometries would be possible. 

4.6. Colonial hydrozoans as models of vascular development 

Colonial hydrozoans grow indeterminately, lack a predefined size or 
shape and continuously develop as they remodel vascular architecture 
to accommodate changes in state associated with growth, competitive 
encounters, substrate availability and environmental changes associated 
with food supply or other factors. Likewise, vascular systems of vertebrates 
grow indeterminately, lack a predefined size or shape and continuously 
develop as they remodel vessel architecture to accommodate internal 
changes in state associated with growth, maturation, reproduction, 
wound healing and various pathologies, despite the relative constancy 
of the organism’s external morphology. These parallels make hydrozoans 
useful models for in vivo study of vascular function. Comparisons 
between organ systems and colonial organisms, may give insight into 
common physiological mechanisms to regulate adaptive plasticity of vas¬ 
cular systems. Buss (2001) suggested that genes homologous to vertebrate 
vascular endothelial growth factors (VEGF), known to act in vascular 
proliferation, be sought in cnidarians. Since then, two VEGF genes and 
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Figure 3 -11 Positive covariation between rates of gastrovascular transport and mitosis 
of stolon cells. (A) Normalized flow rate in two colonies (open or solid circles, 
respectively) was calculated as the maximum amplitude of lumen oscillations 
normalized to external stolon diameter, divided by the period (in seconds) of the 
oscillation cycle. (B) Mitotic index is calculated as the fraction of living cells (indicated 
by staining with DAPI) that are undergoing mitosis (indicated using the commercially 
available BrdU/anti-BrdU staining immunoassay). (C) Composite fluorescent image of 
a stolon showing staining by DAPI (blue; to estimate cell density in a region of stolon) 
and/or FITC (green; to identify mitotic cells). Image constructed from a Z-series stack of 
images from bottom to the top of the stolon lumen at 400 X magnification. Mitotic rates 
were estimated in two colonies (distinguished by solid vs. open symbols) along a linear 
length of stolon proximally near the polyp (circles) and distally near the stolon tip 
(triangles). Symbols and error bars represent means zb 95% confidence intervals of 
n = 4 colonies at each sampling period. (See Plate no.5 in the Color Plate Section.) 







138 


Dudgeon and Kiibler 



Figure 3.11 continued 

a VEGF receptor have been identified in P. carnea (Seipel etal., 2004; D. 
Bridge, personal communication). Expression in stolons has not been 
investigated, but P. carnea VEGF and VEGFR are expressed in the endo- 
dermal tissue of tentacles and of the radial canals of the medusa gastrovas- 
cular system, consistent with VEGF signalling for morphogenetic processes 
leading to tube formation (Seipel etal, 2004). Most interestingly, factors 
known to regulate vertebrate VEGF signalling, redox state, p0 2 and vas¬ 
cular shear stress, are the same as those implicated in the adaptive plasticity 
of hydrozoan colonies living in patchy food, and hypoxic soft sediment, 
environments (Chen etal, 1999; Osawa etal., 2002; Shay-Salit etal, 2002; 
Urbich etal. 2002; Holmes and Zachary, 2005; Milkiewicz etal, 2007). 

4.7. Conclusions 

Physiological and anatomical (architectural) mechanisms regulate adaptive 
plasticity during development of hydrozoan colonies. The potential for 
rapid, large scale environmental changes call for greater understanding of 
physiological processes in ecology. Clonal organisms are particularly good 
models for studies of plasticity in response to environmental change because 
their development and ecological performance are tightly coupled. As 
asserted by both Bonner (1965) and Buss and Dick (1992), unravelling the 
physiological mechanisms that lie on the map between genotype and phe¬ 
notype, in the ‘middle ground’ of biology, still holds the key to understand¬ 
ing the complexity of phenotypic variation in a changing environment. 
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Abstract 

There is general consensus that solar ultraviolet radiation (UVR) negatively 
impacts many marine species. Echinoderms are ubiquitous within the 
marine environment, with members of the phyla often long-lived and 
numerically dominant within the benthic macrofauna, consequently the 
impact of UVR on the population dynamics of these organisms will influence 
marine communities and ecosystems. Research to date has shown that 
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exposure of echinoderms to solar UVR can, affect reproduction and devel¬ 
opment, change behaviour, cause numerous biochemical and physiological 
changes and potentially cause increased mutation rates, by causing DNA 
damage. There is also considerable evidence that echinoderms utilise sev¬ 
eral different mechanisms to protect themselves against excessive UVR and 
subsequent UVR-induced damage. However, these protective mechanisms 
may pose conflicting selection pressures on echinoderms, as UVR is an 
additional stressor in oceans subjected to anthropogenic-induced climate 
change. This review summarises our knowledge of the effects of UVR on the 
Echinodermata. We outline the research conducted to date, highlight key 
studies on UVR that have utilised echinoderms and look to the future of UVR 
research in a rapidly changing ocean. 



Introduction 


1.1. Ultraviolet radiation and the marine environment 


Ultraviolet radiation (UVR) enters the marine environment in the form 
of UVB (290-315 nm) and UVA (315-400 nm). While UVR is effi¬ 
ciently attenuated through the water column, it can penetrate to depths 
of greater than 16 m (10% of surface UVB) and 46 m (10% of surface 
UVA) (reviewed by Tedetti and Sempere, 2006), where it plays a pivotal 
role in shaping the lives of marine organisms. The biological effects of 
being exposed to UVR have been extensively studied and comprehensive 
reviews have outlined the impacts of exposure on marine organisms at the 
physiological and molecular levels (Dahms and Lee, 2010) and on marine 
ecosystems (De Mora etal., 2000; Hiider et al., 2007). This review sum¬ 
marises our knowledge of the effects of UVR on the Echinodermata. We 
outline research to date, highlighting key studies on UVR that have 
utilised echinoderms and look to the future of UVR research in a rapidly 
changing ocean. 


1.2. Phylum Echinodermata 

Echinoderms (Phylum Echinodermata) are a diverse, exclusively marine 
group of invertebrates that consists of over 13,000 extinct species (15 classes) 
and 7,000 extant described species within five classes (Asteroidea, 
Echinoidea, Holothuroidea, Ophiuroidea and Crinoidea). Echinoderms 
are ubiquitous within the marine environment, and members of the phyla 
can often be long-lived, a numerically dominant component of the benthic 
macrofauna, and a keystone species within regions they occupy. 
Consequently, processes that impact on the population dynamics of these 
organisms can, and do, influence marine communities and ecosystems. 
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Echinoderms have a complex life cycle, alternating between the long-lived 
adult stage and the short-lived embryonic and larval stages. Within this 
general reproductive pattern is a range of life-history strategies, including 
brooding and vivipary. Echinoderm gametes are readily obtained and have 
played a pivotal role in enhancing our understanding of reproductive pro¬ 
cesses such as fertilisation and embryonic development, including a signif¬ 
icant contribution to our understanding of the effects of UVR on 
reproduction. 

1.3. The history of UVR research in echinoderms 

Echinoderms have played a pivotal role in our understanding of the role of 
UVR on marine organisms. Research concerning the effects of UVR on 
echinoderms dates back to the early 1900s, with the total number of 
publications reaching ~65 by 2009.The research can be characterised 
by two main phases: (1) a period in the 1950s and 1960s involving 
laboratory studies on the role of UVR in developmental biology and 
(2) during the late 1980s to the present when an appreciation of the 
ecological importance of UVR and ozone depletion drove an increased 
number of physiological and field studies (Fig. 4.1). Of the published 



Figure 4.1 Meta-analysis ((Web of Science (Thomson Reuters, NY) database 
search: key words; Echinoderm* OR sea urchin OR starfish OR sea star OR 
echinoid OR sand dollar OR sea cucumber OR brittlestar OR echinoid OR 
asteroid OR holothuroid OR ophiuroid OR crinoid OR featherstar AND 
ultraviolet) of the number of published research papers on UVR in echinoderms. 
Fig. 4.1 represents the cumulative number of papers from and including the year 
1900 through to 2009. 
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studies, approximately 80% are on echinoids, with remaining classes the 
focus of 14 peer-reviewed studies. Furthermore, within this data set, 75% 
examine UVR responses in the embryological or larval stages, with lesser 
attention given to the post-settlement life-history stages. 



Responses in echinoderms to UVR exposure 


The outcome of exposure to UVR for echinoderms can be in the 
form of avoidance, damage and repair (summarised in Fig. 4.2). These 
outcomes may occur at a number of levels, from the molecular, physio¬ 
logical, developmental and behavioural and ultimately through to the 
population and ecosystem levels. Based on this classification of responses, 
previous research which examines the effects of UVR in echinoderms is 
catalogued in Table 4.1, with responses grouped according to echinoderm 
class. In the following sections, we review the research that has contrib¬ 
uted to an understanding of these responses to UVR in the phylum. 


2.1. Avoidance strategies 

2.1.1. Cellular 

2.1.1.1. Mycosporme-like amino acids Sunscreening is an important 
UVR mitigating strategy, and echinoderms contain a range of 
molecules capable of absorbing solar radiation and acting as 
photoprotectants (Table 4.11). Of these the mycosporine-like amino 
acids (MAAs) are the best studied, and most important group of 
photoprotective pigments in echinoderms (see reviews by Shick and 
Dunlap, 2002; Rastogi et al., 2010). Their occurrence has been 
recorded in all echinoderm classes, and from polar to tropical species 
(Table 4.2). Within the phyla, eight different MAAs have been 
reported, with those found in the highest concentrations being 
mycosporine glycine, Porphyra-334 and shinorine. Concentrations of 
MAAs vary depending on species and tissue, as well as among 
developmental stages. For example, in tropical sea cucumbers, MAAs 
are concentrated in the epidermis and ovaries rather than in the viscera, 
consistent with the allocation of MAAs to tissues that will experience 
greater UVR exposure (Shick etal., 1992; Bandaranayake and Rocher, 
1999). Similarly, Karentz etal. (1997) and McClintockandKarentz (1997) 
noted that the greatest MAA concentrations were found in the ovaries of 
the sea urchin Sterechinus neumayeri, with lower concentrations found in 
the body wall presumably reflecting the greater exposure of embryos and 
free-swimming larvae to UVR at the top of the water column compared 
to the adults found at a greater depth. 
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Figure 4« 2 The major responses to UVR shown by echinoderms. Avoidance and 
protection (at the behavioural and cellular level) are primary mechanisms of defence 
and provide organisms with the ultimate safeguard against UVR. When defence is not 
possible, impairment occurs in the major forms of lipid, protein and DNAdamage. 
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Echinoderms obtain MAAs from their diet (Rastogi etal. , 2010), and a 
number of experiments have indicated the importance of diet on the 
concentration and composition of MAAs within this group. Carroll and 
Shick (1996), by altering the algal diet of the sea urchin Strongylocentrotus 
droebachiensis, were able to produce eggs that varied 25-fold in their 
MAA concentrations and also their composition. Similar findings were 
shown in later work on S. droebachiensis (Adams and Shick, 1996, 2001) 
and for Evechinuschloroticus (Lamare etal., 2004). The interconversion of 
MAAs has been reported in the sea cucumber Thelenotaananas (Shick and 
Dunlap, 2002), and has been proposed to be the result of the action of 
bacteria within the gut tract hydrolysing shinorine and porphyra-334 to 
mycosporine glycine. The activity of the bacteria mirrored changes in 
MAA composition along the gut tract of the animal. 

MAA concentrations may vary depending on ambient UVR irradi- 
ance. Karentz et al. (1997) found significantly higher concentrations of 
MAAs in the tissues of S. neumayeri individuals collected from intertidal 
and shallow subtidal depths than in deeper individuals. Lamare etal. (2004) 
noted that the concentrations of MAAs in the eggs of E. chloroticus varied 
six-fold along a 40 km fiord gradient, and were correlated with ambient 
levels of UVB irradiance found along the fiord. These results likely reflect 
a passive response to changes in MAA concentrations in food sources of 
the sea urchins and not an active increase in the uptake of MAAs by sea 
urchins in response to higher irradiances. This conclusion is supported by 
the observation that, in the laboratory, S. droebachiensis does not accumu¬ 
late MAAs in response to UVR exposure (Adams etal., 2001). 

Embryonic and larval stages of echinoderms also contain MAAs 
(Table 4.2), which are endowed to the eggs maternally (Carroll and 
Shick, 1996). In the laboratory, MAA concentrations decrease during 
initial development (i.e. pre-feeding stages) with no evidence that they 
are supplemented in the larvae by the absorption of MAAs from seawater 
(Adams and Shick, 1996, 2001). Whether decreases in MAAs in labora¬ 
tory-reared embryos reflect patterns in situ is unknown, however, it is 
conceivable that the dietary accumulation of MAAs by feeding larvae 
occurs in the water column, with their phytoplankton and bacterial food 
likely to contain MAAs (Sinha etal., 1998). 


Figure l\.2 (cont.) Echinoderms are equipped with subsequent repair mechanisms 
including high-cost strategies such as protein replacement and low-cost strategies such 
as photoreactivation. "When repair is not sufficient, damage can result in impaired 
embryonic and larval development, slower growth or eventually mortality of 
individuals. Ultimately, population effects can be seen in reduced recruitment and 
complex life-history trade-offs leading to largely undocumented impacts on marine 
communities in which echinoderms play large roles. 





Table 4.1 Summary table of published studies directly looking at impacts of UVR on echinoderms 




1.1 Behavioural responses 



UVR response 

Life-history 

stage 

Collection site/field Location 

Source 

Echinoidea 

Strongylocentrotus 

Phototaxis 

Adults 

Crow Island, Maine, USA (Lab) 

(Adams, 2001) 

droebachiensis 

Paracentrotus liuidus 

Covering 

Adults 

Lough Hyne Marine Nature 
Reserve, County Cork, 

Ireland (Lab/Field) 

(Crook etal., 1999; 
Crook and Barnes, 
2001; Verling etal., 
2002) 

Tripneustes ventricosus 

Covering 

Adults 

Discovery Bay Marine 

Laboratory, Jamaica (Lab/Field) 

Kehas etal. (2005) 

Strongylocentrotus 

droebachiensis 

Covering 

Adults 

Mingan Islands, northern Gulf 
of St. Lawrence, Canada (Lab) 

Dumont etal. (2007) 

Lytechinus 

variegatus 

Covering 

Adults 

Gulf Specimen Marine 

Laboratories, Panacea, Florida, 
USA (Lab) 

Sigg etal. (2007) 

Asteroidea 

Pisasterochraceus 

Habitat selection 

Adults 

Eastsound, Orcas Island, 

Washington, USA (Lab/Field) 

Burnaford and 

Vasquez (2008) 

Ophiuroidea 

Ophiopholis aculeata 

Habitat selection 

Adults 

Mingan Islands, northern Gulf 
of St. Lawrence, Canada 
(Lab/Field) 

Drolet etal. (2004) 

( continued) 
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1.2 Molecular and cellular responses 



UVR response 

Life-history 

stage 

Collection site/field Location 

Source 

Strongylocentrotus 

purpuratus 

Fertilisation, 

Egg cleavage 

Eggs 

nr (Lab) 

Rustad (1959) 

Anthocidaris crassispina 

Cytokinesis 

Eggs 

Misaki Marine Biological 

Station, Japan (Lab) 

Ikeda (1965) 

Clypeasterjaponicus 

Cytokinesis 

Eggs 

Misaki Marine Biological 

Station, Japan (Lab) 

Ikeda (1965) 

Hemicentrotus 

pulcherrimus 

Cytokinesis 

Eggs 

Misaki Marine Biological 

Station, Japan (Lab) 

Ikeda (1965) 

Paracentrotus lividus 

ATPase activity 

Eggs 

Mediterranean (Lab) 

Petzelt (1974) 

Echinus esculent us 

ATPase activity 

Eggs 

North Sea (Lab) 

Petzelt (1974) 

Sphaerechinusgranularis 

Polyspermy 

Gametes 

Roscoff, France (Lab) 

David etal. (1985) 

Hemicentrotus 

pulcherrimus 

Protein 

complement, 

photoreversal 

Embryos 

nr (Lab) 

Akimoto and Shiroya 
(1987a) 

Strongylocentrotus 

droehachiensis 

Sunscreens 

(MAAs) 

Eggs 

Blue Hill Falls, Maine, USA (Lab) 

Adams and Schick (1996) 

Strongylocentrotus 

droehachiensis 

Sunscreens 

(MAAs) 

Gonads/Gut/ 
Test/Spines 

Blue Hill Falls, Maine, USA (Lab) 

Carroll and Shick (1996) 

Sterechinus neumayeri 

Sunscreens 

(MAAs) 

Gonads/ 

Digestive 

tract/Body 

wall 

Hero Inlet, Arthur Harbor, Palmer 
Station, Antarctica (Field) 

Karentz et al. (1997) 
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Abatus shackletoni 

Sunscreens 

Whole 

McMurdo Sound, 


(MAAs) 

animals 

Antarctica (Field) 

Sterechinus neumayeri 

Sunscreens 

Whole 

McMurdo Sound, 


(MAAs) 

animals 

Antarctica (Field) 

Strongylocentrotus 

Sunscreens 

Eggs/ 

Pemaquid Point, Maine, 

droebachiensis 

(MAAs) 

Embryos/ 

Larvae 

USA (Lab) 

Echinarachnis parma 

Sunscreens 

(MAAs) 

Eggs/ 

Embryos/ 

Larvae 

Frazer Point, Maine, USA (Lab) 

Strongylocentrotus 

Sunscreens 

Ovaries, Eggs 

Crow Island & Pemaquid Point, 

droebachiensis 

(MAAs), 


Maine, USA (Lab) 


embryonic 

development 


Evechinus chloroticiis 

Sunscreens 

(MAAs) 

Gonads/ 

Whole 

animals 

Doubtful - Thompson Sound 
complex, Fiordland, New 
Zealand (Field) 

Anthocidaris crassispina 

Oxidative stress, 
fertilisation, 
sperm damage 

Sperm 

Cheung Chau, Hong Kong (Lab) 

Anthocidaris crassispina 

Mitochondrial 

function, 

membrane 

integrity 

Sperm 

Cheung Chau, Hong Kong (Lab) 

Paracentrotus lividns 

Heat-shock 

proteins 

Embryos 

North-Western Sicily, Italy (Lab) 

Paracentrotus Hindus 

HsP expression, 
apoptosis 

Coelomocytes 

Palermo gulf, Italy (Lab) 

Sphaerechinus 

Structural damage. 

Gametes/ 

Bay of Banyuls-sur-Mer, 

granularis 

chromatin 

damage 

Embryos 

Mediterranean coast, 

France (Lab) 


McClintock and 
Karentz (1997) 
McClintock and Karentz 
(1997) 

Adams and Schick (2001) 
Adams and Schick (2001) 
Adams et al. (2001) 


Lamare et al. (2004) 
Lu and Wu (2005a) 
Lu and Wu (2005b) 


Bonaventura et al. (2006) 

Matranga et al. (2006) 

Nahon et al. (2009) 

I |w 
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1.2 Molecular and cellular responses 



UVR response 

Life-history 

stage 

Collection site/field Location 

Source 

Sphaerechinus granulans 

Sperm damage 

Sperm 

Bay of Banyuls-sur-Mer, 

Pruski etal. (2009) 

Strongylocentrotus 

Oxidative stress, 

Embryos 

Mediterranean Coast, 

France (Lab) 

Isles of Shoals, Gulf of Maine, 

Lesser (2010) 

droebachiensis 
Sterechinus neumayeri 

DNA damage 
Oxidative stress 

Embryos 

USA (Field) 

Cape Armitage, McMurdo 

Lister et al. (2010a) 

Tripneustesgratilla 

Oxidative stress 

Larvae 

Sound, Antarctica (Field) 
Aitutaki Lagoon, Cook 

Lister et al. (2010b) 

Paracentrotus lividns 

Stress Proteins 

Embryos 

Islands (Field) 

Palermo, Sicily, Italy (Lab) 

Russo et al. (2010) 

Asteroidea 

Marthasteriasglacialis 

Polyspermy 

Gametes 

Roscoff, France (Lab) 

David et al. (1985) 

Dipiasterias brucei 

Sunscreens 

Whole 

McMurdo Sound, 

McClintock and 


(MAAs) 

animals 

Antarctica (Field) 

Karentz (1997) 

Odontastermeridonalis 

Sunscreens 

Whole 

McMurdo Sound, 

McClintock and 


(MAAs) 

animals 

Antarctica (Field) 

Karentz (1997) 

Odontastervalidus 

Sunscreens 

Whole 

McMurdo Sound, 

McClintock and 


(MAAs) 

animals 

Antarctica (Field) 

Karentz (1997) 

Perknasterfuscus 

Sunscreens 

Whole 

McMurdo Sound, 

McClintock and 


(MAAs) 

animals 

Antarctica (Field) 

Karentz (1997) 


Miles Lamare et at. 






Holothuroidea 

Cucumariaferrari 

Sunscreens 

Whole 

McMurdo Sound, 

McClintock and 


(MAAs) 

animals 

Antarctica (Field) 

Karentz (1997) 

Holotliuria atra 

Sunscreens 

Epidermal 

Lizard Island group, 

Banbaranayake and 


(MAAs) 

tissue/ 

Ovaries/ 

Viscera/Gut 

GBR, Australia (Lab) 

Rocher (1999) 

Stichopusjaponicus 

Body wall 

Body wall 

Zhangzidao Island, 

Beiwei eta /. (2008) 


damage 

tissues 

China (Lab) 


Crinoidea 

Promachocrinus 

Sunscreens 

Whole 

McMurdo Sound, 

McClintock and 

kergttelensis 

(MAAs) 

animals 

Antarctica (Field) 

Karentz (1997) 


1.3 Development/Growth/Mortality 

UVR response Life-history Collection site/field Location Source 

stage 

Echinoidea 


Strongylocentrotus 

purpuratus 

D evelopmental 
abnormality 

Gametes/ 

Embryos/ 

Larvae 

nr (Lab) 

Giese etai. (1955) 

Hemicentrotus 

pulcherrimus 

Developmental 

abnormality 

during 

gastrulation 

Embryos 

nr (Lab) 

Amemiya etal. (1986) 

Cyp east erjap onicus 

Fertilisation, 
gamete fusion 

Eggs 

nr (Lab) 

Hamaguchi and Hiramoto 
(1986) 

( continued) 
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1.3 Development/Growth/Mortality 



UVR response 

Life-history 

stage 

Collection site/field Location 

Source 

Hemicentrotus 

pulcherrimus 

D evelopmental 
abnormality 

Embryos, 

Larvae 

nr (Lab) 

Akimoto etal. (1983) 

Echinometra viridis 

Mortality 

Adults 

Lee Stocking Island, Bahamas (Lab) 

Reaka-Kudla etal. (1993) 

Lytechinus variegatus 

D evelopmental 
abnormality 

Embryos 

Gulf of Mexico, USA 

Steevens etal. (1999) 

Anthocidaris 

crassispina 

Sperm motility, 
fertilisation 

Sperm 

East Peng Chau, China (Lab) 

Au etal. (2002) 

Strongylocentrotus spp. 

Carotenoids 

Gonads/Eggs 

Clallum Bay and the San Juan 
Islands, Washington State, USA 
(Field) 

Lamare and Hoffman (2004) 

Paracentrotus liuidus 

Developmental 
abnormality, 
stress markers 

Embryos 

North-Western Sicily, Italy (Lab) 

Bonaventura et al. (2005) 

Strongylocentrotus 

intermedins 

Ophiuroidea 

Sunscreening, 

Abnormal 

development 

Larvae 

Vostok Bay, Sea of Japan, Russia 

Kipryushina etal. (2010) 

Ophionereis reticulata 

Mortality 

Adults 

Lee Stocking Island, 

Bahamas (Lab) 

Reaka-Kudla et al. (1993) 

Ophiocotna echinata 

Mortality 

Adults 

Lee Stocking Island, 

Bahamas (Lab) 

Reaka-Kudla et al. (1993) 

Ophioderma 

brevispinnm 

Mortality 

Adults 

North Carolina coast, 

USA (Field) 

Johnsen and Kier (1998) 


Miles Lamare et al. 






1.4 DNA damage/repair 




UVR response 

Life-history stage 

Collection site/field 

Location 

Source 

Echinoidea 

Echinarachnius parma 

Developmental delays, 
Photoreactivation 

Gametes/ 

Embryos 

nr (Lab) 

Cook and Rieck (1962) 

Paracentrotus lividus 

Photosensitization 

Sperm 

nr (Lab) 

Colombo (1967) 

Echinarachnius parma 

Photosensitivity, 

Photoreactivation 

Embryos 

nr (Lab) 

Cook (1968) 

Hemicentrotus 

pulcherrimus 

Chromosome 

abnormalities, 

Photoreactivation 

Gametes/ 

Embryos 

nr (Lab) 

Ejima and Shiroya 
(1982a) 

Hemicentrocus 

pulcherrimus 

Developmental 

abnormalities, 

Photoreactivation 

Gametes/ 

Embryos 

nr (Lab) 

Ejima and Shiroya 
(1982b) 

Hemicentrotus 

pulcherrimus 

Action spectrum for 
photoreactivation 

Eggs 

nr (Lab) 

Ejima etal. (1984) 

Hemicentrotus 

pulcherrimus 

Thymine dimmer 
photoreversal, 
developmental 
abnormality 

Embryos, Larvae 

nr (Lab) 

Akimoto and Shiroya 
(1987b) 

Hemicentrotus 

pulcherrimus 

Developmental 

abnormality, 

photoreversal 

Embryos 

nr (Lab) 

Akimoto and Shiroya 
(1987c) 

( continued) 
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Table 4.1 (continued) 


00 


1.4 DNA damage/repair 



UVR response 

Life-history stage 

Collection site/field 

Location 

Source 

Strongylocentrotus 

droebachiensis 

Survivorship, 

development, CPD 
accumulation 

Embryos/Larvae 

Isles of Shoals, 

Gulf of Maine, 

USA (Lab) 

Lesser and Barry (2003) 

Echinarachnius parma 

Survivorship, 

development, CPD 
accumulation 

Embryos/Larvae 

Isles of Shoals, 

Gulf of Maine, 

USA (Lab) 

Lesser and Barry (2003) 

Sterechinus neumayeri 

Developmental 

abnormalities, CPD 
accumulation 

Embryos 

Hero Inlet, Arthur 
Harbor, Palmer 
Station, Antarctica 
(Field) 

Karentz etal. (2004) 

Sterechinus neumayeri 

Developmental 

abnormalities, CPD 
accumulation 

Embryos 

Cape Armitage, 
McMurdo Sound, 
Antarctica (Field) 

Lesser et al. (2004) 

Evechinus chloroticus 

DNA repair rate, 
photoreactivation 

Embryos 

Otago Harbour, 

New Zealand 

Lamare et al. (2006) 

Sterechinus neumayeri 

DNA repair rate, 
photoreactivation 

Embryos 

McMurdo Sound, 
Antarctica (Lab) 

Lamare etal. (2006) 

Diadema setosum 

DNA repair rate, 
photoreactivation 

Embryos 

nr (Lab) 

Lamare etal. (2006) 


Miles Lamare et at. 




Strongylocentrotus 

droebachiensis 

Biological weighting 
functions, DNA 
damage 

Embryos 

Isles of Shoals, 

Gulf of Maine, 

USA (Lab) 

Lesser etal. (2006) 

Sterechinus neumayeri 

Biological weighting 
functions, DNA 
damage 

Embryos 

McMurdo Sound, 
Antarctica (Lab) 

Lesser etal. (2006) 

Evechinus chloroticus 

Biological weighting 
functions, DNA 
damage 

Embryos 

Fiordland, 

New Zealand (Lab) 

Lesser etal. (2006) 

Sterechinus neumayeri 

Developmental 

abnormalities, CPD 
accumulation 

Embryos 

McMurdo Sound, 
Antarctica (Field) 

Lamare etal. (2007) 

Tripneustes gratilla 

Developmental 

abnormalities, CPD 
accumulation 

Embryos 

Aitutaki Lagoon, 

Cook Islands (Field) 

Lamare etal. (2007) 

Evechinus chloroticus 

Developmental 

abnormalities, CPD 
accumulation 

Embryos 

Otago Harbour, 

New Zealand (Field) 

Lamare etal. (2007) 

Diadema savigtiyi 

Developmental 

abnormalities, CPD 
accumulation 

Embryos 

Aitutaki Lagoon, 

Cook Islands (Field) 

Lamare etal. (2007) 

Sterechinus neumayeri 

DNA repair 
(Photolyase 
expression) 

Gonands/Embryos/ 

Larvae 

McMurdo Sound, 
Antarctica (Field) 

Isley etal. (2009) 


Information is divided into four major responses: 1.1) behavioural, 1.2) cellular and molecular, 1.3) developmental and 1.4) DNA. Within each response, references are listed chronologically 
within each class of echinoderm. (nr = location not reported). 
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Table 4.2 Published measurement of Mycosporine-like amino acids in echinoderms. 


Species 

Location 

Habitat 

MAA type 

Tissue(s) 

Total 

concentrations 

Source 

Echinoidea 







Strongylocentrotus 

North-Eastern 

Shallow 

Mg, Sh 

Eggs 

Up to 33.21 

Adams and 

droebachiensis 

USA 

subtidal 



nmol mg -1 DW 

Shick (1996) 

Strongylocentrotus 

North-Eastern 

Subtidal 

Mg, Sh, 

Gonads, Gut, 

Up to 

Carroll and 

droebachiensis 

USA 


P-334, 

Test, Spines 

9.56 nmol mg -1 

Shick (1996) 




M-2-G, 


DW 





PI, 







A-330, Pt 




Sterechinus neumayeri 

Antarctica 

Subtidal 

Mg, Sh, 

Body Wall, 

<0.08 to 1389 

Karentz etal. 




P-334, PI 

Ovaries 

ptg g- 1 DW 

(1997) 

Strongylocentrotus 

North-Eastern 

Shallow 

Sh, P-334, 

Larvae 

~9 nmol mg 1 DW 

Adams and 

droebachiensis 

USA 

subtidal 

M-2-G, 



Shick (2001) 




PI, A-330 




Strongylocentrotus 

North-Eastern 

Subtidal 

Sh, P-334, 

Ovaries 

1.75 to 

Adams etal. 

droebachiensis 

USA 


M-2-G, 


8.86 nmol mg 1 

(2001) 




PI, A-330 


DW 
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Evechinus chloroticus 

Southern New 
Zealand 

Shallow 

subtidal 

Mg, Sh, 
P-334, PI 

Gonads/Eggs 

8.54 to 96.65 
nmol mg -1 DW 

Lamare etal. 
(2004) 

Sterechinus neumayeri 

McMurdo 

Sound, 

Antarctica 

Subtidal 

Sh, P-334, 

PI 

Embryos 

6.2 nmol mg'" 1 
protein 

Lesser etal. 

(2006) 

Strongylocentrotus 

droebachiensis 

North-Eastern 

USA 

Subtidal 

Sh, P-334 

Embryos 

650 nmol mg 1 
protein 

Lesser etal. 

(2006) 

Evechinus chloroticus 

Holothuroidea 

Southern, New 
Zealand 

Subtidal 

Mg, Sh 

Embryos 

155.66 nmol mg' 1 
protein 

Lesser etal. 

(2006) 

Twelve species of 

Great Barrier 


Mg, Sh, 

Principally in 

21.9 to 

Shick etal. 

tropical 

Holothuroids 

Reef, 

Australia 


P-334, 

M-2-G, 

PI, 

A-330, Pt 

the 

epidermis 

2053 nmol mg 1 
protein 

(1992) 

Cucumariaferrari 

McMurdo 

Sound, 

Antarctica 

Subtidal 

Mg, Sh, 
P-334, PI, 
Mg-V 

Ovaries/ 

Brood 

81-156 (§xg g _1 

DW 

McClintock 
and Karentz 
(1997) 

Holothuria atm 

Great Barrier 
Reef, 

Australia 

2 m depth 

Mg, Sh, 
P-334, 

Pt, 

M-2-G, 

As 

Epidermis, 

Viscera, 

Gonads 

0.95 to 

4.93 ptmol g _1 
DW 

Banbaranayake 
and Rocher 
(1999) 

( continued) 
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Table 4.2 (continued) 


Species 

Location 

Habitat 

MAA type 

Tissue(s) 

Total 

concentrations 

Source 

Asteroidea 

Granaster nutrix 

Antarctic 

Peninsula 

Intertidal 

PI, P-334, 

Sh, Mg, 
A-330 

Whole animal 

398 |xg g _1 DW 

Karentz etal. 
(1991) 

Odontaster validus 

McMurdo 

Sound, 

Antarctica 

Subtidal 

Mg, PI 

Ovaries, Body 
wall 

1-13 |JLg g _1 DW 

McClintock 
and Karentz 
(1997) 

Diplasterias brucei 

McMurdo 

Sound, 

Antarctica 

Subtidal 

PI 

Body wall 

30 |jtg g- 1 DW 

McClintock 
and Karentz 
(1997) 

Ophiuroidea 

Amphiplus affinis 

Antarctic 

Peninsula 

Intertidal 

PI 

Whole animal 

5 P-g g -1 DW 

Karentz etal. 
(1991) 

Crinoidea 

Promachocrinus 

kerguelensis 

McMurdo 

Sound, 

Antarctica 

Subtidal 

Mg, Sh, 
P-334, PI 

Whole animal 

775 ptgg -1 DW 

McClintock 
and Karentz 
(1997) 


Note that the total concentrations are given for a range of units depending on the study. (MAA type abbreviations: Mg — Mycosporine glycine, Sh — Shinorine, P-334 — Porphyra-334, 
PI = Palythene, Pt = Palythinol, M-2-G = Mycosprine-2-Glycine, A-330 = Asterina-330, Mg-V = Mycosporine-glycine-valine). 


Miles Lamare et at. 
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The key role ofMAAs as sunscreens for echinoderms, particularly for 
their developmental stages, has been demonstrated in a number of elegant 
experiments undertaken by Adams and Shick (1996, 2001). By dietary 
manipulation, S. droebachiensis eggs with a range of MAA concentrations 
were obtained, and exposed in vitro to UVR. The delay in egg cleavage 
associated with UVR exposure was inversely correlated with MAA con¬ 
centrations in the eggs, clearly demonstrating the protective role of these 
sunscreens (Fig. 4.3). The level of photoprotection was estimated by 
Adams and Shick (1996) based on the size of the eggs and the measured 



I-1-1-1-1-1 

0 0.2 0.4 0.6 0.8 1 

Efficiency factor for self-shading (J) at 334 nm 


Figure 3 Percentage UVR induced cleavage delay in batches of embryos from 
S. droebachiensis as a function of MAA concentration and of J (the efficiency for 
internal self-shading, calculated according to Garcia-Pichel (1994)). Embryos are 
from adults reared in controlled diets of Laminaria saccharina (^) and Mastocarpus 
stellatns (•) or a combination diet of both algae (♦). The logarithmic regression of 
percentage cleavage delay as a function of MAA concentration (y = 27.95 — 13.35 
log(.v)) and linear regression of percentage cleavage delay as a function of J 
(y = 30.79 — 26.17x) are both significant (P < 0.05).Figure is reproduced with 
permission of Dr. Nikki Adams and Professor Malcolm Shick. 
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MAA concentration, suggesting that MAAs could account for up to 86 % 
efficiency in UVR screening (Fig. 4.3). Adams and Shick (2001) noted a 
similar photoprotective role of MAAs in S. droebachiensis embryos, with 
the proportion of cultured embryos developing normally when exposed 
to UVR increasing linearly with increasing MAA concentration. 

2. 1 . 1 .2. Carotenoids Carotenoids have a number of functions in 
marine organisms, including photoprotection (Rastogi et al., 2010) via 
anti oxidation and as sunscreening agents (Mathew-Roth, 1997). 
Carotenoids are widespread in echinoderms where their antioxidant 
role is well established (Matsuno and Tsushima, 2001). However, 
studies are yet to show that they have a significant UV screening role in 
sea urchins, but a positive correlation between carotenoid concentration 
and UVR-induced egg cleavage delay was found to exist in four species of 
Strongylocentrotus (Lamare and Hoffman, 2004). It was concluded that 
this was unlikely due to a self-shading effect that they estimated was 1.5 
and 4.6% for UVB and UVA, respectively, based on egg carotenoid 
concentration and size. 

2.1. 1 . 3 . Melanin AUV screening role of melanin in some echinoderm 
species is possible, although this is yet to be confirmed. Melanin-like 
pigments or enzymes (i.e. phenoloxidase) involved in melanin synthesis 
have been described in echinoderm species across a range of classes, 
including Diadema antillarum (Jacobson and Millott, 1953; Millott, 
1953), Pisaster spp., Stichopus califomicus, Thy one briarens (Millott, 1950) 
and Holothuria spp. (Millott, 1953; Canicatti and Seymour, 1991). 
Vevers (1963) noted that the integument of echinoderms, including 
brittlestars, have ‘melanin’, suggesting the pigment is relatively 
widespread within the phyla. In D. antillarum, melanin occurs in 
epithelial chromatophores that display light sensitivity (Weber and 
Dambach, 1974). The process by which these cells respond to light has 
been described in the Diadematoid, Centrostephanuslongispinus, which also 
has pigmented, light sensitive chromatophores (Weber and Dambach, 
1974). These cells are ovoid when in darkness but extend thin processes 
on exposure to light, increasing the light absorbing area. Adams (2001) 
suggested that this process, along with seeking a cryptic habitat during the 
day, is important for species like D. antillarum to avoid excess solar 
radiation exposure. 

2.1.2. Antioxidants 

Oxidative stress occurs in echinoderms on exposure to UVR (see Section 

3.1.3, Table 4.12). Echinoids contain a range of antioxidant compounds 
including carotenoids and non-enzymatic antioxidants. Their role in 
reducing oxidative stress is discussed in Section 3.1.3. 
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2.1.3. DNA repair mechanisms 

Repair of UV-damaged DNA can be achieved by nucleotide excision 
repair (NER) and photoreactivation (PER), and although both occur in 
echinoderms, only PER has been directly observed within this group 
(Table 4.14). Photoreactivation, a light-dependent process mediated via 
photolyase enzymes (EC 4.1.99.3), is the principal mechanism of repair¬ 
ing DNA lesions (CPDs and 6-4 photoproducts). Marshak (1949) was the 
first to describe photoreactivation in echinoderms, when he noted that 
the effects of UVR on Arbacia eggs could be reduced by exposure to 
visible light. Photoreactivation has been observed in the eggs and embryos 
of number of echinoids (Wells and Giese, 1950; Cook and Rieck, 1962; 
Cook, 1968), although whether photoreactivation occurs in echinoid 
sperm is equivocal (i.e. Wells and Giese, 1950; Cook and Rieck, 1962). 

More recent work showed that photoreactivation prevented the for¬ 
mation of exogastrula or permanent blastulae in Hemicentrotuspulcherrimus 
embryos UV irradiated at the 8- or 16-cell stage (Akimoto and Shiroya, 
1987c). Embryos that were allowed to photoreactivate developed into 
normal plutei, although the repair was dependent on the time that pho¬ 
toreactivating light was delivered was only successful if the embryos were 
illuminated up to the onset of the DNA-synthesis phase of the following 
cell cycle (Akimoto and Shiroya, 1987b,c). Akimoto and Shiroya (1987a) 
also demonstrated that the quantitative and qualitative changes in the 
protein complement of the same embryos was mostly reversed by 
photoreactivation. 

Lamare et al. (2006) established that PER was the primary means of 
reversing DNA dimers, and that rates of photoreactivation in embryos and 
larvae may vary among sea urchin species across latitudes (Lamare etal. , 2006), 
with polar embryos repairing DNA slower than temperature and tropical 
species. The research showed that temperature affected the rate of repair, 
both within and between species, which was attributed to the influence of 
temperature on the photolyase enzyme and on the slower metabolic rates in 
polar larvae. Based on these findings, Lamare etal. (2007) attributed observed 
differences in the in situ accumulation of CPDs in embryos in polar, tem¬ 
perate and tropical species to differences in photoreactivation rates. Isely etal. 
(2009) partially sequenced the photolyase gene from four sea urchins, includ¬ 
ing S. neumayeri, for which the expression of the photolyase gene was 
examined during insitu exposures of embryos under a range of illuminations 
levels (i.e. with depth and sea ice cover). Expression levels were positively 
correlated with ambient UV dose, implicating that the embryos were 
responding to CPD accumulation and initiating repair. 

2.1.4. Behavioural 

Echinoderms demonstrate negative phototaxis on exposure to UVR, 
with the response observed in both larvae and adults (Table 4.11). 
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Pennington and Emlet (1986) noted that echinoid larvae ( Dendraster 
excentricus) , exposed to UVB in vitro, ceased swimming within minutes 
and settled out of the water. When UVB was blocked, the same larvae 
resumed swimming and returned to the water column within 1 h. Similar 
responses were first observed in larvae of Diademasetosum and Paracentrotus 
lividus (Fox, 1925). 

A negative phototaxis response to UVR, as well as covering behaviour 
(placing objects on dorsal body surfaces) has been documented in adult 
echinoderms (Fig. 4.4A). Covering behaviour in echinoids has been 
attributed to a number of stimuli (i.e. predator or desiccation avoidance, 
water movement) with Sharp and Gray (1962) and Lees and Carter 
(1972), the first to attribute the covering behaviour to UVR exposure 
during experimental observations on sea urchins. Adams (2001) found 
that S. droebachiensis kept in aquaria and exposed to artificial irradiation 
(containing both PAR and UVR) sought shade and covered themselves 
significantly more frequently than when irradiated by PAR only. When 
exposed to ambient irradiation in an outdoor aquarium, a similar covering 
response was observed depending on time of day, but this was typically less 
in those individuals exposed to PAR but not UVB. Verling etal. (2002) 
quantified the effect of UVR on the covering behaviour of P. lividus, and 
demonstrated that covering rate and duration increased upon exposure to 
a UVR treatment. The relative importance of UVR in covering 



Figure 4-4 Behavioural responses of echinoderms to exposure to UVR. (A) 
Covering behaviour in the sea urchin S. droebachiensis at a depth of 5—7 m. The 
animal outlined by the dash line is almost entirely covered by coralline algae 
encrusted stones, a behaviour induced by exposure to UVB radiation 
(Dumont et at., 2007). (B) The intertidal seastar, P. orchraceus, seeks shade during 
low tide. These observations demonstrate the role of solar radiation in determining 
these small-scale distributions of P. ochraceus. Confirming this, experimental trials 
showed that the sea stars moved to shade on exposure to solar radiation, although the 
response was to PAR and not UVR (Burnaford and Vasquez, 2008). Part A is 
reproduced from Dumont et at. (2007) with permission of Dr. Clement Dumont. 
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behaviour, compared with other stimuli (predation, wave motion, food, 
size) in S. droebachiensis was tested in the laboratory (Dumont etai, 2007). 
While UVR was again found to induce covering, the response was weaker 
than that caused by other stimuli such as the presence of kelp food and 
waves. 

Field observations are also consistent with solar radiation invoking a 
covering response. Covering has been observed by a number of research¬ 
ers to be greater during daylight hours (Millott, 1956; Crook etal., 1999), 
leading to their conclusion that solar radiation may play a role in their 
behaviour. Kehas et al. (2005) observed covering behaviour in the 
Caribbean sea urchin Tripneustesventricosus, and noted that covering behav¬ 
iour was greatest at shallower depths. Covering behaviour was reduced 
when UVR exposure was experimentally decreased by 66%. 
Interestingly, a small proportion of the population were albinos (4%) 
and these animal were found to have significantly greater rates of covering 
compared with pigmented individuals, possibly indicating greater sensi¬ 
tivity to solar radiation. 

Shade-seeking and covering behaviours also occur in sea stars and 
brittlestars. Drolet etal. (2004) documented shade-seeking in the brittles- 
tar Ophiopholis aculeata, although they noted that filtering UVR did not 
alter this behaviour, indicating a general response to solar radiation. 
Burnaford and Vasquez (2008) observed the intertidal sea star Pisaster 
ochraceus predominantly (85%) occurred in shaded microhabitats during 
low tide (Fig. 4.4). Confirming the role of solar radiation in determining 
these small-scale distributions of P. ochraceus, experimental trials showed 
the sea stars moved to shade on exposure to solar radiation, although the 
response was to PAR not UVR. 

These findings demonstrate the strong behavioural responses that 
UVR can exert on echinoderms, particularly those found in shallow 
water and the intertidal zone. Such behaviours should reduce the dele¬ 
terious effects of UVR, although there may be some indirect costs asso¬ 
ciated with such strategies. For example, covering and seeking shade may 
reduce movement and hence the ability to feed (Kehas et al., 2005; 
Dumont et al., 2007; Burnaford and Vasquez, 2008), may reduce gas 
exchange and could result in an energetic cost associated with moving 
and holding objects (Kehas etai, 2005; Sigg etal., 2007). Consequently, 
reduced movement could decrease the ability of animals to escape preda¬ 
tion. Although yet to be quantified, given the key role that echinoderms 
can play in their environment via movement and feeding, the influence of 
UVR on behaviour may have wider implications for the community. For 
example, Burnaford and Vasquez (2008) noted the feeding behaviour, and 
hence structuring role of the intertidal P. ochraceus may, in part, be con¬ 
trolled by exposure to solar radiation, and at a level greater than has been 
previously appreciated (Fig. 4.4B). 
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UVR-Induced Damage in Echinoderms 


3.1. Cellular 


3.1.1. DNA damage 

Damage to DNA in the form of DNA dimers, strand breaks, linkages, 
DNA cross-overs and changes in hydration have been observed in echi¬ 
noderms upon exposure to UV radiation. DNA dimer (CPDs and 6-4 
photoproducts) formation is a key response to UVR exposure and has 
been observed in echinoderms (Table 4.Id). Despite their importance 
however, DNA dimers have only been quantified in echinoids and only in 
the embryonic and larval stages. Lesser etal. (2006) explored the effects of 
UVR (including CPD formation) in the laboratory on the embryos of 
three echinoid species including two temperate and an Antarctic species. 
CPD concentrations, expressed as relative absorption units, were quan¬ 
tified in embryos exposed under long-pass filters that provided five UVR 
treatments within the UVR spectrum. While all three species showed 
CPD formation in the UVA portion of the spectrum, CPD formation was 
greatest in the UVB wavelengths and decreased progressively as the 
shorter wavelengths were blocked. From these data, biological weighting 
functions (BWFs) were developed from which two important conclu¬ 
sions were drawn: (1) the Antarctic species (S. neumayeri) showed the 
greatest sensitivity to UVR compared with the temperate species 
( E. chloroticus and S. droebachiemis) which (2) corresponded with S. neumayeri 
having the lowest sunscreen compound concentrations. Based on the 
BWFs and measures of ambient spectral irradiances, the authors concluded 
that E. chloroticus may receive the greatest biologically effective UVR in situ, 
although later research by Lamare et al. (2007) showed that CPD accumu¬ 
lation among these species may not reflect these laboratory sensitivities. 
Lamare et al. (2006) also observed significant accumulation of CPDs 
in embryos from four species exposed to UVB for 1 h, during in vitro 
experiments examining DNA photorepair rates. 

DNA dimer accumulation has also been observed in echinoid larvae 
exposed to UVR insitu. In experiments aimed at examining the effects of 
changes in UVB (as a result of ozone depletion) and sea ice coverage on 
sea urchin embryos, Lesser etal. (2004) measured relative levels of CPD 
accumulation in larvae of the Antarctic sea urchin, S. neumayeri exposed to 
ambient UVR under sea ice over two years. CPD accumulation was 
correlated with the UVB dose received, which was found to be a function 
of depth, UV-filter treatment and overhead ozone conditions (which 
varied between the 2 years of the research). Greater rates of embryo 
abnormality were associated with treatments that resulted in accumulated 
CPDs. Karentz etal. (2004) also examined the effects of Antarctic ozone 




Ultraviolet Radiation and Echinoderms: Past, Present and Future Perspectives 


169 


depletion on the embryos of S. neumayeri, including quantifying the 
accumulation of CPDs. CPD concentrations were as high as 
272 CPDs mb 1 DNA after 6 days of exposure. Concentrations of 
CPDs varied with the UVB dose received, although the relationship 
between the two was non-linear with a threshold dose of 25 kj m _ . 
The researchers also identified an ancillary threshold limit of 
17 CPDs mb DNA, above which normal development of embryos 
was unlikely. 

To further understand the relative sensitivity of echinoid embryos to 
UVR, and to quantify natural levels of CPD accumulation in species from 
a range of latitudes, Lamare et al. (2007) undertook a series of in situ 
exposures that included testing the effect of depth and UVR dose on 
CPD accumulation. Maximum CPD concentrations were greatest in the 
Antarctic S. neumayeri embryos (30.8 CPDs mb 1 DNA) compared with 
the temperate ( E. chloroticus) and tropical (Diadema savignyi) species (16.5 
and 3.0 CPDs mb 1 DNA, respectively). Within each species, CPD 
concentration was greater at the shallower depth, and increased with 
UVB dose (but not UVA dose). In a similar set of experiments, Lesser 
(2010) explored depth-dependent DNA damage in embryos of the sea 
urchin, S. droebachiensis in the Gulf of Maine, which included an exam¬ 
ination of the relative levels of CPD accumulation and levels of oxidative 
modification of nucleotide bases (8-hydroxy-2 / -deoxyguanosine (8- 
OHdG)). Lesser (2010) showed that DNA damage occurred in embryos 
exposed to UVR at depths <5 m despite the presence of UVR absorbing 
sunscreen compounds. By examining different forms of DNA damage, 
the research highlighted the differential effects of UVB and UVA, which 
caused direct and indirect damage to DNA, respectively. Importantly, 
Lesser (2010) noted that DNA damage was mostly the result of exposure 
to UVA, reflecting its greater irradiances and depth penetration compared 
to UVB. When biologically weighted for DNA damage (using 
Lesser et al., 2006), the importance of UVA was even more pronounced. 

Research by Lesser (2010) highlights the limited range of measure¬ 
ments of DNA damage that have been made in echinoderms. With the 
exception of 8-OHdG concentrations previously mentioned (Lesser, 
2010) and strand excisions in coelomocytes of the sea urchin P. lividus 
(Schroder et al., 2005), little is known of other indirect effects of UVR 
on DNA such as cross-linking and singe-strand breaks. In the wider 
context, measurements of DNA damage have been restricted to echi- 
noids, and with the exception of coelomocytes, within the embryological 
and larval stages. To fully understand the effects of UVR on DNA 
damage, a greater suite of measurements may be necessary. 

The measurements of DNA damage (CPD concentrations) reported 
in echinoderms are within the range reported for other marine animals 
exposed to ambient light conditions (i.e. 3—30.8 CPD mb -1 DNA, 
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Dahms and Lee, 2010). The exception is the high levels (up to 
272 CPD mb 1 DNA) reported in Antarctic embryos (Karentz et al., 
2004; Lamare etal., 2007). Lamare etal. (2007) concluded that the greater 
accumulation of CPDs in polar embryos may be the result of lower 
concentrations of sunscreens and slower rates of DNA repair within the 
species. 

3.1.2. Protein and lipid damage 

Little is known about the impacts of UVR on protein and lipid biochem¬ 
istry, and metabolism (Table 4.12). Akimoto and Shiroya (1987a) used 
two-dimensional electrophoresis to determine the patterns of accumu¬ 
lating proteins in normal, UV-irradiated and photoreversed sea urchin 
(H. pulcherrimus) embryos. In this study many qualitative and quantitative 
variations in the accumulating proteins were observed. More recently, a 
proteomic investigation of Strongylocentrotuspurpuratus embryos exposed to 
UVR has shown differential protein spot migration between UV-pro- 
tected and UV-treated embryos at 30 and 90 min post-fertilisation 
(Campanale etal., 2011). The study noted that proteins from multiple 
pathways were affected by UVR, including cellular stress, protein synthe¬ 
sis, signal transduction and cytoskeletal dynamics. The further application 
of proteomic, lipidomic and metabolomic profiling techniques could pro¬ 
vide, as they have in studies on human cells exposed to UVR (Huang 
etal., 2005; Reich etal., 2010), useful information on protein and lipid, 
damage, repair and metabolism in echinoderms exposed to UVR. 

3.1.3. Oxidative stress 

Reactive oxygen species (ROS) are a group of highly reactive oxygen 
intermediates that include superoxide (0 2 ~), hydrogen peroxide (H 2 0 2 ) 
and the hydroxyl radical (HO'). ROS are produced by the partial single¬ 
electron reduction of oxygen or by the alteration of oxygen electron spin 
states by photoactivation resulting in the production of singlet oxygen. In 
living cells 0 2 is also produced as a by-product of mitochondrial elec¬ 
tron transport, associated with normal respiration (Halliwell and 
Gutteridge, 2007). Under conditions of normality organisms are exposed 
to relatively low and steady ROS levels and under these conditions the 
cellular antioxidant defences, which include non-enzymatic antioxidants 
such as glutathione and Q'-tocopherol, and enzymatic antioxidants such as 
catalase (CAT), glutathione peroxidase (GPx), glutathione reductase 
(GR) and superoxide dismutase (SOD), are able to regulate ROS levels, 
minimise any potential oxidative damage, and maintain a balanced and 
stable cell redox potential. However, exposure to pollutants, rapid changes 
in environmental condition such as temperature and or oxygen levels, 
high levels of visible light and exposure to UVR can lead to the over 
production of ROS and lead to oxidative damage (Smith et al., 2000; 
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Burritt and MacKenzie, 2003; Lesser, 2006; Halliwell and Gutteridge, 
2007; Burritt, 2008). While most cellular macromolecules are sensitive to 
oxidation, DNA, proteins and lipids are arguably the most critical cellular 
macromolecules that can be damaged due to oxidation. Oxidative damage 
to DNA is often measured as 8-OHdG, to proteins as protein carbonyls 
and to lipids directly as peroxides or indirectly as malonyldialdehyde 
(MDA). All the above are considered to be sensitive biomarkers of oxi¬ 
dative stress and have been used to measure the impact of environmental 
stressors, including UVB, on a wide range on organisms (Lesser, 2006; 
Halliwell and Gutteridge, 2007). 

Numerous studies on organisms have demonstrated that exposure to 
UVR can cause oxidative stress (Lesser, 2006; Halliwell and Gutteridge, 
2007). With respect to echinoderms UVR-induced oxidative stress can in 
theory occur at any stage in the life cycle of the organism, from gametes 
and larvae through to adults, however due to their relative transparency it 
is the gametes and larvae that are probably the most sensitive to UV- 
induced oxidative damage. Lu and Wu (2005a), in a study on the impact of 
environmentally relevant levels of UVA and UVB on the sperm of the sea 
urchin Anthocidariscrassispina, showed that UVA and UVB both enhanced 
the production of ROS and greatly increased LPO. Lu and Wu (2005a) 
also concluded that the increase in ROS production and oxidative damage 
due to UVR exposure could account for the decline in sperm motility 
and fertilisation success that they also observed. They suggested that ROS 
could directly attack the mitochondrion of the sperm, and could cause 
significant membrane disruption. They also speculated that elevated 
UVR levels could pose a significant threat to the ability of natural popula¬ 
tions of sea urchins to reproduce successfully in shallow waters. 

In a field experiment conducted on the embryos of the green sea 
urchin S. droebachiensis. Lesser (2010) demonstrated the importance of 
measuring not only direct damage to DNA caused by UVR exposure, 
in the form of CPDs, but also indirect oxidative damage to DNA, in the 
form of 8-OHdG. In this study despite embryos containing significant 
levels ofmycosporine-like amino acids that absorb UVR, primarily in the 
UVA portion of the spectrum, and increased levels of SOD, significant 
increases in both CPDs and 8-OHdG were observed at depths less than 
5 m. Lesser (2010) suggested that screening alone did not provide ade¬ 
quate protection from DNA damage in the field at depths less than 5 m. 
Perhaps of greater significance was the finding that while exposure to 
UVA did not result in the formation of CPDs, exposure to UVA greatly 
increase oxidative DNA damage at depths of 5 m or less. Other studies 
have also shown UVA can cause oxidative damage to cells (Reinheckel 
etai, 1999; Hoerter etai, 2008). 

Two recent studies have shown that oxidative stress caused by UVR 
occurs irrespective of echinoid species or global location. Lister et al. 
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Figure !\. 5 Changes in oxidative stress in sea urchin embryos in response to ozone 
depletion over the Antarctic (from Lister et al., 2010a). (A) Ozone concentrations 
over the Antarctic during field experiments indicating periods (26 October 2008) of 
normal ozone concentrations overhead and periods of ozone depletion (3 November 
2008) over the experimental site (indicated by a black dot). (B) Experimental racks at 
1 and 4 metres depth on which Sterechinus neumayeri embryos were exposed to 
ambient and filtered (UVB-minus and UVR-minus) light. (C) Oxidative damage 
(mean (d=SE) protein carbonyl and lipid hydroperoxide concentrations) in embryos 
under normal and reduced ozone concentrations. (D) Change in antioxidant enzyme 
activities (mean (d=SE) superoxide dismutase and catalase activity) in embryos under 
normal and reduced ozone concentrations. For all bar graphs, grey bars represent 
high ozone conditions, and white bars indicate low ozone conditions. Figures 
reproduced from Lister et al. (2010a). Ozone maps courtesy of NASA (USA). 
Photograph (B) used with permission of M.F. Barker. 
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(2010a) conducted an investigation into the influence ofincreased UVR, 
caused by the ‘ozone hole’, on the Antarctic sea urchin S. neumayeri over 
two consecutive 4-day periods in the spring of2008. They found that the 
presence of the ozone hole and the corresponding increase in UVB 
exposure caused an increase in oxidative damage to lipids and proteins, 
measured as lipid peroxides and protein carbonyls and developmental 
abnormality in embryos grown in open waters (Fig. 4.5). Their results 
also showed that while the embryos had a limited capacity to increase the 
activities of protective antioxidant enzymes, severe oxidative damage 
occurred following exposure to UVB. More importantly, this study also 
demonstrated that the effect of the ozone hole was largely mitigated by sea 
ice coverage, with embryos under the sea ice largely protected from UV- 
induced oxidative damage. In a second field study Lister et al. (2010b) 
investigated oxidative damage and antioxidant defences in larvae of the 
tropical sea urchin Tripneustesgratilla, at two depths, in Aitutaki in the Cook 
Islands. Compared to those shielded from UVB, larvae exposed to UVB 
held at a depth of 1 m showed increased developmental abnormality and 
increases in oxidative damage to proteins, despite increased levels of CAT, 
GPx, GR and SOD. Increased oxidative damage was also observed in 
larvae held at 4 m, but the amount of damage was much less. However, 
unlike the Lesser (2010) study, both the above studies provided no con¬ 
clusive evidence that UVA caused any significant increases in oxidative 
damage or abnormality. Campanale etal. (2011) provides further support 
for the importance of antioxidant defences in echinoids exposed to UVR, 
by showing that proteins involved in redox regulation (e.g. glutathione 
peroxidase) increased in UV-treated S. purpuratus embryos. 

3.1.4. Cell cycle and stress responses genes 

Lesser and Barry (2003) postulated that cell-cycle arrest at critical DNA 
check points following UVB-induced DNA damage was a component of 
the observed developmental delays in S. droebachensis embryos. Their 
conclusion was based on the observation that survivorship of earlier staged 
embryos was more sensitive to DNA damage than older embryos, which 
may reflect lower repair capacities and subsequently greater levels of 
apoptosis mediated by cell-cycle genes (p53). In support of this, the 
process of apoptosis in response to UVB exposure has been observed in 
sea urchins, with coelomocytes from P. lividus (Matranga et al., 2006) 
displaying apoptotic nuclei after short-term UVB exposure in the 
laboratory. 

In the same study, Matranga etal. (2006) showed that levels of the stress 
protein hsp70 in coelomocytes respond to UVB in a dose-dependent 
manner. An increase in hsp70 expression on exposure to UVR also occurs 
in the embryonic stages of echinoids. Bonaventura et al. (2005) noted a 
dose-dependent increase in the expression of the protein in P. lividus 




174 


Miles Lamare ef at. 


blastula, as well as a second stress marker p38 MAPk. Bonaventura etal. 
(2006) subsequently showed that embryos at the 16-cell stage also 
increased their levels of hsp70 mRNA and the protein itself in a UVB 
dose-dependent manner. P. lividus embryos also increase the expression of 
stress-response proteins from the 14-3-3 protein family after exposure to 
UVB (Russo etal., 2010). Whole-mount in situ hybridisation of the 14-3- 
3 transcripts showed it was locally distributed in the embryonic archen- 
teron and throughout the germ layers. Importantly, these distributions 
mapped cell-cycle disruption (i.e. failure of the archenteron to elongate). 

These studies highlight an emerging appreciation of the range of 
cellular stress responses and regulatory cascades that likely occur in echi- 
noderm tissues and developmental stages on exposure to UVR. These 
responses are dose-dependent, and their initiation will depend on the 
ratio of damage to repair within tissues and developmental stages. While 
studies have quantified levels of damage under ambient conditions (see 
previous sections), the degree that cell regulatory processes are responsible 
for the developmental changes seen during in situ experiments, or in the 
process of mitigating damage, is unknown. 

3.2. Organismal damage 
3.2.1. Developmental biology 

The effects of UVR on developmental biology have been examined 
extensively using echinoderms, and particularly echinoids, as model sys¬ 
tems (Table 4.13). Gametes, fertilisation and embryonic and larval devel¬ 
opment are all sensitive to UVR, with a range of responses observed after 
exposure to UVB and/or UVA, both in the laboratory and under ambient 
field conditions. 

Sperm are the most UVR sensitive gametes in echinoderms with the 
detrimental effects including impaired swimming (Au etal., 2002; Lu and 
Wu, 2005a), oxidative stress (Lu and Wu, 2005a), a loss of mitochondrial 
function and membrane integrity (Lu and Wu, 2005b; Pruski etal ., 2009), 
chromatin and DNA damage (Pruski et al., 2009) and morphological 
alteration (Pruski etal., 2009). A key outcome of the exposure of sperm 
to UVR is a reduction in fertilisation rate, which occurs in a dose- 
dependent manner. Lu and Wu (2005b) demonstrated an exponential 
decrease in fertilisation rates across environmentally relevant doses for 
UVB (91—98% reductions) and UVA (38—72% reductions). Similarly, 
Nahon et al. (2009) observed that rates of fertilisation decreased from 
>90 to 20% when sperm were exposed to a range of environmentally 
relevant UVR irradiances. 

The damage seen in sperm reflects a lack of sunscreens, reduced DNA 
repair and antioxidants, a high surface to volume ratio and high levels of 
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polyunsaturated fatty acids (PUFA) in the membranes in these cells (Lu 
and Wu, 2005a) making them vulnerable to lipid peroxidation (Lu and 
Wu, 2005a; Pruski etal., 2009). In contrast to sperm, initial fertilisation 
rates do not decrease when eggs are exposed to UVR prior to insemina¬ 
tion (Ikeda, 1965; Nahon et al., 2008, 2009), although the processes 
during fertilisation can be altered by irradiation. Rustad (1959) noted 
the lifting of the fertilisation envelope was suppressed on the UV-irradi- 
ated side of S.purpuratus eggs, and was inhibited under highest doses, with 
Whalley etal. (1995) suggesting this represents an interruption exocytosis 
of the cortical granule. 

Pronuclear fusion may also be influenced by UVR, with sperm aster 
and female pronuclei migration responding to regional UV illumination 
in the fertilised egg (Hamaguchi and Hiramoto, 1986). Cook (1968) 
noted in earlier work that pronuclear fusion timing was not affected by 
UVR exposure, although they noted that prophase was delayed. 
Irradiation can increase rates of polyspermy in sea urchin and sea star eggs 
(David et al., 1985), which was attributed to the destruction of Na + 
channels and a loss of fast-block capacity. 

Post-fertilisation developmental effects of exposing eggs to UVR are 
evident and well documented. Outcomes of laboratory irradiances 
include delays in cleavage rate (Giese, 1939; Cook, 1968; Rustad, 
1971; Adams and Shick, 2001; Lamare and Hoffman, 2004; Nahon 
etal., 2009), reduced survivorship (Lesser and Barry, 2003; Nahon etal., 
2009) and increased rates of developmental abnormalities in eggs and 
larval stages (Akimoto and Shiroya, 1987b,c; Adams and Shick, 2001; 
Lesser and Barry, 2003; Bonaventura et al., 2006; Lamare et al., 2006; 
Lesser et al., 2006; Nahon etal., 2009) (Fig. 4.6). 

Eggs, embryos and larvae undergo a range of detrimental changes on 
exposure to UVR which may include damage and reduced migration of 
primary mesenchyme cells (Akimoto and Shiroya, 1987c), impaired skel¬ 
etal formation (Bonaventura etal., 2005; Nahon etal., 2009) and apoptosis 
(Lesser etal., 2003). The effects on embryonic development of UV treat¬ 
ment are demonstrable in echinoderms, characterised by reduced devel¬ 
opment, abnormal cell division, thickening of body walls and cellular 
packing of the blastocoel (Fig. 4.6). At the molecular level, damage 
includes DNA lesions (Lesser and Barry, 2003; Lesset et al., 2006; 
Lamare etal., 2006; Nahon etal., 2008), oxidative stress (Lesser, 2006), 
reduced expression of a gene (P/-SM30) coding for skeletal matrix pro¬ 
teins (Bonaventura etal., 2005), expression of heat shock proteins and 
stress markers (Bonaventura etal., 2005, 2006) and changes in protein 
profiles (Campanale et al., 2011) including a decrease in actin and tubulin 
(Akimota and Shiroya, 1987). 

Many of these processes have been observed in embryos exposed to 
ambient UVR in field experiments across a range of latitudes. Eggs and 
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Figure 4.6 The effects of UVR on the development of P. Hindus embryos 24 h after 
exposure to UVB doses of200 J m 2 (B), 400 J m -2 (C), 800 J m -2 (D) and no dose 
(A). Photographs are from Bonaventura et al. (2006), and indicate the typical 
responses to UVR observed during development, namely reduced size, delayed 
development, abnormal cell division patterns, thickening of the blastoderm and 
the accumulation of cells in the blastocoel. Scale bar = 50 |xm. Images reproduced 
with permission of Dr. Rosa Bonaventura. 


embryos of the Antarctic sea urchin S. neumayeri show abnormal devel¬ 
opment (Lesser et al., 2004; Karentz et al., 2004; Lamare et al., 2007), 
oxidative stress and an increase in antioxidant activities (Lister et al., 
2010a), DNA damage (Karentz et al., 2004; Lesser et al., 2004; Lamare 
etal., 2007) and expression of DNA repair enzymes (Isely et al., 2010). 
Oxidative stress and DNA damage occur in temperate echinoid embryos 
and larvae exposed to UVR (Lamare et al., 2007; Lesser, 2010), with 
similar responses observed in the developmental stages of tropical species 
exposed in lagoonal environments (Lamare et al., 2007; Lister et al., 
2010 b). 
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These findings point to the important role that UVR may play during 
reproduction and development of echinoderms. This is despite eggs and 
embryos being well equipped to mitigate the effects of solar radiation, 
either via repair and protection mechanisms (DNA repair pathways, heat- 
shock proteins, sunscreens and antioxidants) or possibly by behavioural 
responses (see previous sections). The often shallow distribution of 
embryos and larvae in the water column will also contribute to their 
vulnerability to UVR (Rumrill et al., 1985; Ernlet, 1986; Pennington 
and Ernlet, 1986; Pedrotti, 1990; Sewell and Watson, 1993; Miller and 
Emlet, 1997; Lamare, 1998), with echinoderm larvae probably not able to 
control their vertical location and therefore more likely to be exposed to 
damaging levels of UVR in their natural environment. Similarly, repro¬ 
duction in intertidal and shallow subtidal species will almost inevitably 
involve exposure of gametes to potentially damaging levels of solar 
radiation. 

3.2.2. Mortality and population level effects 

Very little is known of the effects of UVR on post-settlement echino¬ 
derms in terms of mortality (Table 4.1c). Exposure of the brittlestar 
Ophiodermabrevispinum to ambient solar radiation in experimental aquaria 
resulted in almost 100% mortality after 4 days (Johnsen and Kier, 1998). In 
treatments where UVR was blocked by filters, no mortality occurred. In 
contrast, Raeka-Kudla etal. (1993) found no effect of ambient UVR on 
the survival of two brittlestars, Ophio nereis reticulata and Ophiocomaechinata, 
or in the sea urchin Echinometra viridis exposed in experimental aquaria. 
Species from both studies were tropical, with no information available for 
how temperate and polar species survive on exposure to UVR. 



. Future Research 


In reviewing research on the effects of UVR on echinoderms, a 
number of patterns in relation to the scope and focus of these investiga¬ 
tions have emerged. These patterns are insightful, particularly as a means 
of providing a guide to future research for echinoderm photobiologists. 
We give three of the potential directions for this research. 


4.1. A greater diversity of UVR research within 
the echinoderm phyla 

While a body of literature exists on UVR effects on echinoderms, there 
are a number of limitations within it. Firstly, of the published studies, 
approximately 80% are on echinoids, with remaining classes the focus of 
only 14 peer-reviewed studies. Furthermore, within this literature, 75% 
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of studies examined UVR responses in the embryological or larval stages, 
with little attention given to the post-settlement life-history stages. This 
bias is for good reason; sea urchins are a widely used model system in 
developmental research and findings from research on echinoid embryos 
are generally applicable to a wide range of species, beyond the echino- 
derm phyla. It does mean, however, that we have a limited understanding 
of the effects of UVR within the Echinodermata for both non-echinoid 
species and for their post-settlement stages. 

An understanding of the effects of UVR on these echinoderms at the 
population level, both in terms of their dynamics and ecology is also 
worthy of exploration. This is important as many post-settlement echi- 
noderm populations have a keystone role within their ecosystem. The sea 
star P. ochraceus is such a species as it can control the vertical distribution of 
intertidal communities (Paine, 1974), and Burnaford and Vasquez (2008) 
have shown in the laboratory, that UVR can influence the distribution 
and movement of this organism. P. ochraceus seeks shelter from solar 
radiation in the field and, hence by influencing the seastar, UVR con¬ 
tributes to key prey/predator interactions in this habitat. Similar interac¬ 
tions in echinoids and sea cucumbers (important grazers and bioturbators) 
probably occur (i.e. Dumont et al., 2007), along with a suite of other 
ecological processes associated with movement for which UVR may 
influence yet we know little of. 

Similarly, there has been little linkage or application of developmental 
studies to post-settlement processes. While many laboratory and field 
studies have shown that UVR is detrimental to development, the impli¬ 
cations for recruitment and the population dynamics of echinoderm 
species is yet to be quantified. A good example is the Antarctic, where 
research has focused on the effects of ozone depletion and increased UVR 
on embryos and larvae of S. neumayeri. While this work has shown that the 
ozone hole results in DNA damage, oxidative stress, mortality and devel¬ 
opmental delays, there is, as yet, no evidence to suggest whether there 
have been changes at the population level during 40 years of ozone 
depletion. Researchers are aware of these issues; Lesser et al. (2004) dis¬ 
cussed the potential for UVR-induced recruitment failure of S. neumayeri 
in McMurdo Sound over this period and suggested monitoring of popu¬ 
lations over the long term would shed light on this interesting question. 

A second limitation is in the range of responses examined in echino¬ 
derms. As model organisms, echinoderms can readily be used to quantify 
the effects of UVR across a suite of molecular and cellular processes. For 
example, research by Lesser (2010) highlights the limited range of mea¬ 
surements of DNA damage that have been made in echinoderms and 
notes that damage can mostly be the result of exposure to UVA and 
8 -OHdG formation. Similarly, repair of DNA can be via photoreactiva¬ 
tion, DNA glycosylases and NER, the latter two of which have not been 
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examined in echinoderms. The same limitations could be noted for a 
range of possible molecular and cellular responses to UVR such as cell 
signalling, stress responses and apoptosis. Insight into these processes 
could be forthcoming with the application of modern molecular techni¬ 
ques such as proteomics (Campanale et al., 2011), genomics and DNA 
microarrays (i.e. Hofmann and Gaines, 2008) now possible with the DNA 
sequencing of species within the Echinodermata (i.e. S. purpuratus). 

4.2. Synergistic effects 

Research on the effects of UVR on echinoderms has generally been 
conducted in isolation and does not consider the potential impacts of 
other environment parameters. In reality, UVR acts in concert with a 
range of abiotic and biotic factors that could be independent of, or 
interactive with UVR, either in an additive, synergistic or antagonistic 
manner. The cumulative impacts on marine systems from human activ¬ 
ities have shown these interactions occur at approximately equal propor¬ 
tions when two or more stressors are examined (Crain et al., 2008). 
Understanding such interactions is especially relevant when considering 
future climate change scenarios (Todgham and Hofmann, 2009), where 
increases in UVR may be accompanied by ocean acidification, rising sea 
temperatures, reduced sea ice cover, altered salinities and changes in 
productivity (reviewed by Hoegh-Gulberg and Bruno, 2010). Hader 
et al. (2007) provided a perspective on the effects of UV radiation on 
aquatic systems and interactions with climate change, noting that at the 
organism level UVR often acts as an additional stressor while at the 
ecosystem level the effect of UVR and climate change processes can be 
more pronounced, particularly on community and trophic structure. 
Indeed, research on the interactive effect of UVR and ocean acidification 
is underway on other marine groups (Gao etal., 2009; Gao and Zheng, 
2010), and past studies have shown that exposure to UVR exacerbates the 
detrimental effects of reduced pH (and vice versa) on marine algae. 
Echinoderms may provide excellent model systems to further examine 
these processes in marine metazoans. 

Equally important is the interaction of UVR with pollutants. 
Steevens etal. (1999) examined the combined effects of UVB and poly¬ 
cyclic aromatic hydrocarbons (PCBs) on Lytechinus variegatus embryos, 
noting an additive toxicological effect. In contrast, the amount of DNA 
strand breakage in sea urchin coelomocytes following UVR exposure 
was, unexpectedly, reduced when combined with toxic concentrations of 
cadmium (Schroder et al., 2005). The interactive effect of UVR and 
pollutants may be intrinsic (i.e. cadmium-induced inhibition of the 
DNA excision repair enzyme, endonucleases in the Schroder et al., 
2005 study) or extrinsic (i.e. changes in the pollutant chemistry and its 
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toxicity by UVR-induced photochemical processes). The toxicological 
effects of a range of pollutants such as metals, hydrocarbons and PCBs 
have been studied extensively in echinoderm gametes and embryos 
(Pagano et al., 1997; Coteur ef al., 2004; Manzo ef al., 2010). Such 
approaches, combined with UVR as a factor, provide an opportunity 
to expand our understanding of the interaction of these pollutants under 
real-world conditions were UVR could play an important synergistic 
role. 


4.3. UVR and life histories 

UVR has had a major influence on the evolution of life, including its role 
in promoting oxygen toxicity (i.e. oxidative stress) in metazoans (Lesser, 
2006). In recent years ecological research has focused on the relevance of 
antioxidants and oxidative stress in the evolution of life-history strategies 
and physiological trade-offs in animals. Managing oxidative stress is likely 
to be a major determinant of life histories, as virtually all metabolic 
activities generate ROS, the levels of which also increase in organisms 
exposed to abiotic and/or biotic stressors (Monaghan et al., 2009). 
Important considerations are that the levels of oxidative stress incurred 
by organisms are not constant but vary with developmental stage, envi¬ 
ronmental conditions and levels of activity (Monaghan et al., 2009). 
Echinoderms have been used extensively to theorise on life-history strat¬ 
egies (Hart, 2002) and could provide ideal model systems to address 
questions relating to the role of UVR and oxidative stress in shaping 
developmental patterns. 



. Conclusions 


The influence of UVR on echinoderms has been researched for 
more than 100 years and has increased our understanding of the role of 
UVR in the marine environment. While focused mainly on echinoids, 
and particularly their developmental stages, this research has added to our 
knowledge in a number of areas: 


• The effect of UVR on developmental biology, and specifically on 
gametes, embryos and larval stages. 

• The physiological responses and mitigating strategies occurring in ani¬ 
mals exposed to UVR, particularly the role of sunscreens, DNA damage 
and repair and oxidative damage. 

• The types and level of UVR damage occurring in situ, including an 
understanding of the effects of ozone depletion on marine ecosystems in 
polar regions. 
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Looking forward, future research on echinoderms has the potential to 
make a significant contribution to photobiology at a number of levels: 

• Echinoderms, particularly echinoids, as model systems for developmen¬ 
tal biology will continue to add to our understanding of the effects of 
UVR on development, particularly with the application of new molec¬ 
ular techniques. 

• In a changing marine environment, the role of UVR and its synergistic 
relationship with other climate change agents, will need to be addressed. 
With a good general knowledge of the effects of UVR on echinoderms 
coupled with a broad knowledge of developmental and physiological 
processes in the phyla, echinoderms could be at the forefront of this 
research. 

• Fundamental and theoretical biological questions relating to photobi¬ 
ology (such as the role of UVR shaping life histories) would be ideally 
asked in the Echinodermata. 
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Plate 1 Phylogenetic tree of the p53 superfamily. Sequences corresponding to the 
conserved DNA binding domain were used to generate an Unrooted Phylogenetic 
Tree using the Splits Tree 4 program (Huson and Bryant, 2006). Abbreviations used 
are Ag, Anopheles gambiae ; Ap, Acyrthosiphon pisum; Bf, Branchiostoma jloridae 
(Amphioxus); Cbr, Caenorhabditis brenneri; Cbri, Caenorhabditis briggsae ; Ce, 
Caenorhabditis elegans ; Ci, Ciona intestinalis ; Co, Capsaspora owczarzaki ; Cp, Capitella 
sp.; Cq, Culex quinquefasciatus; Cre, Caenorahbditis remanei; Cs, Ciona savignyi ; Dp, 
Daphnia pulex; Dm, Drosophila melanogaster ; Dr, Danio rerio (Zebra fish); Hr, Helobdella 
robusta; Hs, Homo sapiens ; Lf, L. forbesi; Lg, Lottia gigantea; Ma, M. arenaria ; Mb, 
Monosiga brevicollis ; Me, M. edulis; Nav, Nasonia vitripenes; Nv, N. vectensis; Ph, 
Pediculus humanus ; Pp, Pristionchus pacificus (Elephant shark); Leuconraja erinacea 
(Little skate); Sp, Strongylocentrotus purpuratus ; Ta, Tnchoplax sp.; Tc, Tribolium 
castaneum ; Tp, Trichonella sp.; TET, tetramerization domain; SAM, sterile alpha 
motif; TS, tetramerization domain and SAM domains. Reprinted from Rutkowski 
et al. (2010) with permission from CSH Press. 
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MapS3/63/73 Core: 


MeONp€3/73: MiwtRTGmYR 
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QIVPTLEKDNQVLGRRCVEVRICACPGRDRKADERGSLPPMVSGGVKKSQMPKPSMGTEITTVSSGKKRKPEDDEQTrTLTVRGREMYOMLCKIRDS 


LErAALLPOHQLQSLKQKQVEVQRO 


Map63/73 Extension: 


EQNLYSLLQLODFS LDDLAKHKIGNSRRNXIWKSLLELRNQGFTTAESQD5LAKQAS5TSTIS LASQGSISQNS 7YNPGF YEV7RHTFKHTISLTKE 
ERHVOTRSSASKAD 


Map53 Extension: 


WLTNVLAKEGKSRLIKKKHRPGKIIRHPLK 


Map53/63/73 intron 3 ..ataaataaaaatgattaagtttgagagaactggcttcacaacttatagggtgagtacgc... 

, mu n miimiimi n mu mm 

McDNp63/73 mRNA ATGATCAAATTTGAGAGAACTGGATTTACAACCTATAGG 

MIKFERTGFTTYR 


Plate 2 (a) Genomic structure ofp53/63/73 in M. arenaria. Boxes represent exons 

and V indicates intron regions (not drawn to scale). The 5 , -UTR is grey, the core 
sequence for both clam p53 and p63/73 proteins is green, the p53 extension is blue 
and the p63/73 extension is red. The yellow arrow denotes the location of a 
predicted DN isoform; (b) amino acid sequences for clam p53 and p63/73 with 
intron/exon boundaries marked (▼, clam p53; human p53). Functional 
domains are indicated by horizontal lines above the sequence. Conserved amino 
acids implicated in Mdm2 binding are underlined. Sequence highlighted in yellow 
denotes the predicted DN isoform. Genomic DNA was purified from flash-frozen 
gill tissue of M. arenaria as previously described (Barker et al ., 1997). DNA template 
was prepared using the Universal Genome Walker Kit according to the 
manufacturer’s instructions (Clontech). Clam sequence-specific PCR. primers 
were designed based on the clam p53 and p63/73 coding sequences (GenBank 
Accession numbers AF253323 and AF253324, respectively) and used in 
conjunction with adapter primers (Clontech). Amplified products were cloned 
into the pCRII vector (Invitrogen) and sequenced in both directions at the 
University of Maine DNA Sequencing Facility. Contiguous sequences were 
assembled using Sequencher version 4.1.2 software (Gene Codes Corp., Ann 
Arbor, MI) and the open reading flames identified; sequence for genomic DNA 
was deposited under GenBank Accession number FJ041332 and (c) sequence for M. 
edulis DNp63/73 (DQ865151.1) was aligned to FJ041332 using Spidey (http:// 
www.ncbi.nlm.nih.gov/spidey/) to elucidate a sequence and intron/exon 
boundaries for a predicted DN isoform in M. arenaria , and was found to be 87% 
identical to M. edulis at the nucleic acid level and 100% at the protein level and located 
in intron 3. 












Plate 3 Structural analysis of M. arenaria p53 shows conservation of a putative Mdm2 
binding site. Representation of a portion of the human Mdm2 protein (ribbon 
structure) bound to the minimal binding site of the human p53 protein (stick 
structure). Modified from Kussie et al. (1996) using RasMol software. Four residues 
critical for hydrophobic binding interactions (F19, L22, "W23 and L26, shown in 
green) are conserved between human p53 and clam p53. Two potentially important 
changes are the substitutions in M. arenaria p53 of glutamic acid and tyrosine, for S20 
(red) and D21 (blue) in human p53, respectively (Holbrook et al., 2009). 
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Plate 4 Clam haemocyte cancer: (a) interference image of freshly collected 
cancerous clam haemocytes; notice dumb-bell shaped haemocyte involved in 
cytokinesis near the centre; (b) Cytospin of CCH stained with oil red O; notice 
mitotic division figure on the left of the image; (c) Quantum-dot 
immunocytochemical localization in cytoplasmically sequestered p53 in an 
individual CCH; (d) quantum-dot immunocytochemical localization in cytoplasmic 
mortalin in the same CCH; (e) dual quantum-dot immunocytochemical localization of 
cytoplasmically sequestered p53 and mortalin in the same CCH—yellow indicates co¬ 
localization of the two proteins (scale bar in (a) and (b) = 20 |Jim and in (c—e) 
= 10 (Jim) and (f) co-immunoprecipitation of clam p53 and mortalin in normal 
(NCH) and CCH. Std, protein standard; Cl and C2, first and second elute from 
CCH lysates; N1 and N2, first and second elude from NCH lysates; NA, first elude 
from gel loaded with NCH and CCH lysate mix (negative control); Q, first elude from 
quenched gel loaded with NCH and CCH lysate mix (negative control); Ccyt, 
cytoplasmic protein extract from CCH and Cnuc, nuclear protein extract from CCH. 
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Plate 5 Positive covariation between rates of gastrovascular transport and mitosis of 
stolon cells. (A) Normalized flow rate in two colonies (open or solid circles, respectively) 
was calculated as the maximum amplitude of lumen oscillations normalized to external 
stolon diameter, divided by the period (in seconds) of the oscillation cycle. (B) Mitotic 
index is calculated as the fraction of living cells (indicated by staining wtith DAPI) that are 
undergoing mitosis (indicated using the commercially available BrdU/anti-BrdU 
staining immunoassay). (C) Composite fluorescent image of a stolon showing staining 
by DAPI (blue; to estimate cell density in a region of stolon) and/or FITC (green; to 
identify mitotic cells). Image constructed from a Z-series stack of images from bottom to 
the top of the stolon lumen at 400 x magnification. Mitotic rates were estimated in two 
colonies (distinguished by solid vs. open symbols) along a linear length of stolon 
proximally near the polyp (circles) and distally near the stolon tip (triangles). Symbols 
and error bars represent means =b 95% confidence intervals of n = 4 colonies at each 
sampling period. 
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